Finite element derived femoral strength is a better predictor of hip fracture risk than aBMD in the AGES Reykjavik study cohort.

Hip fractures associated with a high economic burden, loss of independence, and a high rate of post-fracture mortality, are a major health concern for modern societies. Areal bone mineral density is the current clinical metric of choice when assessing an individual's future risk of fracture. However, this metric has been shown to lack sensitivity and specificity in the targeted selection of individuals for preventive interventions. Although femoral strength derived from computed tomography based finite element models has been proposed as an alternative based on its superior femoral strength prediction ex vivo, such predictions have only shown marginal or no improvement for assessing hip fracture risk. This study compares finite element derived femoral strength to aBMD as a metric for hip fracture risk assessment in subjects (N = 601) from the AGES Reykjavik Study cohort and analyses the dependence of femoral strength predictions and classification accuracy on the material model and femoral loading alignment. We found hip fracture classification based on finite element derived femoral strength to be significantly improved compared to aBMD. Finite element models with non-linear material models performed better at classifying hip fractures compared to finite element models with linear material models and loading alignments with low internal rotation and adduction, which do not correspond to weak femur alignments, were found to be most suitable for hip fracture classification.

Correction: On the internal reaction forces, energy absorption, and fracture in the hip during simulated sideways fall impact.

[This corrects the article DOI: 10.1371/journal.pone.0200952.].

On the internal reaction forces, energy absorption, and fracture in the hip during simulated sideways fall impact.

The majority of hip fractures have been reported to occur as a result of a fall with impact to the greater trochanter of the femur. Recently, we developed a novel cadaveric pendulum-based hip impact model and tested two cadaveric femur-pelvis constructs, embedded in a soft tissue surrogate. The outcome was a femoral neck fracture in a male specimen while a female specimen had no fracture. The aim of the present study was, first, to develop a methodology for constructing and assessing the accuracy of explicit Finite Element Models (FEMs) for simulation of sideways falls to the hip based on the experimental model. Second, to use the FEMs for quantifying the internal reaction forces and energy absorption in the hip during impact. Third, to assess the potential of the FEMs in terms of separating a femoral fracture endpoint from a non-fracture endpoint. Using a non-linear, strain rate dependent, and heterogeneous material mapping strategy for bone tissue in these models, we found the FEM-derived results to closely match the experimental test results in terms of impact forces and displacements of pelvic video markers up to the time of peak impact force with errors below 10%. We found the internal reaction forces in the femoral neck on the impact side to be approximately 35% lower than the impact force measured between soft tissue and ground for both specimens. In addition, we found the soft tissue to be the component that absorbed the largest part of the energy of the tissue types in the hip region. Finally, we found surface strain patterns derived from FEM results to match the fracture location and extent based on post testing x-rays of the specimens. This is the first study with quantitative data on the energy absorption in the pelvic region during a sideways fall.

Strain rate dependency of bovine trabecular bone under impact loading at sideways fall velocity.

There is currently a knowledge gap in scientific literature concerning the strain rate dependent properties of trabecular bone at intermediate strain rates. Meanwhile, strain rates between 10 and 200/s have been observed in previous dynamic finite element models of the proximal femur loaded at realistic sideways fall speeds. This study aimed to quantify the effect of strain rate (ε̇) on modulus of elasticity (E), ultimate stress (σu), failure energy (Uf), and minimum stress (σm) of trabecular bone in order to improve the biofidelity of material properties used in dynamic simulations of sideways fall loading on the hip. Cylindrical cores of trabecular bone (D = 8 mm, Lgauge = 16 mm, n = 34) from bovine proximal tibiae and distal femurs were scanned in µCT (10 µm), quantifying apparent density (ρapp) and degree of anisotropy (DA), and subsequently impacted within a miniature drop tower. Force of impact was measured using a piezoelectric load cell (400 kHz), while displacement during compression was measured from high speed video (50,000 frames/s). Four groups, with similar density distributions, were loaded at different impact velocities (0.84, 1.33, 1.75, and 2.16 m/s) with constant kinetic energy (0.4 J) by adjusting the impact mass. The mean strain rates of each group were significantly different (p < 0.05) except for the two fastest impact speeds (p = 0.09). Non-linear regression models correlated strain rate, DA, and ρapp with ultimate stress (R2 = 0.76), elastic modulus (R2 = 0.63), failure energy (R2 = 0.38), and minimum stress (R2 = 0.57). These results indicate that previous estimates of σu could be under predicting the mechanical properties at strain rates above 10/s.

Interactive graph-cut segmentation for fast creation of finite element models from clinical ct data for hip fracture prediction.

In this study, we propose interactive graph cut image segmentation for fast creation of femur finite element (FE) models from clinical computed tomography scans for hip fracture prediction. Using a sample of N = 48 bone scans representing normal, osteopenic and osteoporotic subjects, the proximal femur was segmented using manual (gold standard) and graph cut segmentation. Segmentations were subsequently used to generate FE models to calculate overall stiffness and peak force in a sideways fall simulations. Results show that, comparable FE results can be obtained with the graph cut method, with a reduction from 20 to 2-5 min interaction time. Average differences between segmentation methods of 0.22 mm were not significantly correlated with differences in FE derived stiffness (R2 = 0.08, p = 0.05) and weakly correlated to differences in FE derived peak force (R2 = 0.16, p = 0.01). We further found that changes in automatically assigned boundary conditions as a consequence of small segmentation differences were significantly correlated with FE derived results. The proposed interactive graph cut segmentation software MITK-GEM is freely available online at https://simtk.org/home/mitk-gem .

Mapping anisotropy of the proximal femur for enhanced image based finite element analysis.

Finite element (FE) models of bone derived from quantitative computed tomography (QCT) rely on realistic material properties to accurately predict bone strength. QCT cannot resolve bone microarchitecture, therefore QCT-based FE models lack the anisotropy apparent within the underlying bone tissue. This study proposes a method for mapping femoral anisotropy using high-resolution peripheral quantitative computed tomography (HR-pQCT) scans of human cadaver specimens. Femur HR-pQCT images were sub-divided into numerous overlapping cubic sub-volumes and the local anisotropy was quantified using a 'direct-mechanics' method. The resulting directionality reflected all the major stress lines visible within the trabecular lattice, and provided a realistic estimate of the alignment of Harvesian systems within the cortical compartment. QCT-based FE models of the proximal femur were constructed with isotropic and anisotropic material properties, with directionality interpolated from the map of anisotropy. Models were loaded in a sideways fall configuration and the resulting whole bone stiffness was compared to experimental stiffness and ultimate strength. Anisotropic models were consistently less stiff, but no statistically significant differences in correlation were observed between material models against experimental data. The mean difference in whole bone stiffness between model types was approximately 26%, suggesting that anisotropy can still effect considerable change in the mechanics of proximal femur models. The under prediction of whole bone stiffness in anisotropic models suggests that the orthotropic elastic constants require further investigation. The ability to map mechanical anisotropy from high-resolution images and interpolate information into clinical-resolution models will allow testing of new anisotropic material mapping strategies.

Multi-unit sustained vibration loading platform for biological tissues: design, validation and experimentation.

The relationships between mechanical inputs and resulting biological tissue structure, composition, and metabolism are critical to detailing the nuances of tissue mechanobiology in both healthy and injured tissues. Developing a model system to test the mechanobiology of tissues ex-vivo is a complex task, as controlling chemical and mechanical boundary layers in-vitro are difficult to replicate. A novel multi-unit vibration loading platform for intervertebral discs was designed and validated with both independent electronic data and experimental loading of 6 bovine intervertebral discs (IVDs) and an equal number of unloaded controls. Sustained vibration was applied using closed-loop positional control of pushrods within four independent bioreactors with circulating phosphate buffered saline. The bioreactors were designed to be modular with removable components allowing for easy cleaning and replacement. The loading regime was chosen to maximize target mRNA expression as reported in previous research. Aggrecan, decorin, and versican mRNA all reported statistically significant increases above control levels. Biglycan, collagen type I and II showed no significant difference from the control group. Further study is required to determine the resulting effect of increased mRNA expressions on long-term disc health. However these results indicate that this research is past the proof of concept stage, supporting future studies of mechanobiology utilizing this new device. The next stage in developing this novel loading platform should consider modifying the tissue grips to explore the effects of different directional loading on different gene expression, and also loading different types of tissues.