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.

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.

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.

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.

Comparison of explicit finite element and mechanical simulation of the proximal femur during dynamic drop-tower testing.

Current screening techniques based on areal bone mineral density (aBMD) measurements are unable to identify the majority of people who sustain hip fractures. Biomechanical examination of such events may help determine what predisposes a hip to be susceptible to fracture. Recently, drop-tower simulations of in-vitro sideways falls have allowed the study of the mechanical response of the proximal human femur at realistic impact speeds. This technique has created an opportunity to validate explicit finite element (FE) models against dynamic test data. This study compared the outcomes of 15 human femoral specimens fractured using a drop tower with complementary specimen-specific explicit FE analysis. Correlation coefficient and root mean square error (RMSE) were found to be moderate for whole bone stiffness comparison (R(2)=0.3476 and 22.85% respectively). No correlation was found between experimentally and computationally predicted peak force, however, energy absorption comparison produced moderate correlation and RMSE (R(2)=0.4781 and 29.14% respectively). By comparing predicted strain maps to high speed video data we demonstrated the ability of the FE models to detect vulnerable portions of the bones. Based on our observations, we conclude that there exists a need to extend the current apparent level material models for bone to cover higher strain rates than previously tested experimentally.

Development of a balanced experimental-computational approach to understanding the mechanics of proximal femur fractures.

The majority of people who sustain hip fractures after a fall to the side would not have been identified using current screening techniques such as areal bone mineral density. Identifying them, however, is essential so that appropriate pharmacological or lifestyle interventions can be implemented. A protocol, demonstrated on a single specimen, is introduced, comprising the following components; in vitro biofidelic drop tower testing of a proximal femur; high-speed image analysis through digital image correlation; detailed accounting of the energy present during the drop tower test; organ level finite element simulations of the drop tower test; micro level finite element simulations of critical volumes of interest in the trabecular bone. Fracture in the femoral specimen initiated in the superior part of the neck. Measured fracture load was 3760N, compared to 4871N predicted based on the finite element analysis. Digital image correlation showed compressive surface strains as high as 7.1% prior to fracture. Voxel level results were consistent with high-speed video data and helped identify hidden local structural weaknesses. We found using a drop tower test protocol that a femoral neck fracture can be created with a fall velocity and energy representative of a sideways fall from standing. Additionally, we found that the nested explicit finite element method used allowed us to identify local structural weaknesses associated with femur fracture initiation.

Influence of cement stiffness and bone morphology on the compressive properties of bone-cement composites in simulated vertebroplasty.

Vertebroplasty is widely used to treat vertebral compression fractures. Little is known about the influence of morphological parameters of the bone on the properties of the bone-cement composite. Furthermore, although generic finite element (FE) models have been suggested as a way to compute the values of these properties, their accuracy has not been established. In the experimental part of this study, we tested bovine cancellous bone and three different polymethylmethacrylate bone cements and determined six quasi-static uniaxial compressive properties of bone-cement composite specimens and 10 morphological parameters of the bone. For the FE work, we used two simulations, one being µFE and the other unit cell FE. In conclusion, we found that (1) for composite specimens, that relative contribution of the cement to the overall response of the composite increases with increasing cement stiffness; (2) the anisotropy ratio is the bone morphological property that exerts the most significant influence on the experimentally obtained compressive properties of the bone-cement composites determined; (3) the accuracy of the computed compressive properties of the composites ranged from low to high, depending on simulation method used. The largest errors, however, can partially be explained by difference in boundary conditions between the experimental testing and the simulation techniques. The lattermost finding points to the potential for simplified FE models being incorporated into automatic material mapping schemes in whole bone vertebra FE simulations.

Experimental and computational models to investigate the effect of adhesion on the mechanical properties of bone-cement composites.

A generic finite element approach was developed to study the effect of adhesion on the mechanical response of bone cement composites and validated against literature data. The results showed that a zero friction bone-cement (PMMA) interface conditions captured the results of the experimental testing better than assuming a fully bonded interface. An experimental model for studying the effect of interface adhesion in a bone-cement like composite was also developed in the present study. The results using this model indicate that the difference in Young's modulus and ultimate strength between a fully bonded interface and unbonded interface is approximately 30% for bone volume fraction similar to what can be found in osteoporotic vertebrae. Apart from concluding that bone to cement adhesion is a major contributor to the mechanical response of bone-cement composites, our studies based on the generic FE approach also indicate that the mechanical properties of the cement is the most important contributor to the resulting mechanical properties of the composite at bone volume fraction relevant in terms of vertebroplasty treatment.