Clare K. Fitzpatrick

Dynamic finite element knee simulation for evaluation of knee replacement mechanics

In vitro pre-clinical testing of total knee replacement (TKR) devices is a necessary step in the evaluation of new implant designs. Whole joint knee simulators, like the Kansas knee simulator (KKS), provide a controlled and repeatable loading environment for comparative evaluation of component designs or surgical alignment under dynamic conditions. Experimental testing, however, is time and cost prohibitive for design-phase evaluation of tens or hundreds of design variations. Experimentally-verified computational models provide an efficient platform for analysis of multiple components, sizes, and alignment conditions. The purpose of the current study was to develop and verify a computational model of a dynamic, whole joint knee simulator. Experimental internal-external and valgus-varus laxity tests, followed by dynamic deep knee bend and gait simulations in the KKS were performed on three cadaveric specimens. Specimen-specific finite element (FE) models of posterior-stabilized TKR were created from magnetic resonance images and CAD geometry. The laxity data was used to optimize mechanical properties of tibiofemoral soft-tissue structures on a specimen-specific basis. Each specimen was subsequently analyzed in a computational model of the experimental KKS, simulating both dynamic activities. The computational model represented all joints and actuators in the experimental setup, including a proportional-integral-derivative (PID) controller to drive quadriceps actuation. The computational model was verified against six degree-of-freedom patellofemoral (PF) and tibiofemoral (TF) kinematics and actuator loading during both deep knee bend and gait activities, with good agreement in trends and magnitudes between model predictions and experimental kinematics; differences were less than 1.8 mm and 2.2° for PF and TF translations and rotations. The whole joint FE simulator described in this study can be applied to investigate a wide range of clinical and research questions.

Development of a statistical shape model of the patellofemoral joint for investigating relationships between shape and function

Patellofemoral (PF)-related pathologies, including joint laxity, patellar maltracking, cartilage degradation and anterior knee pain, affect nearly 25% of the population. Researchers have investigated the influence of articular geometry on kinematics and contact mechanics in order to gain insight into the etiology of these conditions. The purpose of the current study was to create a three-dimensional statistical shape model of the PF joint and to characterize relationships between PF shape and function (kinematics and contact mechanics). A statistical shape model of the patellar and femoral articular surfaces and their relative alignment was developed from magnetic resonance images. Using 15 shape parameters, the model characterized 97% of the variation in the training set. The first three shape modes primarily described variation in size, patella alta-baja and depth of the sulcus groove. A previously verified finite element model was used to predict kinematics and contact mechanics for each subject. Combining the shape and joint mechanics data, a statistical shape-function model was developed that established quantitative relations of how changes in the shape of the PF joint influence mechanics. The predictive capability of the shape-function model was evaluated by comparing statistical model and finite element predictions, resulting in kinematic root mean square errors of less than 3° and 2.5 mm. The key results of the study are dually in the implementation of a novel approach linking statistical shape and finite element models and the relationships elucidated between PF articular geometry and mechanics.

The role of patient, surgical, and implant design variation in total knee replacement performance

Clinical studies demonstrate substantial variation in kinematic and functional performance within the total knee replacement (TKR) patient population. Some of this variation is due to differences in implant design, surgical technique and component alignment, while some is due to subject-specific differences in joint loading and anatomy that are inherently present within the population. Combined finite element and probabilistic methods were employed to assess the relative contributions of implant design, surgical, and subject-specific factors to overall tibiofemoral (TF) and patellofemoral (PF) joint mechanics, including kinematics, contact mechanics, joint loads, and ligament and quadriceps force during simulated squat, stance-phase gait and stepdown activities. The most influential design, surgical and subject-specific factors were femoral condyle sagittal plane radii, tibial insert superior-inferior (joint line) position and coronal plane alignment, and vertical hip load, respectively. Design factors were the primary contributors to condylar contact mechanics and TF anterior-posterior kinematics; TF ligament forces were dependent on surgical factors; and joint loads and quadriceps force were dependent on subject-specific factors. Understanding which design and surgical factors are most influential to TKR mechanics during activities of daily living, and how robust implant designs and surgical techniques must be in order to adequately accommodate subject-specific variation, will aid in directing design and surgical decisions towards optimal TKR mechanics for the population as a whole.

The influence of total knee arthroplasty geometry on mid-flexion stability: An experimental and finite element study

Fluoroscopic evaluation of total knee arthroplasty (TKA) has reported sudden anterior translation of the femur relative to the tibia (paradoxical anterior motion) for some cruciate-retaining designs. This motion may be tied to abrupt changes in the femoral sagittal radius of curvature characteristic of traditional TKA designs, as the geometry transitions from a large load-bearing distal radius to a smaller posterior radius which can accommodate femoral rollback. It was hypothesized that a gradually reducing radius may attenuate sudden changes in anterior-posterior motion that occur in mid-flexion with traditional discrete-radius designs. A combined experimental and computational approach was employed to test this hypothesis. A previously developed finite element (FE) model of the Kansas knee simulator (KKS), virtually implanted with multiple implant designs, was used to predict the amount of paradoxical anterior femoral slide during a simulated deep knee bend. The model predicted kinematics demonstrated that incorporating a gradually reducing radius in mid-flexion reduced the magnitude of paradoxical anterior translation between 21% and 68%, depending on the conformity of the tibial insert. Subsequently, both a dual-radius design and a modified design incorporating gradually reducing radii were tested in vitro in the KKS for verification. The model-predicted and experimentally observed kinematics exhibited good agreement, while the average experimental kinematics demonstrated an 81% reduction in anterior translation with the modified design. The FE model demonstrated sufficient sensitivity to appropriately differentiate kinematic changes due to subtle changes in implant design, and served as a useful pre-clinical design-phase tool to improve implant kinematics.

Relative contributions of design, alignment, and loading variability in knee replacement mechanics

Substantial variation in total knee replacement (TKR) outcomes exists within the patient population. Some of this variability is due to differences in the design of the implanted components and variation in surgical alignment, while other variability is due to differences in the applied forces and torques due to anatomic and physiological differences within a patient population. We evaluated the relative contributions of implant design, surgical alignment, and patient‐specific loading variability to overall tibiofemoral joint mechanics to provide insight into which measures can be influenced through design and surgical decisions, and which are inherently dependent on variation within the patient population and should be considered in the robustness of the implant design and surgical procedure. Design, surgical, and loading parameters were assessed using probabilistic finite element methods during simulated stance‐phase gait and squat activities. Patient‐specific loading was found to be the primary contributor to joint loading and kinematics during low flexion, particularly under conditions of high external loads (for instance, the gait cycle with high internal–external torque), while design and surgical factors, particularly femoral posterior radius and posterior slope of the tibial insert became increasingly important in TKR performance in deeper flexion.

Comparison of patellar bone strain in the natural and implanted knee during simulated deep flexion

Instances of anterior knee pain and patellar fracture are significant complications following total knee replacement (TKR). Bone strain measured in the patella can provide an indication of patellar fracture risk and may also be related to anterior knee pain. The objective of this study was to develop subject‐specific finite element models of the patellofemoral (PF) joint including density‐mapped material properties to gain insight into the patellar bone strain distribution in the natural and implanted knee. In eight subjects, the volume of bone experiencing strains >0.5% in the implanted condition was ∼200% larger, on average, than the natural condition. An inverse relationship with a correlation of −0.74 was established between postoperative bone volume and strain in the implanted specimens, suggesting that patellar geometry may be a useful indicator of postoperative strain. Comparing strains between regions (superior, inferior, medial, and lateral), it was found that although highly strained bone was evenly distributed between medial and lateral regions in the natural case, the implanted specimens demonstrated significantly larger volumes of highly strained bone medially as a result of substantially lower modulus bone in the medial compartment. Understanding distributions of PF strain may aid in preoperative identification of those patients at risk for patellar fracture or anterior knee pain, guidance regarding altered component placement for at‐risk patients, and design of components considering the implications of PF load transfer and patellar strain distribution.

Anatomic variation of the clavicle: A novel three‐dimensional study

An understanding of the complex anatomy of the clavicle is helpful in the treatment of clavicular fractures. Using three‐dimensional (3D) statistical shape analysis, the author presents a novel method to assess geometric morphology of the clavicle. Fifteen fresh frozen shoulder specimens were scanned using high‐resolution computerized tomography (CT) but four were excluded from the study. A further 16 high‐resolution CT scans of the clavicle were obtained by searching the hospital database. All 27 scans were reconstructed and subsequently imported into and analyzed using a specifically developed statistical software package. Using statistical shape analysis, geometric parameters were then measured. Both gender as well as side specific geometric morphology were observed. Clavicles in men were longer, wider, and thicker than in women. Right clavicles had a greater medial depth than left clavicles, especially in women. Clavicles in men had a greater lateral depth than in women. The sternal angle in women was larger than in men. Using 3D statistical shape analysis and applying it to the clavicle standardizes the study of its anatomy, rules out any variability, and calculates morphological parameters that are accurate, precise, and reproducible. This unique approach provides information that is useful not only to the clinician but also in the modification of current or design of future clavicle fixation devices. More importantly, from an anatomy standpoint, implementation of this novel approach in anatomical studies would eliminate intra‐ and interobserver variation and allow all studies to be standardized and thus more comparable. Clin. Anat. 23:199–209, 2010.

Computationally Efficient Finite Element Evaluation of Natural Patellofemoral Mechanics

Finite element methods have been applied to evaluate in vivo joint behavior, new devices, and surgical techniques but have typically been applied to a small or single subject cohort. Anatomic variability necessitates the use of many subject-specific models or probabilistic methods in order to adequately evaluate a device or procedure for a population. However, a fully deformable finite element model can be computationally expensive, prohibiting large multisubject or probabilistic analyses. The aim of this study was to develop a group of subject-specific models of the patellofemoral joint and evaluate trade-offs in analysis time and accuracy with fully deformable and rigid body articular cartilage representations. Finite element models of eight subjects were used to tune a pressure-overclosure relationship during a simulated deep flexion cycle. Patellofemoral kinematics and contact mechanics were evaluated and compared between a fully deformable and a rigid body analysis. Additional eight subjects were used to determine the validity of the rigid body pressure-overclosure relationship as a subject-independent parameter. There was good agreement in predicted kinematics and contact mechanics between deformable and rigid analyses for both the tuned and test groups. Root mean square differences in kinematics were less than 0.5 deg and 0.2 mm for both groups throughout flexion. Differences in contact area and peak and average contact pressures averaged 5.4%, 9.6%, and 3.8%, respectively, for the tuned group and 6.9%, 13.1%, and 6.4%, respectively, for the test group, with no significant differences between the two groups. There was a 95% reduction in computational time with the rigid body analysis as compared with the deformable analysis. The tuned pressure-overclosure relationship derived from the patellofemoral analysis was also applied to tibiofemoral (TF) articular cartilage in a group of eight subjects. Differences in contact area and peak and average contact pressures averaged 8.3%, 11.2%, and 5.7% between rigid and deformable analyses in the tibiofemoral joint. As statistical, probabilistic, and optimization techniques can require hundreds to thousands of analyses, a viable platform is crucial to component evaluation or clinical applications. The computationally efficient rigid body platform described in this study may be integrated with statistical and probabilistic methods and has potential clinical application in understanding in vivo joint mechanics on a subject-specific or population basis.

Computational analysis of factors contributing to patellar dislocation.

Treatment for an initial incidence of patellar dislocation is usually conservative management; however, almost half of patients experience a subsequent, or multiple, dislocation(s). Patients often undergo multiple procedures which fail to treat the underlying anatomic abnormalities. The objective of this study was to evaluate interactions between key predisposing anatomic factors to patellar dislocation and identify combinations of abnormal factors which increase the risk of recurrent lateral dislocation. Four factors associated with lateral patellar dislocation were identified (sulcus angle, Insall‐Salvati ratio, tibial tubercle‐trochlear groove distance, and femoral anteversion). A finite element model of the patellofemoral joint was developed and parameterized so that a value for each factor could be applied and the model geometry/alignment would be modified accordingly. 100 combinations of the four factors were generated in separate computational simulations and resulting kinematics and forces of the patellofemoral joint were recorded. Sulcus angle was the most impactful factor on constraint. Multiple abnormal factors were generally required to produce the extremes of patellar alignment observed in this analysis. Understanding the underlying anatomic factors, and their effect on joint mechanics, for patients with recurrent lateral patellar dislocation will aid in determining optimal treatment pathways on a patient‐specific basis.

Identifying alignment parameters affecting implanted patellofemoral mechanics

Complications of the patellofemoral (PF) joint remain a common cause for revision of total knee replacements. PF complications, such as patellar maltracking, subluxation, and implant failure, have been linked to femoral and patellar component alignment. In this study, a dynamic finite element model of an implanted PF joint was applied in conjunction with a probabilistic simulation to establish relationships between alignment parameters and PF kinematics, contact mechanics, and internal stresses. Both traditional sensitivity analysis and a coupled probabilistic and principal component analysis approach were applied to characterize relationships between implant alignment and resulting joint mechanics. Critical alignment parameters, and combinations of parameters, affecting PF mechanics were identified for three patellar designs (dome, modified dome, and anatomic). Femoral internal–external (I‐E) alignment was identified as a critical alignment factor for all component designs, influencing medial–lateral contact force and anterior–posterior translation. The anatomic design was sensitive to patellar flexion–extension (F‐E) alignment, while the dome, as expected, was less influenced by rotational alignment, and more by translational position. The modified dome was sensitive to a combination of superior–inferior, F‐E, and I‐E alignments. Understanding the relationships and design‐specific dependencies between alignment parameters can aid preoperative planning, and help focus instrumentation design on those alignment parameters of primary concern.