Validation of 3D finite element models from simulated DXA images for biofidelic simulations of sideways fall impact to the hip

Predicting Tibia-Fibula Geometry and Density From Anatomical Landmarks Via Statistical Appearance Model: Influence of Errors on Finite Element-Calculated Bone Strain

State-of-the-art participant-specific finite element models require advanced medical imaging to quantify bone geometry and density distribution; access to and cost of imaging is prohibitive to the use of this approach. Statistical appearance models may enable estimation of participants' geometry and density in the absence of medical imaging. The purpose of this study was to: (1) quantify errors associated with predicting tibia-fibula geometry and density distribution from skin-mounted landmarks using a statistical appearance model and (2) quantify how those errors propagate to finite element-calculated bone strain. Participant-informed models of the tibia and fibula were generated for thirty participants from height and sex and from twelve skin-mounted landmarks using a statistical appearance model. Participant-specific running loads, calculated using gait data and a musculoskeletal model, were applied to participant-informed and CT-based models to predict bone strain using the finite element method. Participant-informed meshes illustrated median geometry and density distribution errors of 4.39-5.17 mm and 0.116-0.142 g/cm3, respectively, resulting in large errors in strain distribution (median RMSE = 476-492 με), peak strain (limits of agreement =±27-34%), and strained volume (limits of agreement =±104-202%). These findings indicate that neither skin-mounted landmark nor height and sex-based predictions could adequately approximate CT-derived participant-specific geometry, density distribution, or finite element-predicted bone strain and therefore should not be used for analyses comparing between groups or individuals.

3D-DXA Based Finite Element Modelling for Femur Strength Prediction: Evaluation Against QCT

Osteoporosis is characterised by the loss of bone density resulting in an increased risk of fragility fractures. The clinical gold standard for diagnosing osteoporosis is based on the areal bone mineral density (aBMD) used as a surrogate for bone strength, in combination with clinical risk factors. Finite element (FE) analyses based on quantitative computed tomography (QCT) have been shown to estimate bone strength better than aBMD. However, their application in the osteoporosis clinics is limited due to exposure of patients to increased X-rays radiation dose. Statistical modelling methods (3D-DXA) enabling the estimation of 3D femur shape and volumetric bone density from dual energy X-ray absorptiometry (DXA) scan have been shown to improve osteoporosis management. The current study used 3D-DXA based FE analyses to estimate femur strength from the routine clinical DXA scans and compared its results against 151 QCT based FE analyses, in a clinical cohort of 157 subjects. The linear regression between the femur strength predicted by QCT-FE and 3D-DXA-FE models correlated highly (coefficient of determination R2 = 0.86) with a root mean square error (RMSE) of 397 N. In conclusion, the current study presented a 3D-DXA-FE modelling tool providing accurate femur strength estimates noninvasively, compared to QCT-FE models.

3D Finite Element Models Reconstructed From 2D Dual-Energy X-Ray Absorptiometry (DXA) Images Improve Hip Fracture Prediction Compared to Areal BMD in Osteoporotic Fractures in Men (MrOS) Sweden Cohort

Bone strength is an important contributor to fracture risk. Areal bone mineral density (aBMD) derived from dual-energy X-ray absorptiometry (DXA) is used as a surrogate for bone strength in fracture risk prediction tools. 3D finite element (FE) models predict bone strength better than aBMD, but their clinical use is limited by the need for 3D computed tomography and lack of automation. We have earlier developed a method to reconstruct the 3D hip anatomy from a 2D DXA image, followed by subject-specific FE-based prediction of proximal femoral strength. In the current study, we aim to evaluate the method's ability to predict incident hip fractures in a population-based cohort (Osteoporotic Fractures in Men [MrOS] Sweden). We defined two subcohorts: (i) hip fracture cases and controls cohort: 120 men with a hip fracture (<10 years from baseline) and two controls to each hip fracture case, matched by age, height, and body mass index; and (ii) fallers cohort: 86 men who had fallen the year before their hip DXA scan was acquired, 15 of which sustained a hip fracture during the following 10 years. For each participant, we reconstructed the 3D hip anatomy and predicted proximal femoral strength in 10 sideways fall configurations using FE analysis. The FE-predicted proximal femoral strength was a better predictor of incident hip fractures than aBMD for both hip fracture cases and controls (difference in area under the receiver operating characteristics curve, ΔAUROC = 0.06) and fallers (ΔAUROC = 0.22) cohorts. This is the first time that FE models outperformed aBMD in predicting incident hip fractures in a population-based prospectively followed cohort based on 3D FE models obtained from a 2D DXA scan. Our approach has potential to notably improve the accuracy of fracture risk predictions in a clinically feasible manner (only one single DXA image is needed) and without additional costs compared to the current clinical approach. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

The accuracy of statistical shape models in predicting bone shape: A systematic review

BACKGROUND

This systematic review aims to ascertain how accurately 3D models can be predicted from two-dimensional (2D) imaging utilising statistical shape modelling.

METHODS

A systematic search of published literature was conducted in September 2022. All papers which assessed the accuracy of 3D models predicted from 2D imaging utilising statistical shape models and which validated the models against the ground truth were eligible.

RESULTS

2127 papers were screened and a total of 34 studies were included for final data extraction. The best overall achievable accuracy was 0.45 mm (root mean square error) and 0.16 mm (average error).

CONCLUSION

Statistical shape modelling can predict detailed 3D anatomical models from minimal 2D imaging. Future studies should report the intended application domain of the model, the level of accuracy required, the underlying demographics of subjects, and the method in which accuracy was calculated, with root mean square error recommended if appropriate.

The efficacy of femoral augmentation for hip fracture prevention using ceramic-based cements: A preliminary experimentally-driven finite element investigation

Femoral fractures due to sideways falls continue to be a major cause of concern for the elderly. Existing approaches for the prevention of these injuries have limited efficacy. Prophylactic femoral augmentation systems, particularly those involving the injection of ceramic-based bone cements, are gaining more attention as a potential alternative preventative approach. We evaluated the mechanical effectiveness of three variations of a bone cement injection pattern (basic ellipsoid, hollow ellipsoid, small ellipsoid) utilizing finite element simulations of sideways fall impacts. The basic augmentation pattern was tested with both high- and low-strength ceramic-based cements. The cement patterns were added to the finite element models (FEMs) of five cadaveric femurs, which were then subject to simulated sideways falls at seven impact velocities ranging from 1.0 m/s to 4.0 m/s. Peak impact forces and peak acetabular forces were examined, and failure was evaluated using a strain-based criterion. We found that the basic HA ellipsoid provided the highest increases in both the force at the acetabulum of the impacted femur ("acetabular force", 55.0% ± 22.0%) and at the force plate ("impact force", 37.4% ± 15.8%). Changing the cement to a weaker material, brushite, resulted in reduced strengthening of the femur (45.2% ± 19.4% acetabular and 30.4% ± 13.0% impact). Using a hollow version of the ellipsoid appeared to have no effect on the fracture outcome and only a minor effect on the other metrics (54.1% ± 22.3% acetabular force increase and 35.3% ± 16.0% impact force increase). However, when the outer two layers of the ellipsoid were removed (small ellipsoid), the force increases that were achieved were only 9.8% ± 5.5% acetabular force and 8.2% ± 4.1% impact force. These results demonstrate the importance of supporting the femoral neck cortex to prevent femoral fractures in a sideways fall, and provide plausible options for prophylactic femoral augmentation. As this is a preliminary study, the surgical technique, the possible effects of trabecular bone damage during the augmentation process, and the effect on the blood supply to the femoral head must be assessed further.

2D-3D reconstruction of the proximal femur from DXA scans: Evaluation of the 3D-Shaper software

Introduction: Osteoporosis is currently diagnosed based on areal bone mineral density (aBMD) computed from 2D DXA scans. However, aBMD is a limited surrogate for femoral strength since it does not account for 3D bone geometry and density distribution. QCT scans combined with finite element (FE) analysis can deliver improved femoral strength predictions. However, non-negligible radiation dose and high costs prevent a systematic usage of this technique for screening purposes. As an alternative, the 3D-Shaper software (3D-Shaper Medical, Spain) reconstructs the 3D shape and density distribution of the femur from 2D DXA scans. This approach could deliver a more accurate estimation of femoral strength than aBMD by using FE analysis on the reconstructed 3D DXA. Methods: Here we present the first independent evaluation of the software, using a dataset of 77 ex vivo femora. We extend a prior evaluation by including the density distribution differences, the spatial correlation of density values and an FE analysis. Yet, cortical thickness is left out of this evaluation, since the cortex is not resolved in our FE models. Results: We found an average surface distance of 1.16 mm between 3D DXA and QCT images, which shows a good reconstruction of the bone geometry. Although BMD values obtained from 3D DXA and QCT correlated well (r ² = 0.92), the 3D DXA BMD were systematically lower. The average BMD difference amounted to 64 mg/cm³, more than one-third of the 3D DXA BMD. Furthermore, the low correlation (r ² = 0.48) between density values of both images indicates a limited reconstruction of the 3D density distribution. FE results were in good agreement between QCT and 3D DXA images, with a high coefficient of determination (r ² = 0.88). However, this correlation was not statistically different from a direct prediction by aBMD. Moreover, we found differences in the fracture patterns between the two image types. QCT-based FE analysis resulted mostly in femoral neck fractures and 3D DXA-based FE in subcapital or pertrochanteric fractures. Discussion: In conclusion, 3D-Shaper generates an altered BMD distribution compared to QCT but, after careful density calibration, shows an interesting potential for deriving a standardized femoral strength from a DXA scan.

The Influence of Fall Direction and Hip Protector on Fracture Risk: FE Model Predictions Driven by Experimental Data

Hip fractures in older adults, which often lead to lasting impairments and an increased risk of mortality, are a major public health concern. Hip fracture risk is multi-factorial, affected by the risk of falling, the load acting on the femur, and the load the femur can withstand. This study investigates the influence of impact direction on hip fracture risk and hip protector efficacy. We simulated falls for 4 subjects, in 7 different impact directions (15° and 30° anterior, lateral, and 15°, 30°, 60°, and 90° posterior) at two different impact velocities (2.1 and 3.1 m/s), all with and without hip protector, using previously validated biofidelic finite element models. We found the highest number of fractures and highest fragility ratios in lateral and 15° posterior impacts. The hip protector attenuated femur forces by 23–49 % for slim subjects under impact directions that resulted in fractures (30° anterior to 30° posterior). The hip protector prevented all fractures (6/6) for 2.1 m/s impacts, but only 10% of fractures for 3.1 m/s impacts. Our results provide evidence that, regarding hip fracture risk, posterior-lateral impacts are as dangerous as lateral impacts, and they support the efficacy of soft-shell hip protectors for anterior- and posterior-lateral impacts.

Finite Element Assessment of Bone Fragility from Clinical Images

PURPOSE OF REVIEW

We re-evaluated clinical applications of image-to-FE models to understand if clinical advantages are already evident, which proposals are promising, and which questions are still open.

RECENT FINDINGS

CT-to-FE is useful in longitudinal treatment evaluation and groups discrimination. In metastatic lesions, CT-to-FE strength alone accurately predicts impending femoral fractures. In osteoporosis, strength from CT-to-FE or DXA-to-FE predicts incident fractures similarly to DXA-aBMD. Coupling loads and strength (possibly in dynamic models) may improve prediction. One promising MRI-to-FE workflow may now be tested on clinical data. Evidence of artificial intelligence usefulness is appearing. CT-to-FE is already clinical in opportunistic CT screening for osteoporosis, and risk of metastasis-related impending fractures. Short-term keys to improve image-to-FE in osteoporosis may be coupling FE with fall risk estimates, pool FE results with other parameters through robust artificial intelligence approaches, and increase reproducibility and cross-validation of models. Modeling bone modifications over time and bone fracture mechanics are still open issues.

Statistical Shape and Appearance Models: Development Towards Improved Osteoporosis Care

PURPOSE OF REVIEW

Statistical models of shape and appearance have increased their popularity since the 1990s and are today highly prevalent in the field of medical image analysis. In this article, we review the recent literature about how statistical models have been applied in the context of osteoporosis and fracture risk estimation.

RECENT FINDINGS

Recent developments have increased their ability to accurately segment bones, as well as to perform 3D reconstruction and classify bone anatomies, all features of high interest in the field of osteoporosis and fragility fractures diagnosis, prevention, and treatment. An increasing number of studies used statistical models to estimate fracture risk in retrospective case-control cohorts, which is a promising step towards future clinical application. All the reviewed application areas made considerable steps forward in the past 5-6 years. Heterogeneities in validation hinder a thorough comparison between the different methods and represent one of the future challenges to be addressed to reach clinical implementation.

Prophylactic augmentation implants in the proximal femur for hip fracture prevention: An in silico investigation of simulated sideways fall impacts

Femoral fractures from sideways falls in the elderly are associated with significant rates of morbidity and mortality. Approaches to prevent these catastrophic injuries include pharmacological treatments, which have limited efficacy. Prophylactic femoral augmentation systems are a promising alternative that are gaining prominence by addressing the most debilitating osteoporosis-related fracture. We have developed finite element models (FEMs) of a novel experimental sideways fall simulator for cadavers. By virtue of the range of specimens and injury outcomes, these FEMs are well-suited to the evaluation of such implants. The purpose of this study was to use the FEMs to evaluate the mechanical effectiveness of three different prophylactic femoral augmentation systems. Models of the Y-Strut® (Hyprevention®, Pessac, France), Gamma Nail® (Stryker, Kalamazoo, USA), and a simple lag screw femoral fracture implant systems were placed into FEMs of five cadaveric pelvis-femur constructs embedded in a soft tissue surrogate, which were then subject to simulated sideways falls at seven impact velocities. Femur-only FEMs were also evaluated. Peak impact forces and peak acetabular forces were examined, and failure was evaluated using a strain-based criterion. We found that the femoral augmentation systems increased the peak forces prior to fracture, but were unable to prevent fracture for severe impacts. The Gamma Nail® system consistently produced the largest strength increases relative to the unaugmented femur for all five specimens in both the pendulum-drop FEMs and the femur-only simulations. In some cases, the same implant appeared to cause fractures in the acetabulum. The femur-only FEMs showed larger force increases than the pendulum-drop simulations, which suggests that the results of the femur-only simulations may not represent sideways falls as accurately as the soft tissue-embedded pendulum-drop simulations. The results from this study demonstrate the ability to simulate a high energy phenomenon and the effect of implants in an in silico environment. The results also suggest that implants could increase the force applied to the proximal femur during impact. Fracture outcomes from the tested implants can be used to inform the design of future devices, which reaffirms the value of modelling with biofidelic considerations in the implant design process. To the authors' knowledge, this is the first paper to use more complex biofidelic FEMs to assess prophylactic femoral augmentation methods.