Biofidelic finite element models for accurately classifying hip fracture in a retrospective clinical study of elderly women from the AGES Reykjavik cohort.

Clinical retrospective studies have only reported limited improvements in hip fracture classification accuracy using finite element (FE) models compared to conventional areal bone mineral density (aBMD) measurements. A possible explanation is that state-of-the-art quasi-static models do not estimate patient-specific loads. A novel FE modeling technique was developed to improve the biofidelity of simulated impact loading from sideways falling. This included surrogate models of the pelvis, lower extremities, and soft tissue that were morphed based on subject anthropometrics. Hip fracture prediction models based on aBMD and FE measurements were compared in a retrospective study of 254 elderly female subjects from the AGES-Reykjavik study. Subject fragility ratio (FR) was defined as the ratio between the ultimate forces of paired biofidelic models, one with linear elastic and the other with non-linear stress-strain relationships in the proximal femur. The expected end-point value (EEV) was defined as the FR weighted by the probability of one sideways fall over five years, based on self-reported fall frequency at baseline. The change in maximum volumetric strain (ΔMVS) on the surface of the femoral neck was calculated between time of ultimate femur force and 90% post-ultimate force in order to assess the extent of tensile tissue damage present in non-linear models. After age-adjusted logistic regression, the area under the receiver-operator curve (AUC) was highest for ΔMVS (0.72), followed by FR (0.71), aBMD (0.70), and EEV (0.67), however the differences between FEA and aBMD based prediction models were not deemed statistically significant. When subjects with no history of falling were excluded from the analysis, thus artificially assuming that falls were known a priori with no uncertainty, a statistically significant difference in AUC was detected between ΔMVS (0.85), and aBMD (0.74). Multivariable linear regression suggested that the variance in maximum elastic femur force was best explained by femoral head radius, pelvis width, and soft tissue thickness (R2 = 0.79; RMSE = 0.46 kN; p < 0.005). Weighting the hip fracture prediction models based on self-reported fall frequency did not improve the models' sensitivity, however excluding non-fallers lead to significant differences between aBMD and FE based models. These findings suggest that an accurate assessment of fall probability is necessary for accurately identifying individuals predisposed to hip fracture.

Larger vertebral endplate concavities cause higher failure load and work at failure under high-rate impact loading of rabbit spinal explants.

Vertebral fractures are among the most common of all osteoporosis related fracture types and its risk assessment is largely based on bone quality measures. Morphometric parameters are not yet considered, although endplate thickness and concavity shape were found to be important in fracture prediction in low-rate tests. We hypothesized that, under high-rate impact loading, the shape and size of the central endplate concavity are of key importance for fracture prediction. Therefore, we tested rabbit spinal segment explants in vitro under high-rate impact loading. With a combination of microCT to describe endplate morphometry, high-speed video imaging, and impact force measurement, endplate morphometry was correlated to the mechanical response. We found that endplate concavity shape and volume were important in describing the mechanical response: larger concavities caused higher failure load. We suggest a model for the fracture mechanism under high-rate impact loading, considering the morphometry of the endplates: wider and more voluminous concavities are protective whereas steeper slopes of the concavity edges and increasing bone volume fraction of the central endplate moiety are disadvantageous. Therefore, the shape and size of endplate morphometry are important in vertebral fracture prediction and should be considered included in vertebral fracture risk assessment.

On the Failure Initiation in the Proximal Human Femur Under Simulated Sideways Fall.

The limitations of areal bone mineral density measurements for identifying at-risk individuals have led to the development of alternative screening methods for hip fracture risk including the use of geometrical measurements from the proximal femur and subject specific finite element analysis (FEA) for predicting femoral strength, based on quantitative CT data (qCT). However, these methods need more development to gain widespread clinical applications. This study had three aims: To investigate whether proximal femur geometrical parameters correlate with obtained femur peak force during the impact testing; to examine whether or not failure of the proximal femur initiates in the cancellous (trabecular) bone; and finally, to examine whether or not surface fracture initiates in the places where holes perforate the cortex of the proximal femur. We found that cortical thickness around the trochanteric-fossa is significantly correlated to the peak force obtained from simulated sideways falling (R 2 = 0.69) more so than femoral neck cortical thickness (R 2 = 0.15). Dynamic macro level FE simulations predicted that fracture generally initiates in the cancellous bone compartments. Moreover, our micro level FEA results indicated that surface holes may be involved in primary failure events.

Material mapping strategy to improve the predicted response of the proximal femur to a sideways fall impact.

Sideways falls are largely responsible for the highly prevalent osteoporotic hip fractures in today's society. These injuries are dynamic events, therefore dynamic FE models validated with dynamic ex vivo experiments provide a more realistic simulation than simple quasi-static analysis. Drop tower experiments using cadaveric specimens were used to identify the material mapping strategy that provided the most realistic mechanical response under impact loading. The present study tested the addition of compression-tension asymmetry, tensile bone damage, and cortical-specific strain rate dependency to the material mapping strategy of fifteen dynamic FE models of the proximal femur, and found improved correlations and reduced error for whole bone stiffness (R2 = 0.54, RSME = 0.87kN/mm) and absolute maximum force (R2 = 0.56, RSME =0.57kN), and a high correlation in impulse response (R2 = 0.82, RSME =12.38kg/s). Simulations using fully bonded nodes between the rigid bottom plate and PMMA cap supporting the femoral head had higher correlations and less error than simulations using a frictionless sliding at this contact surface. Strain rates over 100/s were observed in certain elements in the femoral neck and trochanter, indicating that additional research is required to better quantify the strain rate dependencies of both trabecular and cortical bone at these strain rates. These results represent the current benchmark in dynamic FE modeling of the proximal femur in sideways falls. Future work should also investigate improvements in experimental validation techniques by developing better displacement measurements and by enhancing the biofidelity of the impact loading wherever possible.

Morphology based anisotropic finite element models of the proximal femur validated with experimental data.

Finite element analysis (FEA) of bones scanned with Quantitative Computed Tomography (QCT) can improve early detection of osteoporosis. The accuracy of these models partially depends on the assigned material properties, but anisotropy of the trabecular bone cannot be fully captured due to insufficient resolution of QCT. The inclusion of anisotropy measured from high resolution peripheral QCT (HR-pQCT) could potentially improve QCT-based FEA of the femur, although no improvements have yet been demonstrated in previous experimental studies. This study analyzed the effects of adding anisotropy to clinical resolution femur models by constructing six sets of FE models (two isotropic and four anisotropic) for each specimen from a set of sixteen femurs that were experimentally tested in sideways fall loading with a strain gauge on the superior femoral neck. Two different modulus-density relationships were tested, both with and without anisotropy derived from mean intercept length analysis of HR-pQCT scans. Comparing iso- and anisotropic models to the experimental data resulted in nearly identical correlation and highly similar linear regressions for both whole bone stiffness and strain gauge measurements. Anisotropic models contained consistently greater principal compressive strains, approximately 14% in magnitude, in certain internal elements located in the femoral neck, greater trochanter, and femoral head. In summary, anisotropy had minimal impact on macroscopic measurements, but did alter internal strain behavior. This suggests that organ level QCT-based FE models measuring femoral stiffness have little to gain from the addition of anisotropy, but studies considering failure of internal structures should consider including anisotropy to their models.

The influence of the modulus-density relationship and the material mapping method on the simulated mechanical response of the proximal femur in side-ways fall loading configuration.

Contributing to slow advance of finite element (FE) simulations for hip fracture risk prediction, into clinical practice, could be a lack of consensus in the biomechanics community on how to map properties to the models. Thus, the aim of the present study was first, to systematically quantify the influence of the modulus-density relationship (E-ρ) and the material mapping method (MMM) on the predicted mechanical response of the proximal femur in a side-ways fall (SWF) loading configuration and second, to perform a model-to-model comparison of the predicted mechanical response within the femoral neck for all the specimens tested in the present study, using three different modelling techniques that have yielded good validation outcome in terms of surface strain prediction and whole bone response according to the literature. We found the outcome to be highly dependent on both the E-ρ relationship and the MMM. In addition, we found that the three modelling techniques that have resulted in good validation outcome in the literature yielded different principal strain prediction both on the surface as well as internally in the femoral neck region of the specimens modelled in the present study. We conclude that there exists a need to carry out a more comprehensive validation study for the SWF loading mode to identify which combination of MMMs and E-ρ relationship leads to the best match for whole bone and local mechanical response. The MMMs tested in the present study have been made publicly available at https://simtk.org/home/mitk-gem.