Azhar A. Ali

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.

Specimen-specific modeling of hip fracture pattern and repair

Hip fracture remains a major health problem for the elderly. Clinical studies have assessed fracture risk based on bone quality in the aging population and cadaveric testing has quantified bone strength and fracture loads. Prior modeling has primarily focused on quantifying the strain distribution in bone as an indicator of fracture risk. Recent advances in the extended finite element method (XFEM) enable prediction of the initiation and propagation of cracks without requiring a priori knowledge of the crack path. Accordingly, the objectives of this study were to predict femoral fracture in specimen-specific models using the XFEM approach, to perform one-to-one comparisons of predicted and in vitro fracture patterns, and to develop a framework to assess the mechanics and load transfer in the fractured femur when it is repaired with an osteosynthesis implant. Five specimen-specific femur models were developed from in vitro experiments under a simulated stance loading condition. Predicted fracture patterns closely matched the in vitro patterns; however, predictions of fracture load differed by approximately 50% due to sensitivity to local material properties. Specimen-specific intertrochanteric fractures were induced by subjecting the femur models to a sideways fall and repaired with a contemporary implant. Under a post-surgical stance loading, model-predicted load sharing between the implant and bone across the fracture surface varied from 59%:41% to 89%:11%, underscoring the importance of considering anatomic and fracture variability in the evaluation of implants. XFEM modeling shows potential as a macro-level analysis enabling fracture investigations of clinical cohorts, including at-risk groups, and the design of robust implants.

Effects of resection thickness on mechanics of resurfaced patellae.

Patellar resection thickness during total knee replacement (TKR) has been cited as a contributor to patellar fracture, anterior knee pain and quadriceps efficiency; however, optimal thickness required to minimize clinical complications remains unclear. The objectives of the current study were to determine how patellar resection thickness and bone quality impacts patellar bone strain, kinematics, and quadriceps efficiency. A series of specimen-specific finite element models of the knee joint with distributed patellar bone material properties were developed. Each specimen was virtually implanted with a TKR system. Each specimen was analyzed with patellar bone resected to thicknesses which varied from 9 to 14 mm. Simulations with reduced modulus bone were also performed. Each model perturbation was evaluated during a dynamic squat cycle, and bone strain, quadriceps force and six-degree-of-freedom kinematics were predicted. Highest peak bone strain was predicted in the thinnest patellae, indicating greatest risk of patellar fracture; highest median bone strain was predicted in the thickest patellae. Consistent differences in quadriceps efficiency were predicted; in early flexion the thickest patellae required the least quadriceps force. Greater sagittal plane tilt was observed for the thinnest patellae. Reduced modulus models (50% lower modulus) demonstrated an increase in peak bone strain of up to seven times the original modulus models. Understanding the complex interactions between patellar resection thickness, muscle requirements, kinematics, bone quality, and bone property distribution may aid in developing an understanding of which patients are most at risk from patellar fracture and anterior knee pain and how best to treat individuals to reduce potential complications.

Validation of predicted patellofemoral mechanics in a finite element model of the healthy and cruciate-deficient knee.

Healthy patellofemoral (PF) joint mechanics are critical to optimal function of the knee joint. Patellar maltracking may lead to large joint reaction loads and high stresses on the articular cartilage, increasing the risk of cartilage wear and the onset of osteoarthritis. While the mechanical sources of PF joint dysfunction are not well understood, links have been established between PF tracking and abnormal kinematics of the tibiofemoral (TF) joint, specifically following cruciate ligament injury and repair. The objective of this study was to create a validated finite element (FE) representation of the PF joint in order to predict PF kinematics and quadriceps force across healthy and pathological specimens. Measurements from a series of dynamic in-vitro cadaveric experiments were used to develop finite element models of the knee for three specimens. Specimens were loaded under intact, ACL-resected and both ACL and PCL-resected conditions. Finite element models of each specimen were constructed and calibrated to the outputs of the intact knee condition, and subsequently used to predict PF kinematics, contact mechanics, quadriceps force, patellar tendon moment arm and patellar tendon angle of the cruciate resected conditions. Model results for the intact and cruciate resected trials successfully matched experimental kinematics (avg. RMSE 4.0°, 3.1mm) and peak quadriceps forces (avg. difference 5.6%). Cruciate resections demonstrated either increased patellar tendon loads or increased joint reaction forces. The current study advances the standard for evaluation of PF mechanics through direct validation of cruciate-resected conditions including specimen-specific representations of PF anatomy.

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 5mm anterior tibial translation) while increased TF joint loading was shown in the other (up to 400N).

A Combined Experimental and Computational Approach to Subject-Specific Analysis of Knee Joint Laxity

Modeling complex knee biomechanics is a continual challenge, which has resulted in many models of varying levels of quality, complexity, and validation. Beyond modeling healthy knees, accurately mimicking pathologic knee mechanics, such as after cruciate rupture or meniscectomy, is difficult. Experimental tests of knee laxity can provide important information about ligament engagement and overall contributions to knee stability for development of subject-specific models to accurately simulate knee motion and loading. Our objective was to provide combined experimental tests and finite-element (FE) models of natural knee laxity that are subject-specific, have one-to-one experiment to model calibration, simulate ligament engagement in agreement with literature, and are adaptable for a variety of biomechanical investigations (e.g., cartilage contact, ligament strain, in vivo kinematics). Calibration involved perturbing ligament stiffness, initial ligament strain, and attachment location until model-predicted kinematics and ligament engagement matched experimental reports. Errors between model-predicted and experimental kinematics averaged <2 deg during varus-valgus (VV) rotations, <6 deg during internal-external (IE) rotations, and <3 mm of translation during anterior-posterior (AP) displacements. Engagement of the individual ligaments agreed with literature descriptions. These results demonstrate the ability of our constraint models to be customized for multiple individuals and simultaneously call attention to the need to verify that ligament engagement is in good general agreement with literature. To facilitate further investigations of subject-specific or population based knee joint biomechanics, data collected during the experimental and modeling phases of this study are available for download by the research community.

In vivo comparison of medialized dome and anatomic patellofemoral geometries using subject-specific computational modeling: SUBJECT-SPECIFIC COMPUTATIONAL MODELING

Successful outcome following total knee arthroplasty (TKA) with patella resurfacing is partly determined by the restoration of patellofemoral (PF) function and recovery of the quadriceps mechanism. The current study compared two patellar TKA geometries (medialized dome and anatomic) to determine their impact on PF mechanics and quadriceps function. In‐vivo, subject‐specific patellar mechanics were evaluated using a sequential experimental and modeling approach. First, stereo radiography, marker‐based motion capture, and force plate data were collected for TKA patients (10 dome, 10 anatomic) performing a knee extension and lunge. Second, subject‐specific, whole‐body, musculoskeletal models, including 6 degrees‐of‐freedom (DOF) knee joint kinematics, were created for each subject and activity to predict quadriceps forces. Last, finite element models of each subject and activity were created to predict PF kinematics, patellar loading, moment arm, and patellar tendon angle. Differences in mechanics between dome and anatomic patients were highlighted during load‐bearing (lunge) activity. Anatomic subjects demonstrated greater PF flexion angles (avg. 11 ± 3°) compared to dome subjects during lunge. Similar to the natural knee, contact locations on the patella migrated inferior to superior as the knee flexed in anatomic subjects, but remained relatively superior in dome subjects. Differences in kinematics and contact location likely contributed to altered mechanics with anatomic subjects presenting greater load transfer from the quadriceps to the patellar tendon in deep flexion (>75°), and dome subjects demonstrating larger contact forces during lunge. Although there was substantial patient variability, evaluations of PF mechanics suggested improved quadriceps function and more natural kinematics in the anatomic design.

The interaction of muscle moment arm, knee laxity, and torque in a multi-scale musculoskeletal model of the lower limb

INTRODUCTION Musculoskeletal modeling allows insight into the interaction of muscle force and knee joint kinematics that cannot be measured in the laboratory. However, musculoskeletal models of the lower extremity commonly use simplified representations of the knee that may limit analyses of the interaction between muscle forces and joint kinematics. The goal of this research was to demonstrate how muscle forces alter knee kinematics and consequently muscle moment arms and joint torque in a musculoskeletal model of the lower limb that includes a deformable representation of the knee. METHODS Two musculoskeletal models of the lower limb including specimen-specific articular geometries and ligament deformability at the knee were built in a finite element framework and calibrated to match mean isometric torque data collected from 12 healthy subjects. Muscle moment arms were compared between simulations of passive knee flexion and maximum isometric knee extension and flexion. In addition, isometric torque results were compared with predictions using simplified knee models in which the deformability of the knee was removed and the kinematics at the joint were prescribed for all degrees of freedom. RESULTS Peak isometric torque estimated with a deformable knee representation occurred between 45° and 60° in extension, and 45° in flexion. The maximum isometric flexion torques generated by the models with deformable ligaments were 14.6% and 17.9% larger than those generated by the models with prescribed kinematics; by contrast, the maximum isometric extension torques generated by the models were similar. The change in hamstrings moment arms during isometric flexion was greater than that of the quadriceps during isometric extension (a mean RMS difference of 9.8 mm compared to 2.9 mm, respectively). DISCUSSION The large changes in the moment arms of the hamstrings, when activated in a model with deformable ligaments, resulted in changes to flexion torque. When simulating human motion, the inclusion of a deformable joint in a multi-scale musculoskeletal finite element model of the lower limb may preserve the realistic interaction of muscle force with knee kinematics and torque.