Identification of elastic properties of human patellae using micro-finite element analysis

Patellar bone strain after total knee arthroplasty is correlated with bone mineral density and body mass index

Patella-related complications after total knee arthroplasty (TKA) remain a major clinical concern. Previous studies have suggested that increased postoperative patellar bone strain could be related to such complications, but there is limited knowledge on patellar strain after TKA. The objective of this study was to predict patellar bone strain after TKA and evaluate correlations with various preoperative data. Fourteen TKA patients with a minimum follow-up of one year were included in this study. Using preoperative CT datasets, preoperative planning, and postoperative X-rays, a method is presented to generate patient-specific finite element models after virtual TKA. Patellar kinematics and forces were predicted during a squat movement, and patellar bone strain was evaluated at 60° of knee flexion. Strain varied greatly among patients, but was strongly negatively correlated (r = -0.85, p < 0.001) with bone mineral density (BMD) and moderately positively (r  = 0.54, p  =  0.05) with body mass index (BMI). The BMI/BMD ratio explained 87% of strain, and should be further investigated as a potential risk factor for clinical complications. This study represents a preliminary step towards the identification of patients at risk of patellar complications after TKA.

A comparison of density-modulus relationships used in finite element modeling of the shoulder

Subject- and site-specific modeling techniques greatly improve the accuracy of computational models derived from clinical-resolution quantitative computed tomography (QCT) data. The majority of shoulder finite element (FE) studies use density-modulus relationships developed for alternative anatomical locations. As such, the objectives of this study were to compare the six most commonly used density-modulus relationships in shoulder finite element (FE) studies. To achieve this, ninety-eight (98) virtual trabecular bone cores were extracted from uCT scans of scapulae from 14 cadaveric specimens (7 male; 7 female). Homogeneous tissue moduli of 20 GPa, and heterogeneous tissue moduli scaled by CT-intensity were considered. Micro finite element models (µ-FEMs) of each virtual core were compressively loaded to 0.5% apparent strain and apparent strain energy density (SED_app) was collected. Each uCT virtual core was then co-registered to clinical QCT images, QCT-FEMs created, and each of the 6 density-modulus relationships applied (6 × 98 = 588 QCT-FEMs). The loading and boundary conditions were replicated and SED_app was collected and compared to µ-FEM SED_app. When a homogeneous tissue modulus was considered in the µ-FEMs, SED_app was best predicted in QCT-FEMs with the density-modulus relationship developed from pooled anatomical locations (QCT-FEM SED_app = 0.979µ-FEM SED_app + 0.0066, r² = 0.933). A different density-modulus relationship best predicted SED_app (QCT-FEM SED_app = 1.014µ-FEM SED_app + 0.0034, r² = 0.935) when a heterogeneous tissue modulus was considered. This study compared density-modulus relationships used in shoulder FE studies using an independent computational methodology for comparing these relationships.

μCT-based trabecular anisotropy can be reproducibly computed from HR-pQCT scans using the triangulated bone surface

The trabecular structure can be assessed at the wrist or tibia via high-resolution peripheral quantitative computed tomography (HR-pQCT). Yet on this modality, the performance of the existing methods, evaluating trabecular anisotropy is usually overlooked, especially in terms of reproducibility. We thus proposed to compare the TRI routine used by SCANCO Medical AG (Brüttisellen, Switzerland), the classical mean intercept length (MIL), and the grey-level structure tensor (GST) to the mean surface length (MSL), a new method for evaluating a second-order fabric tensor based on the triangulation of the bone surface. The distal radius of 24 fresh-frozen human forearms was scanned three times via HR-pQCT protocols (61μm, 82μm nominal voxel size), dissected, and imaged via micro computed tomography (μCT) at 16μm nominal voxel size. After registering the scans, we compared for each resolution the fabric tensors, determined by the mentioned techniques for 182 trabecular regions of interest. We then evaluated the reproducibility of the fabric information measured by HR-pQCT via precision errors. On μCT, TRI and GST were respectively the best and worst surrogates for MIL_μCT (MIL computed on μCT) in terms of eigenvalues and main direction of anisotropy. On HR-pQCT, however, MSL provided the best approximation of MIL_μCT. Surprisingly, surface-based approaches (TRI, MSL) also proved to be more precise than both MIL and GST. Our findings confirm that MSL can reproducibly estimate MIL_μCT, the current gold standard. MSL thus enables the direct mapping of the fabric-dependent material properties required in homogenised HR-pQCT-based finite element models.