Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 23 October 2025

Sec. Biomechanics

Volume 13 - 2025 | https://doi.org/10.3389/fbioe.2025.1673339

Inelastic mechanical descriptors for osteoporotic hip fracture discrimination with 3D-DXA-based nonlinear finite element models

  • 1BCN MedTech, Department of Engineering, Pompeu Fabra University, Barcelona, Spain
  • 2CETIR Ascires, Barcelona, Spain
  • 33D-Shaper Medical, Barcelona, Spain

Introduction: Dual-energy X-ray absorptiometry (DXA) is the gold standard for diagnosing osteoporosis. Advances in 2D-3D modelling to generate patient-specific 3D-DXA models out of DXA images enable accurate volumetric representations of the femur, with potential for fracture risk prediction when combined with finite element (FE) analyses. This study evaluates the ability of 3D-DXA-based FE models to discriminate hip fractures under side-fall loading.

Methods: We used a retrospective case-control study including 128 women, 64 of whom suffered a hip fracture. Mechanical descriptors, including strength, nonlinear deformation, residual displacement, and energy absorption under elastic-plastic assumptions, were derived from force-displacement curves.

Results: The area under the receiver operating characteristic curve (AUROC) showed that strength and trabecular volumetric bone mineral density (vBMD) equally discriminated between fracture and control subjects. Residual displacement due to plastic strain accumulation at failure emerged as a key descriptor which, when combined with strength, significantly improved fracture discrimination (ΔAUROC = 0.11 vs. areal bone mineral density (aBMD); ΔAUROC = 0.08 vs. trabecular vBMD).

Discussion: These findings highlight the potential of 3D-DXA and FE modelling to improve fracture assessment within current DXA-based clinical workflows.

1 Introduction

Osteoporosis is identified by a loss of bone mass leading to bone strength reduction. It evolves silently, and its late diagnosis and/or lack of prevention may lead to fragile fractures, loss of life quality, as well as a significant excess of mortality for the patients (Kanis et al., 2021).

Osteoporosis is currently diagnosed by using the areal bone mineral density (aBMD) computed from 2D dual-energy X-ray absorptiometry (DXA) images of the proximal femur and lumbar spine (Stone et al., 2003). However, aBMD is an insufficient metric to estimate the risk of a fragile fracture since it captures neither the distribution of volumetric bone mineral density (vBMD) nor other important factors that control the concentration of mechanical loads in bone tissues, such as bone geometry or external loads. Predictive tools such as FRAX and Garvan (Adler, 2014; Watts et al., 2009; McCloskey et al., 2012) were proposed to estimate the fracture risk, considering different clinical risk factors (age, weight, height, sex, BMI, current smoking status, etc.) combined with aBMD or alone. Although these tools consider multiple risk factors for fracture assessment, they still do not account for bone geometry, spatial distribution of bone density, and mechanical loads. Developing accurate measures of bone strength and mechanical behavior is an essential step toward enhancing fracture risk prediction tools.

In recent studies, subject-specific finite element (FE) models based on quantitative computed tomography (QCT) or DXA scans were proposed, as promising tools to improve fracture risk evaluation, as they allow bone strength estimations by integrating physics-based rationales (Cody et al., 1999). Mechanical experimental models of the proximal femur load and fracture under side-way fall conditions (Dragomir-Daescu et al., 2011; Kok et al., 2021; Courtney et al., 1995; Grassi et al., 2012; Dall’Ara et al., 2016) were developed to validate the aforementioned FE models.

QCT-based subject-specific FE models led to improved fracture risk evaluation compared to aBMD in different studies (Schileo and Taddei, 2021; Fleps and Morgan, 2022). Fleps et al. (2022) used nonlinear QCT-based FE models to assess the fracture risk of the proximal femur under multiple loading conditions by developing an FE modeling pipeline previously established by Enns-Bray et al. (2019). They used nonlinear and linear material models in their study. The results demonstrate that the non-linear FE model was a better approach to discriminate the fracture risk compared to linear FE models. Loading mechanisms with low adduction and internal rotation were found to be more suitable for fracture classification. Cao et al. (2024) presented a QCT-based FE model and computed the yielding load, strength, and energy to failure. They combined the FE parameters with aBMD and six demographical and clinical parameters to develop an FE-based fracture risk index using four loading mechanisms. The results showed that mechanical enrichment of aBMD information through FE parameters improved femoral fracture risk assessment compared to aBMD adjusted for six covariates. Yet, the high cost and radiation dose of QCT limit the clinical applicability of QCT-based FE approaches.

Since DXA is the current clinical standard for osteoporosis management, DXA-based 2D FE models were proposed to predict bone strength (Dall’Ara et al., 2016) and fracture risk in clinical cohorts (Terzini et al., 2019; Leslie et al., 2019; Yang et al., 2018). Dall’Ara et al. (2016) combined experimental investigation on the proximal femur in a standing configuration and a fall onto the greater trochanter, with subject-specific DXA-based 2D linear FE models predicting the bone strength and failure criteria in good agreement with the experimental results. Ulivieri and Rinaudo (2022) introduced the Bone Strain Index (BSI) as a metric for assessing the load resistance of the femur and lumbar spine using DXA-based 2D FE models. Their approach employed linear elastic FE models, where the forces acting on the greater trochanter were determined by the patient’s weight and height. The BSI values represent the average bone strain at the analyzed site (neck BSI and total hip BSI), with higher BSI values indicating increased bone strain and, consequently, a higher fracture risk. Since 2018, several clinical studies have explored the clinical utility of BSI in identifying osteoporotic patient subgroups having a particular tendency to fragility fractures (Rinaudo et al., 2024; Ulivieri et al., 2022; Sornay-Rendu et al., 2022). Naylor et al. (2013) used a DXA-based 2D FE model to assess hip fractures in 728 female subjects. The area under the receiver operating characteristic curve (AUROC) calculated using the load-to-strength ratio (LSR) combined with femoral neck (FN) aBMD was significantly greater compared to FN aBMD alone. Although DXA-based 2D FE simulations could reflect mechanical behaviors (Dall’Ara et al., 2016), their ability for osteoporosis fracture classification remains limited (Naylor et al., 2013). One of the limitations of the DXA-based 2D FE method is the inability to simulate sideways fall models using different loading configurations (Falcinelli et al., 2014). Moreover, in these DXA-based 2D FE models, the volumetric information is missing, and planar information may not sufficiently encompass the out-of-plane distribution of shapes and bone mineral density, depending on the specific characteristics of a patient’s femur.

3D-DXA-based FE models have the potential to overcome the limitations of DXA-based 2D FE models (Thevenot et al., 2014; Väänänen et al., 2015; Grassi et al., 2017; Grassi et al., 2021; Wills et al., 2019). Grassi et al. employed statistical shape and appearance models (SSAM) to reconstruct the three-dimensional anatomy of individual femurs from DXA images, subsequently conducting 3D FE simulations (Grassi et al., 2023). The AUROC values computed using bone strength values normalized to the subject’s weight in 10 different side fall configurations revealed significant improvement compared to aBMD + body mass index (BMI) criteria for hip fracture discrimination (Grassi et al., 2023; Grassi et al., 2025). Although most QCT- and 3D-DXA-based FE models improved fracture risk discrimination compared to aBMD, their implementation in clinical practices remains challenging due to high FE simulation cost, lack of automation (Viceconti et al., 2018), or absence of regulatory approval for use in clinical practice (Grassi et al., 2023). Moreover, assessing a robust mechanical criterion that incorporates the effective mechanical parameters associated with femur fracture could enhance the discrimination of hip fractures in a case-control study.

3D-Shaper® is a software solution that implements the 3D-DXA technology and has obtained regulatory approvals for clinical use (Humbert et al., 2016). Hip DXA scans are analyzed automatically to provide a 3D analysis of the cortical and trabecular compartments. A 3D-DXA-based FE pipeline by using the 3D-Shaper software was used in our previous study (Wills et al., 2019) in which a static peak load that depended on patient mass and height was applied to the femur head. Discrimination of fracture subjects through the local mechanical descriptors was performed in the study, which showed that the major principal stress was a better discriminator compared to the vBMD measured at the same region of interest. Moreover, high correlations were reported between femur strength estimated using the 3D-DXA-based FE model and QCT-based FE models (Dudle et al., 2023; Qasim et al., 2024). The mean computation time per subject for the 3D-DXA algorithm was 1 min 30 s (Humbert et al., 2016) and 15 min for the FE simulation (Qasim et al., 2024), making this approach potentially suitable for use in clinical settings, compared to the other 3D-DXA-based FE algorithms (Grassi et al., 2023).

In the current study, this 3D-DXA-based nonlinear FE model (Qasim et al., 2024) is used to assess the ability of global mechanical descriptors to discriminate between hip fractures and controls. These descriptors integrate the effects of organ geometry, distribution of bone mineral density and mechanical properties, and tissue-dependent organ structure, in 3D patient-specific representations of the proximal femur.

2 Materials and methods

2.1 Cohort and DXA images

We used data from a retrospective clinical cohort collected in a previous study (Humbert et al., 2020). The retrospective clinical cohort study was performed at CETIR Group Mèdic (Barcelona, Spain) and included 64 subjects with incident hip fractures and 64 controls. The study was conducted under the latest version of the Declaration of Helsinki. Ethical approval was granted by the CETIR Group Mèdic Scientific Committee. Informed consent was obtained from all participants. Subject anonymity was maintained by assigning subject-specific numeric codes to all records. The cohort comprised Caucasian post-menopausal women with baseline visits between 2000 and 2011. Subjects with any prior osteoporotic fractures or diseases other than osteoporosis that affect bone metabolism were excluded from the study. Hip fracture was diagnosed through clinical history or a phone call in case of missing clinical information. A review of clinical history and DXA scans of subjects confirmed the absence of fracture at the baseline and follow-up period in sex- and age-matched subjects in the control group. There was no significant difference (P-value ≥0.05) between age, height, weight, and body mass index (BMI) between subjects in the hip fracture and control groups (Table 1).

Table 1
www.frontiersin.org

Table 1. Mean ± standard deviation characteristics in control and hip fracture groups (Humbert et al., 2020).

DXA scans were conducted using a Prodigy scanner (GE Healthcare, Madison, WI, USA) following the recommendations provided by the manufacturer. aBMD (g/cm2), bone mineral content (g), and the area (cm2), were measured for the entire scanned region of the proximal femur by a trained radiologist using enCORE software v14.10 (GE Healthcare) (Humbert et al., 2020).

2.2 3D-DXA analysis

3D-Shaper software v2.14 (3D-Shaper Medical, Barcelona, Spain) was used to analyze the DXA scan and provide patient-specific 3D-DXA geometrical and bone density models. In brief, the software uses an SSAM of the proximal femur, built out of a database of QCT images and registered onto the DXA scan, to estimate a patient-specific 3D femur shape and volumetric density distribution (appearance) (Humbert et al., 2016). This approach allows for representation of both the cortical and trabecular compartments (Humbert et al., 2020), which were subsequently used to generate FE structural meshes with vBMD-based material properties specific to each tissue (Qasim et al., 2024).

Integral vBMD (mg/cm3) was calculated as the mean vBMD of the cortical and trabecular compartments in the total hip (TH). Trabecular vBMD mg/cm3, was calculated as the average vBMD of the trabecular compartment at the TH (Humbert et al., 2016). Cortical surface BMD (sBMD) (mg/cm2), was computed as the multiplication of the cortical vBMD by the cortical thickness, both averaged over the TH region of interest. The accuracy of the 3D-Shaper modeling methods was evaluated in different studies (Humbert et al., 2016; Dudle et al., 2023). It is noteworthy that in 3D-DXA analyses, the proximal femur, excluding the pelvis and surrounding structures, was considered as the primary region of interest, as it more directly corresponds to the anatomical site relevant for hip fracture risk assessment.

2.3 3D-DXA-based nonlinear FE modeling

A 3D-DXA-based nonlinear FE model previously described (Qasim et al., 2024) was utilized to simulate the mechanical behavior of the femur under sideways fall. In brief, the 3D proximal femur anatomy derived from the 3D-Shaper software (Humbert et al., 2016) served for the automatic personalization of FE models. A generic structural FE mesh was morphed into the subject-specific femur shape using thin-plate splines (TPS) transformations. The TPS transformations ensured that separate FE mesh components were consistently morphed into the cortical and trabecular compartments of each subject-specific model. Linear hexahedral eight-node brick elements were used in the FE model.

Using linear interpolations of the bone density values in the image provided by the 3D-Shaper software, the subject-specific bone densities were mapped onto the FE meshed model. In each tissue compartment, the elastic modulus of each element was determined by using the density-elasticity relationships (Ho et al., 1992; Keller, 1994) presented in Equation 1 where the elastic modulus E was in MPa, the ash density ρash in gcm3, and the apparent density ρapp in kgm3.

ETrabecular=0.003715×ρapp1.96,ECortical=10200×ρash2.01(1)

The Poisson’s ratio was assumed to be 0.3. The ash density of the bone, in gcm3, was calculated as follows (Equation 2) (Schileo et al., 2008a):

ρash=0.87ρQCT+0.079(2)

where ρQCT refers to the radiological density, which was calculated using the vBMD values (gcm3) derived from the 3D-Shaper modeling process. To access the apparent density in Equation 1, we assumed ρashρapp ratio of 0.6, in the range of values (0.55–0.63) provided by available studies, where both ρash and ρapp were in gcm3 (Schileo et al., 2008a; Goulet et al., 1994; Keyak et al., 1994).

To simulate the nonlinear behavior of bone, an elastic-perfect plastic model was applied (Johannesdottir et al., 2017; Bayraktar et al., 2004; Bevill et al., 2007; Morgan and Keaveny, 2001):

σycTrabecular=38.5×ρapp1.48,σytTrabecular=22.6×ρapp1.26σycCortical=041+0.0062×Ecortical,σytCortical=0.33+.0039×Ecortical

where σyc and σyt represent the compressive and tensile yield stress, respectively.

Boundary conditions that mimic a state-of-the-art experimental model of lateral fall with femur cadaveric experiments (Dragomir-Daescu et al., 2011; Kok et al., 2021; Courtney et al., 1995; Grassi et al., 2012; Dall’Ara et al., 2016), recurrently used in FE models for femur strength calculation (Grassi et al., 2012; Grassi et al., 2023; Abe et al., 2022) were employed in the 3D-DXA-based FE simulations. The end of the modelled shaft region was fully constrained, and the surface of the greater trochanter was fixed only in the fall direction (Figure 1a). A rigid-body constraint imposed a displacement on the femur head through a reference point (RP), until the reaction force started to drop with increasing displacements. Such a drop was considered as a surrogate of mechanical failure, following the nonlinear apparent elastic-plastic response of the femur model. The RP and the corresponding slave nodes on the femur head are shown in Figure 1a. The strength (F0) was defined as the corresponding maximum force, and both F0 and the corresponding displacement, D0, were used as the main output parameters of the 3D-DXA-based nonlinear FE model (Figure 1b). Furthermore, to account for variability in impact direction, different sideways fall configurations covering 0°–30° of adduction and internal rotation were also modeled in the study (Falcinelli et al., 2014; Grassi et al., 2023).

Figure 1
Image consists of two parts: (a) A green 3D mesh model representing a bone, with an arrow indicating the direction of fall labeled

Figure 1. 3D-DXA FE femur model under side-way fall direction (a) and FE Strength and displacement on force-displacement plot (b).

Figure 2 provides a schematic representation of the conceptual framework adopted in this study.

Figure 2
Flowchart depicting a process for hip fracture discrimination, starting with DXA scans for aBMD, followed by SSAM 2D-to-3D reconstruction to obtain a subject-specific 3D shape. The process involves TPS-based mesh morphing using a 3D generic FE mesh. Material properties, using the volumetric density distribution and density-elasticity relationships, are assigned, leading to a subject-specific FE model. Outputs include FE-derived mechanical descriptors under side fall simulations. A graph compares true positive and false positive rates between FEM and aBMD, with AUROC values of 0.80 and 0.69, respectively.

Figure 2. Schematic representation of the concept underlying this study.

2.4 Mechanical descriptors for fracture and non-fracture discrimination

The strength (F0), displacement at failure (D0), and the nonlinear energy up to the failure point (Figure 1b) were used as nonlinear mechanical descriptors for fracture group discrimination.

To quantify the different characteristics of the elastic-plastic behavior, the nonlinear FE-derived curves were post-processed using an optimized piecewise linear method, described in detail in the Supplementary Appendix A. The piece-wise approximation provided a series of additional mechanical descriptors (Figure 3). These included: the linear elastic force (F1); the linear elastic displacement (D1); the area under the linear elastic region of the bilinear curve (linear elastic energy); the nonlinear deformation (D2) (Figure 3a). A loading/unloading analysis was also conducted on the nonlinear FE models to assess the mechanical energy dissipated by the accumulated plastic deformations. Key parameters from this analysis, including residual displacement (D3), dissipated energy per se, and residual energy (elastic energy), are illustrated in Figure 3b. The full list of mechanical descriptors and combinations thereof, used for the discrimination of fragile fracture in the cohort, is provided in Table 2.

Figure 3
Two graphs illustrate force versus displacement. Graph (a) shows a bilinear curve with an area shaded in purple representing linear elastic energy. The curve initially increases sharply then flattens. Graph (b) displays a similar curve under loading/unloading process, with areas shaded in pink and purple, representing dissipated and residual energy, respectively. Connecting dots indicates nonlinear FE simulation results. Horizontal axes show displacement in millimeters, and vertical axes show force in newtons.

Figure 3. Piecewise linear analysis, (a) bilinear force-displacement curve of the femur, (b) energy evolution characteristic of the FE analysis.

Table 2
www.frontiersin.org

Table 2. FE-based mechanical descriptors obtained by the nonlinear analysis, piecewise linear analysis, and their combination.

Among the analysed parameters, strength is the maximum load sustained by the femur before global failure, reflecting overall bone structural integrity. Displacement at maximum reaction force denotes the volume fraction of elements that have transitioned into the plastic region, indicating the bone’s capacity to deform prior to fracture. Nonlinear energy quantifies the total mechanical energy absorbed up to failure, integrating elastic and plastic contributions, and reflects intrinsic bone toughness. Greater energy absorption is protective, as bone capable of dissipating higher mechanical energy during a fall is less likely to fracture under equivalent strength.

Within the elastic regime, linear elastic force corresponds to the elastic mechanical energy of elements remaining in the elastic region, whereas linear elastic displacement represents their displacement, indicating initial and functional bone deformability. These two parameters characterize bone stiffness; lower force at a given displacement, or greater displacement at a given force, reflects reduced stiffness, linked to thinner cortices or lower trabecular density. Linear elastic energy quantifies the fraction of total mechanical energy absorbed elastically, highlighting the bone’s capacity to resist loading prior to plastic deformation, and reflects its ability to store reversible energy, closely associated with stiffness and load-bearing capacity.

Beyond the elastic region, nonlinear deformation, calculated as the difference between maximum displacement at the failure point and linear elastic displacement, represents the post-yield ductility of the model, potentially mitigating fracture risk despite low bone strength. Residual displacement serves as a surrogate for the elements, i.e., the volume, that have undergone irreversible failure or plastic deformation due to local stress concentrations. Dissipated energy corresponds to mechanical energy due to plastic deformations. Residual energy, i.e., the energy retained after unloading, provided a quantification of how much mechanical energy could remain at failure, according to the strength model.

The maximum load sustained by the femur, captured as strength in the nonlinear analysis and as linear elastic force in the piecewise linear analysis, was further combined with deformation- and energy-based metrics. It aimed to represent a more comprehensive assessment of bone mechanical behaviour, reflecting a personalised balance of stored and dissipated energy along the way to failure. The mechanical description and clinical relevance of Strength-Based- Combinations and Linear Elastic Force- Based Combination are presented in Supplementary Appendix B (Supplementary Tables SB1, SB2).

2.5 Statistical analysis

Logistic regression and AUROC were used in this study to evaluate how the descriptors of the 3D-DXA-based FE model discriminate between fracture subjects and controls, compared to aBMD and 3D-DXA cortical and trabecular parameters. Furthermore, for variables that followed a normal distribution, we applied the unpaired two-sample t-test and for those that did not meet the assumption of normality, we used the non-parametric Mann–Whitney U test to assess whether the mean values of the mechanical parameters in the case group differed significantly from those in the control group. The mechanical descriptors were combined in a multi-logistic regression analysis for fracture assessment. A leave-one-out cross-validation (Airola et al., 2011; Stone, 1974) was used in all logistic regressions, and the statistical significance of the differences among the obtained AUROC values was evaluated using the DeLong method (DeLong et al., 1988).

3 Results

3.1 Mechanical behavior of the femur under side-way fall configuration

Figure 4 illustrates the force-displacement curves for 3D-DXA-based femur models simulated under side-fall loading, in randomly selected subjects. The curves show that subjects with comparable bone strength could reach different failure points, characterized by varying displacement (D0) and nonlinear energy values.

Figure 4
Line graph showing the relationship between force (N) and displacement (mm). Multiple colored lines represent different subjects, each labeled with a unique number. Two subjects shown in the same color indicate identical strength but different final displacement and corresponding energy.

Figure 4. Force-displacement curves of a few subjects under sideway fall FE simulation.

3D-DXA-based FE analysis results, depicting the distribution of maximum principal plastic strain and maximum principal stress at maximum reaction force, in a randomly selected femur model under side-fall loading is provided in the Supplementary Appendix C. High values of both maximum principal plastic strain and stress were observed at the femoral neck, suggesting this region as the most probable site of structural failure.

3.2 FE-predicted mechanical descriptors and their ability to discriminate fracture and non-fracture subjects

The mean ± standard deviation of TH aBMD, 3D-DXA measurements, and mechanical descriptors calculated using the FE analysis for the hip fracture and control groups are presented in Table 3. Parameters including aBMD, integral vBMD, trabecular vBMD, cortical sBMD, strength, linear elastic force, linear elastic energy, and residual energy, were significantly lower in hip fracture cases compared to controls. In contrast, displacement, linear elastic displacement, nonlinear deformation, and residual displacement were significantly higher in hip fracture cases. Nonlinear energy was higher in controls, while dissipated energy was greater in hip fracture cases, although none of these differences were statistically significant.

Table 3
www.frontiersin.org

Table 3. Mean ± standard deviation of DXA, 3D-DXA measurement, and FE-based mechanical descriptors in the control and hip fracture groups.

The AUROCs achieved with the FE–based mechanical descriptors were compared to those achieved with aBMD and 3D-DXA measurements, i.e., integral vBMD, trabecular vBMD, and cortical sBMD (Table 4). The FE-based mechanical descriptors were classified into four categories: (i) Nonlinear FE Analysis; (ii) Piecewise Linear Analysis; (iii) Combination of Bone Strength with Mechanical Descriptors; (iv) Combination of Linear Elastic Force with Mechanical Descriptors.

1. Nonlinear FE Analysis: Nonlinear FE analysis identified displacement as a key mechanical descriptor for fracture discrimination. However, strength alone showed slightly higher discrimination power compared to aBMD and comparable performance to trabecular vBMD.

2. Piecewise Linear Analysis: Piecewise linear analysis generated two categories of mechanical descriptors based on elastic-plastic characteristics: a) linear elastic-based descriptors, including elastic force, displacement, and energy; b) nonlinear plastic-based descriptors, such as nonlinear deformation, residual displacement, and residual energy. Mechanical descriptors derived from these analyses, particularly nonlinear deformation and residual displacement, demonstrated a statistically significant enhancement in discriminating hip fracture cases from controls compared to the aBMD method [AUROC: 0.79 versus 0.69]. Additionally, linear elastic force and linear elastic displacement provided a discrimination capacity similar to trabecular vBMD [AUROC:0.72].

3. Combination of Bone Strength with Mechanical Descriptors: Combining FE strength with corresponding displacement, nonlinear energy, and nonlinear plastic-based descriptors significantly enhanced fracture discrimination compared to BMD-based methods. Notable combinations that significantly over-performed any aBMD-based descriptors included: Strength + Displacement, Strength + Nonlinear Energy, Strength + Nonlinear Deformation, Strength + Residual Displacement, and Strength + Dissipated Energy. Among these, FE strength + Residual Displacement emerged as the most effective combination, achieving a notable improvement in AUROC compared to both aBMD (0.80 vs. 0.69, P = 0.001) and trabecular vBMD (0.80 vs. 0.72, P = 0.016).

4. Combination of Linear Elastic Force with Mechanical Descriptors: Integrating linear elastic force with the corresponding linear elastic displacement, energy, and nonlinear plastic-based descriptors further improved the ability to discriminate fracture cases from controls across most combinations. Particularly effective combinations included: Linear Elastic Force + Linear Elastic Displacement, Linear Elastic Force + Linear Elastic Energy, Linear Elastic Force + Nonlinear Deformation, Linear Elastic Force + Residual Displacement, and Linear Elastic Force + Dissipated Energy. Notably, Linear Elastic Force + Residual Displacement showed a statistically significant improvement in fracture discrimination compared to aBMD and 3D-DXA measurements [0.79 versus 0.69 and 0,72, respectively], further highlighting the importance of incorporating nonlinear-plastic based descriptors into fracture risk models.

Table 4
www.frontiersin.org

Table 4. AUROC calculated using leave-one-out cross-validation.

The ROC curve obtained using strength + residual displacement, which resulted in the highest AUROC, is presented in Figure 5, together with the ROC curves obtained using aBMD and trabecular vBMD. The FE-derived mechanical descriptors for the control and fracture subjects are presented in Supplementary Appendix D (Supplementary Tables SD1, SD2). It is noteworthy that the influence of different fall scenarios and the impact thereof on the discrimination ability is summarized in Supplementary Appendix E. The presented scenarios correspond to those most commonly applied in previous investigations (Falcinelli et al., 2014; Grassi et al., 2023; Pinilla et al., 1996; Abe et al., 2022). Remarkably, the capacity of the femur strength parameter derived from FE simulations to overcome the capacity of the discrimination by aBMD (AUROC = 0.69) depended on the loading direction. Strength AUROC values ranged from 0.68 to 0.72. The best performance of the strength descriptor was achieved with a load aligned with the vertical when the fall models had no internal rotation and no adduction. In every case, Strength AUROC values were overperformed when mechanical information about the inelastic model response was considered for discrimination, according to the best-performing descriptors achieved with the default load.

Figure 5
Two graphs show Receiver Operating Characteristic (ROC) curves. Graph (a) compares FE Strength+FE Residual Displacement (pink line, AUROC = 0.80) and aBMD (blue line, AUROC = 0.69). Graph (b) compares FE Strength+FE Residual Displacement (pink line, AUROC = 0.80) and Trabecular vBMD (blue line, AUROC = 0.72). Both graphs plot True Positive Rate against False Positive Rate.

Figure 5. ROC curves for the mechanical descriptors compared with (a) aBMD and (b) Trabecular vBMD.

4 Discussion

This study explored the added value of the global mechanical descriptors to discriminate the occurrence of femur fragile fracture in a case-control study using DXA-based 3D modeling and FE methods.

Three sets of FE-derived descriptors were chosen to capture distinct mechanical and possibly clinically meaningful aspects of hip fracture risk, as far as fracture risk groups can be compared. Importantly, in contrast to sophisticated FE models that describe bone failure mechanics (Zysset et al., 2015), we propose mechanical surrogates as per the capacity of the femoral bone to absorb mechanical energy in a binary partition of osteoporosis patients: fracture cases, and controls. In our models, nonlinear FE-based descriptors, extracted from the load-displacement curve up to failure, capture the volumetric accumulation of plastic strain up to the point at which the femur model stops to increasingly resist the load. We see from our results (Table 4, Strength-based + nonlinear descriptors) that the nonlinearity associated with this simulated process nicely complements the state-of-the-art bone strength, to eventually overperform the discrimination capacity of both DXA-based aBMD and 3D-DXA-based vBMD. The piecewise modelling of the nonlinear mechanical response further informed about how a specific femur model was reaching its maximum ability to store mechanical energy: it was used to evaluate the balance of elastic energy stored vs. dissipating plastic deformations until the maximum reaction force was reached. Especially, it gave access to the residual displacement descriptor that, once combined with the strength, led to the best discrimination capacity, and proved to be, per se a promising descriptor for fracture risk assessment. Overall, these results, achieved with surrogates of the inelastic behaviour of bone, suggest that quantifying the bone’s post-yield behaviour, according to measurable bone properties at baseline, could improve significantly the clinical assessment of the risk of bone fragile fracture.

The proposed analyses underscore two critical aspects of fracture risk assessment. First, relying on bone strength alone is insufficient for patient-specific fracture discrimination, as it overlooks critical mechanical parameters combined with the elastic-plastic assumptions, hereby subjacent to the concept of strength. Second, the FE analyses in this study use a static loading scenario, whereas fatigue effects from repetitive or dynamic loads shall influence fracture risk through micro-damage accumulation (Mirzaali et al., 2020), especially when the bone turnover balance is altered, as in osteoporosis. Fatigue-related effects cannot be addressed by using bone strength alone, but we hypothesized that their representation can be improved through parameters such as residual displacement, nonlinear deformation, and energy absorption. Such assessment of the FE-predicted mechanical parameters, in the present study, considers cross-validation, which is otherwise too often overlooked in the literature.

To assess the discrimination power and clinical applicability of the proposed mechanical descriptors, we compare those with existing FE-based methods, including QCT-based, 2D-DXA-based, and 3D-DXA-based models. The discrimination capability of QCT-based FE analysis for hip fractures was first investigated in the MrOS study by Orwoll et al. (2009), which revealed that femoral strength (AUROC:0.83) and the load-to-strength ratio (AUROC:0.79) were key predictors, showing comparable performance with aBMD (AUROC:0.85), without any significant difference. Subsequent studies by Amin et al. (2011) in community-dwelling adults and Kopperdahl et al. (2014) used the AGES Reykjavik cohort and confirmed these findings, showing no significant improvement using FE-based parameters compared to aBMD. In contrast, our study demonstrates that the proposed FE strength + residual displacement descriptor from the 3D-DXA-based model provides superior fracture discrimination compared to aBMD and outperforms previous QCT-based FE studies (Orwoll et al., 2009; Amin et al., 2011; Kopperdahl et al., 2014). A large retrospective study by Adams et al. (2018), following the MrOS (Orwoll et al., 2009) and AGES Reykjavik (Kopperdahl et al., 2014) protocols, provided the first evidence of a significant advantage for QCT-based FE analysis (hazard ratio (HR) = 2.8, 95% CI 2.2–3.5) over aBMD (HR = 2.1, 95% CI 1.7–2.5), validating the clinical value of the Biomechanical Computed Tomography (BCT) approach. Keaveny et al. (2023) introduced BCT+, an advanced method incorporating nine variables, including bone strength, age, bone volume, fat content, and body composition. BCT+ outperformed BMD and FRAX, achieving an AUROC of 0.79 vs. 0.72 and 0.73, respectively. Our study showed a greater improvement in incident hip fracture assessment (ΔAUROC = 0.11 vs. 0.07), and importantly, our 3D-DXA-based FE models are verified against QCT-based FE models (Dudle et al., 2023; Qasim et al., 2024). Bessho et al. (2007) introduced a novel concept, later expanded by Falcinelli et al. (2014), who assessed femoral strength under various loading conditions to identify the weakest scenario for each patient. Their findings showed that minimum femoral strength better-discriminated fractures than standard sideways fall strength, achieving a higher AUROC than aBMD (0.88 vs. 0.79). However, the study had limitations, as QCT images were acquired post-fracture by scanning the contralateral femur, potentially affecting FE model accuracy in fracture risk assessment. Fleps et al. (2022) demonstrated that nonlinear FE models improved fracture prediction over aBMD (AUROC: 0.78 versus 0.72 in men, 0.75 versus 0.69 in women) within the AGES Reykjavik cohort. In their study, the most discriminative fall direction had low internal rotation and adduction, differing from Falcinelli et al. (2014). The AUROC values in our study were slightly higher than those reported by Fleps et al. (2022), further emphasizing the clinical potential of the 3D-DXA-based FE model alternative. Although QCT-based FE models show strong fracture discrimination (Adams et al., 2018; Keaveny et al., 2023), their high cost and radiation exposure limit widespread use. In addition, while most previous QCT-based FE studies focused on single metrics such as femoral strength, the load-to-strength ratio, or the multivariate BCT+ combining strength with clinical parameters, none have explored novel mechanical descriptors from FE simulations. Our work introduces descriptors integrating strength with nonlinear behaviour (e.g., residual displacement), capturing aspects of bone fragility beyond traditional measures. By rigorously assessing their discrimination performance, our findings extend existing QCT-based approaches and highlight the potential of global mechanical descriptors to improve fracture-risk prediction.

As an alternative, DXA-based FE models have been proposed to estimate femoral strength for hip fracture assessment (Wills et al., 2019; Grassi et al., 2023). DXA-based 2D FE models have been developed to assess bone strength and hip fracture risk. Naylor et al. (2013) found only marginal improvement in fracture discrimination using strength (AUROC = 0.67) and load-to-strength ratio (AUROC = 0.68) compared to FN aBMD (AUROC = 0.66), with a slight increase after adjustment for FN aBMD (AUROC = 0.68 and 0.69, for FN aBMD + strength and FN aBMD + load-to-strength ratio, respectively). Yang et al. (2014) reported femoral strength as a strong predictor of hip fracture risk after adjusting for age and BMI (HR = 2.2, 95% CI 1.95–2.50), outperforming TH aBMD (HR 1.9, 95% CI 1.7–2.1) but not FN aBMD (HR 2, 95% CI 1.8–2.3). More recently, DXA-derived indexes, such as the Bone Strain Index (BSI), have been developed. Sornay-Rendu et al. (2022) utilized a DXA-based 2D FE model to assess FN and TH BSI, alongside aBMD, in the OFELY study, which focuses on all major osteoporotic fractures instead of purely on incident hip fractures like in the present study. The AUROC for TH BSI (0.64) was not significantly different from that of TH aBMD (0.65) or Trabecular Bone Score (TBS) (0.67). In the current study, the discrimination power of aBMD, with an AUROC of 0.69, aligns with those results, and 3D-DXA vBMD tends to improve the discrimination (AUROC = 0.72). These outcomes reveal the likely importance of 3D information, especially if the latter is extended to mechanical descriptors (Table 4).

Indeed, since 2D models do not capture variations in bone geometry, density distribution, or impact forces, DXA-based 3D modeling has been developed to enable 3D-DXA-based FE models, as utilized in the current study (Wills et al., 2019; Grassi et al., 2023). Grassi et al. (2023) employed a 3D-DXA-based linear FE model in a MrOS cohort, demonstrating a significant improvement in the AUROC values for bone strength normalized by body weight (BW) across 10 different fall scenarios compared to BMD combined with BMI (0.78 vs. 0.72). This study is the only 3D-DXA-based investigation directly comparable to our approach and the first to show a statistically significant improvement in hip fracture discrimination over TH aBMD using subject-specific 2D-to-3D FE models. The current study exhibits slightly improved discrimination power compared to Grassi et al. (2023) (ΔAUROC = 0.11 vs. 0.06), likely due to key methodological differences. In contrast to the linear elastic approach adopted by the Authors (Grassi et al., 2023), we incorporated nonlinear material properties, providing a computationally efficient way to reflect the inelastic femur mechanics under side-fall loading, consistent with experimental findings 7,9. Fleps et al. 14 also demonstrated the superiority of nonlinear FE models in fracture discrimination. Besides, (Grassi et al., 2023) estimated femoral strength using a surface-strain criterion (Schileo et al., 2014; Schileo et al., 2008b) originally established from experimental observations (Grassi et al., 2012). While this approach restricts failure assessment to the bone surface, our methodology captures the volumetric behavior of the femur, not excluding that failure prediction might be controlled by the mechanics of internal compartments of the model (Wills et al., 2019). This provides a more comprehensive representation of bone mechanics, including the mechanics of the entire neck and intertrochanteric regions, facilitates clinical translation, and is consistent with the framework proposed by Keaveny et al. (2023), Orwoll et al. (2009), Amin et al. (2011).

Furthermore, our mechanical descriptors enhance patient-specific assessments by integrating organ geometry, BMD distribution, mechanical properties, and structural characteristics. In particular, residual displacement, nonlinear deformation, and energy absorption seemed to capture key determinants of bone fragility between the two patient groups, beyond bone strength alone, ultimately improving fracture prediction.

Extending the above observations, the effect of different fall scenarios was evaluated in the present study, and our findings are consistent with those of Fleps et al. (2022), who reported that fall models with low internal rotation and adduction are the most suitable for discriminating hip fracture cases from controls. Remarkably, the discrimination achieved by incorporating descriptors of the inelastic behaviors of the femur models was consistently superior to that achieved with only the state-of-the-art Strength parameter. Nevertheless, we found a high variability as per the performance of the mechanical descriptor to discriminate the fracture cases, versus the BMD parameters, when load directions were varied (Supplementary Appendix E, Supplementary Table SE1). Therefore, we recommend that probabilistic fracture risk be assessed across a range of fall directions (Falcinelli et al., 2014; Grassi et al., 2023).

The proposed 3D-DXA-based approach offers a clinically viable solution for fracture risk assessment by utilizing standard DXA scans, eliminating additional radiation exposure, and reducing costs compared to QCT-based FE methods (Fleps et al., 2022; Adams et al., 2018). Its computational efficiency, implemented using the commercially available 3D-Shaper software (Humbert et al., 2016), further enhances its clinical feasibility. Processing time per subject to compute 3D-DXA analysis using 3D-Shaper software is 1 min 30 s on an Intel Core i7-4790K CPU, which is significantly faster than the methods used by Grassi et al. (2023), which required approximately 6 hours per scan using 20 CPU cores. Additionally, our FE simulation takes about 15 min per subject using 4 CPU cores, compared to 2 hours reported by Fleps et al. (2022). The existing regulatory approvals for the 3D-Shaper software are also an asset for the future clinical implementation of the 3D-DXA-based FE approach evaluated in this study.

Despite its advantages, the proposed approach has some limitations. The study was based on a retrospective cohort of 64 hip fracture cases and 64 controls. Larger, prospective cohorts are needed for further validation. The elastic-plastic model used to capture inelastic bone mechanics was a rough approximation of the real bone failure mechanics. It allowed computationally affordable simulations and indicated a robust improvement compared to BMD measure or strength calculations alone. Yet, whether macroscopic or multiscale bone failure mechanics models can further improve fracture risk predictions remains to be explored. A challenge also arises in femoral head density reconstruction, as clinical DXA scans often suffer from acetabulum superimposition. This region is excluded from the 2D-to-3D reconstruction using 3D-Shaper software (Humbert et al., 2016), and the femoral head is reconstructed by mapping QCT-derived density onto the DXA image (Humbert et al., 2016). This approximation may affect the accuracy of the femoral geometry and material properties, but rigid-body constraints with a reference point were applied during the FE simulation to mitigate the effect, prevent stress concentration in the femur head, and lead to a more realistic model. The present study prioritized global over local mechanical descriptors (Wills et al., 2019). While this can be considered as a limitation, as local descriptors might achieve improved discriminations (Wills et al., 2019), we hereby avoid the need for site-specific fracture analyses. This strategy reduces sensitivity to mesh convergence, enhances computational efficiency, and reduces the complexity of the mechanical analysis (only one surrogate descriptor for bone fracture risk at the organ level). It gathers, therefore, important characteristics to make the clinical translation accessible. Indeed, the current cohort does not include annotations of fracture location (neck vs. trochanteric). Hence, the possibility to adopt global descriptors fits well with the reality of the data.

5 Conclusion

This case-control study demonstrates a significant improvement in hip fracture discrimination using 3D-DXA-based FE models compared to standard aBMD and 3D-DXA measurements. The consideration of the global balance of elastic and inelastic femur mechanics during side-fall simulations was a cornerstone to achieve improved fragile fracture discrimination at reduced computational cost. These findings suggest that 3D-DXA-based FE models could enhance osteoporosis management while maintaining DXA as the clinical standard of care. Among the mechanical descriptors evaluated, the combination of FE strength and residual displacement emerged as the most effective for discriminating hip fractures, as it incorporates key mechanical parameters relevant to hip fracture risk.

Data availability statement

The data analyzed in this study is subject to the following licenses/restrictions: The data are not publicly available due to privacy or ethical restrictions. Requests to access these datasets should be directed to Ludovic Humbert, bHVkb3ZpYy5odW1iZXJ0QDNkLXNoYXBlci5jb20=.

Ethics statement

The studies involving humans were approved by CETIR Group Mèdic Scientific Committee. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

EA: Conceptualization, Formal Analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review and editing. CR: Conceptualization, Methodology, Writing – review and editing, Validation. LR: Data curation, Writing – review and editing. SD: Data curation, Writing – review and editing. LH: Conceptualization, Investigation, Methodology, Software, Supervision, Validation, Writing – review and editing. JN: Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was funded by the Spanish Government (CPP2021-008574 and CPP2023-010970) and AGAUR (2023 BP 00224).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2025.1673339/full#supplementary-material

References

Abe, S., Kouhia, R., Nikander, R., Narra, N., Hyttinen, J., and Sievänen, H. (2022). Effect of fall direction on the lower hip fracture risk in athletes with different loading histories: a finite element modeling study in multiple sideways fall configurations. Bone 158, 116351. doi:10.1016/j.bone.2022.116351

PubMed Abstract | CrossRef Full Text | Google Scholar

Adams, A. L., Fischer, H., Kopperdahl, D. L., Lee, D. C., Black, D. M., Bouxsein, M. L., et al. (2018). Osteoporosis and hip fracture risk from routine computed tomography scans: the fracture, osteoporosis, and CT utilization study (FOCUS). J. Bone Min. Res. 33, 1291–1301. doi:10.1002/jbmr.3423

PubMed Abstract | CrossRef Full Text | Google Scholar

Adler, R. A. (2014). Osteoporosis in men: a review. Bone Res. 2, 14001–14008. doi:10.1038/boneres.2014.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Airola, A., Pahikkala, T., Waegeman, W., De Baets, B., and Salakoski, T. (2011). An experimental comparison of cross-validation techniques for estimating the area under the ROC curve. Comput. Stat. \and Data Anal. 55, 1828–1844. doi:10.1016/j.csda.2010.11.018

CrossRef Full Text | Google Scholar

Amin, S., Kopperdhal, D. L., Melton, L. J., Achenbach, S. J., Therneau, T. M., Riggs, B. L., et al. (2011). Association of hip strength estimates by finite-element analysis with fractures in women and men. J. bone Min. Res. 26, 1593–1600. doi:10.1002/jbmr.347

PubMed Abstract | CrossRef Full Text | Google Scholar

Bayraktar, H. H., Morgan, E. F., Niebur, G. L., Morris, G. E., Wong, E. K., and Keaveny, T. M. (2004). Comparison of the elastic and yield properties of human femoral trabecular and cortical bone tissue. J. Biomech. 37, 27–35. doi:10.1016/s0021-9290(03)00257-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Bessho, M., Ohnishi, I., Matsuyama, J., Matsumoto, T., Imai, K., and Nakamura, K. (2007). Prediction of strength and strain of the proximal femur by a CT-based finite element method. J. Biomech. 40, 1745–1753. doi:10.1016/j.jbiomech.2006.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Bevill, G., Easley, S. K., and Keaveny, T. M. (2007). Side-artifact errors in yield strength and elastic modulus for human trabecular bone and their dependence on bone volume fraction and anatomic site. J. Biomech. 40, 3381–3388. doi:10.1016/j.jbiomech.2007.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, X., Keyak, J. H., Sigurdsson, S., Zhao, C., Zhou, W., Liu, A., et al. (2024). A new hip fracture risk index derived from FEA-Computed proximal femur fracture loads and energies-to-failure. Osteoporos. Int. 1–10. doi:10.1007/s00198-024-07015-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Cody, D. D., Gross, G. J., J. Hou, F., Spencer, H. J., Goldstein, S. A., and P. Fyhrie, D. (1999). Femoral strength is better predicted by finite element models than QCT and DXA. J. Biomech. 32, 1013–1020. doi:10.1016/s0021-9290(99)00099-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Courtney, A. C., Wachtel, E. F., Myers, E. R., and Hayes, W. C. (1995). Age-related reductions in the strength of the femur tested in a fall-loading configuration. JBJS 77, 387–395. doi:10.2106/00004623-199503000-00008

PubMed Abstract | CrossRef Full Text | Google Scholar

Dall’Ara, E., Eastell, R., Viceconti, M., Pahr, D., and Yang, L. (2016). Experimental validation of DXA-Based finite element models for prediction of femoral strength. J. Mech. Behav. Biomed. Mat. 63, 17–25. doi:10.1016/j.jmbbm.2016.06.004

CrossRef Full Text | Google Scholar

DeLong, E. R., DeLong, D. M., and Clarke-Pearson, D. L. (1988). Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics 44, 837–845. doi:10.2307/2531595

PubMed Abstract | CrossRef Full Text | Google Scholar

Dragomir-Daescu, D., Op Den Buijs, J., McEligot, S., Dai, Y., Entwistle, R. C., Salas, C., et al. (2011). Robust QCT/FEA models of proximal femur stiffness and fracture load during a sideways fall on the hip. Ann. Biomed. Eng. 39, 742–755. doi:10.1007/s10439-010-0196-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Dudle, A., Gugler, Y., Pretterklieber, M., Ferrari, S., Lippuner, K., and Zysset, P. (2023). 2D-3D reconstruction of the proximal femur from DXA scans: evaluation of the 3D-Shaper software. Front. Bioeng. Biotechnol. 11, 1111020. doi:10.3389/fbioe.2023.1111020

PubMed Abstract | CrossRef Full Text | Google Scholar

Enns-Bray, W. S., Bahaloo, H., Fleps, I., Pauchard, Y., Taghizadeh, E., Sigurdsson, S., et al. (2019). Biofidelic finite element models for accurately classifying hip fracture in a retrospective clinical study of elderly women from the AGES reykjavik cohort. Bone 120, 25–37. doi:10.1016/j.bone.2018.09.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Falcinelli, C., Schileo, E., Balistreri, L., Baruffaldi, F., Bordini, B., Viceconti, M., et al. (2014). Multiple loading conditions analysis can improve the association between finite element bone strength estimates and proximal femur fractures: a preliminary study in elderly women. Bone 67, 71–80. doi:10.1016/j.bone.2014.06.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Fleps, I., and Morgan, E. F. (2022). A review of CT-based fracture risk assessment with finite element modeling and machine learning. Curr. Osteoporos. Rep. 20, 309–319. doi:10.1007/s11914-022-00743-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Fleps, I., Pálsson, H., Baker, A., Enns-Bray, W., Bahaloo, H., Danner, M., et al. (2022). Finite element derived femoral strength is a better predictor of hip fracture risk than aBMD in the AGES reykjavik study cohort. Bone 154, 116219. doi:10.1016/j.bone.2021.116219

PubMed Abstract | CrossRef Full Text | Google Scholar

Goulet, Rw., Goldstein, S., Ciarelli, M., Kuhn, J., Brown, M., and Feldkamp, L. (1994). The relationship between the structural and orthogonal compressive properties of trabecular bone. J. Biomech. 27, 375–389. doi:10.1016/0021-9290(94)90014-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Grassi, L., Väänänen, S. P., Jehpsson, L., Ljunggren, Ö., Rosengren, B. E., Karlsson, M. K., et al. (2023). 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. J. Bone Min. Res. 38, 1258–1267. doi:10.1002/jbmr.4878

PubMed Abstract | CrossRef Full Text | Google Scholar

Grassi, L., Schileo, E., Taddei, F., Zani, L., Juszczyk, M., Cristofolini, L., et al. (2012). Accuracy of finite element predictions in sideways load configurations for the proximal human femur. J. Biomech. 45, 394–399. doi:10.1016/j.jbiomech.2011.10.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Grassi, L., Väänänen, S. P., Ristinmaa, M., Jurvelin, J. S., and Isaksson, H. (2017). Prediction of femoral strength using 3D finite element models reconstructed from DXA images: validation against experiments. Biomech. Model. Mechanobiol. 16, 989–1000. doi:10.1007/s10237-016-0866-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Grassi, L., Fleps, I., Sahlstedt, H., Väänänen, S. P., Ferguson, S. J., Isaksson, H., et al. (2021). Validation of 3D finite element models from simulated DXA images for biofidelic simulations of sideways fall impact to the hip. Bone 142, 115678. doi:10.1016/j.bone.2020.115678

PubMed Abstract | CrossRef Full Text | Google Scholar

Grassi, L., Väänänen, S. P., Voss, A., Nissinen, T., Sund, R., Kröger, H., et al. (2025). DXA-Based 3D finite element models predict hip fractures better THAN areal BMD in elderly women. Bone 195, 117457. doi:10.1016/j.bone.2025.117457

PubMed Abstract | CrossRef Full Text | Google Scholar

Hodgskinson, R., and Currey, J. D. (1992). Young’s modulus, density and material properties in cancellous bone over a large density range. J. Mat. Sci. Mat. Med. 3, 377–381. doi:10.1007/bf00705371

CrossRef Full Text | Google Scholar

Humbert, L., Martelli, Y., Fonolla, R., Steghofer, M., Di Gregorio, S., Malouf, J., et al. (2016). 3D-DXA: assessing the femoral shape, the trabecular macrostructure and the cortex in 3D from DXA images. IEEE Trans. Med. Imaging 36, 27–39. doi:10.1109/tmi.2016.2593346

PubMed Abstract | CrossRef Full Text | Google Scholar

Humbert, L., Bagué, A., Di Gregorio, S., Winzenrieth, R., Sevillano, X., González Ballester, M. Á., et al. (2020). DXA-Based 3D analysis of the cortical and trabecular bone of hip fracture postmenopausal women: a case-control study. J. Clin. Densitom. 23, 403–410. doi:10.1016/j.jocd.2018.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Johannesdottir, F., Thrall, E., Muller, J., Keaveny, T. M., Kopperdahl, D. L., and Bouxsein, M. L. (2017). Comparison of non-invasive assessments of strength of the proximal femur. Bone 105, 93–102. doi:10.1016/j.bone.2017.07.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanis, J. A., Norton, N., Harvey, N. C., Jacobson, T., Johansson, H., Lorentzon, M., et al. (2021). SCOPE 2021: a new scorecard for osteoporosis in Europe. Arch. Osteoporos. 16, 82. doi:10.1007/s11657-020-00871-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Keaveny, D. L., Bouxsein, M., Orwoll, E., Sundeep Khosla, M., McClung, E. S., and Kopperdahl, D. (2023). Improved real-world prediction of hip fracture using CT-Based measurements of bone, muscle, and soft tissue in women and men. J. Bone Min. Res. 38, 1545–2024.

Google Scholar

Keller, T. S. (1994). Predicting the compressive mechanical behavior of bone. J. Biomech. 27, 1159–1168. doi:10.1016/0021-9290(94)90056-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Keyak, J. H., Lee, I. Y., and Skinner, H. B. (1994). Correlations between orthogonal mechanical properties and density of trabecular bone: use of different densitometric measures. J. Biomed. Mat. Res. 28, 1329–1336. doi:10.1002/jbm.820281111

PubMed Abstract | CrossRef Full Text | Google Scholar

Kok, J., Grassi, L., Gustafsson, A., and Isaksson, H. (2021). Femoral strength and strains in sideways fall: validation of finite element models against bilateral strain measurements. J. Biomech. 122, 110445. doi:10.1016/j.jbiomech.2021.110445

PubMed Abstract | CrossRef Full Text | Google Scholar

Kopperdahl, D. L., Aspelund, T., Hoffmann, P. F., Sigurdsson, S., Siggeirsdottir, K., Harris, T. B., et al. (2014). Assessment of incident spine and hip fractures in women and men using finite element analysis of CT scans. J. Bone Min. Res. 29, 570–580. doi:10.1002/jbmr.2069

PubMed Abstract | CrossRef Full Text | Google Scholar

Leslie, W. D., Luo, Y., Yang, S., Goertzen, A. L., Ahmed, S., Delubac, I., et al. (2019). Fracture risk indices from DXA-based finite element analysis predict incident fractures independently from FRAX: the Manitoba BMD Registry. J. Clin. Densitom. 22, 338–345. doi:10.1016/j.jocd.2019.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

McCloskey, E., Johansson, H., Oden, A., and Kanis, J. A. (2012). Fracture risk assessment. Clin. Biochem. 45, 887–893. doi:10.1016/j.clinbiochem.2012.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Mirzaali, M. J., Libonati, F., Böhm, C., Rinaudo, L., Cesana, B., Ulivieri, F., et al. (2020). Fatigue-caused damage in trabecular bone from clinical, morphological and mechanical perspectives. Int. J. Fatigue 133, 105451. doi:10.1016/j.ijfatigue.2019.105451

CrossRef Full Text | Google Scholar

Morgan, E. F., and Keaveny, T. M. (2001). Dependence of yield strain of human trabecular bone on anatomic site. J. Biomech. 34, 569–577. doi:10.1016/s0021-9290(01)00011-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Naylor, K. E., McCloskey, E. V., Eastell, R., and Yang, L. (2013). Use of DXA-based finite element analysis of the proximal femur in a longitudinal study of hip fracture. J. Bone Min. Res. 28, 1014–1021. doi:10.1002/jbmr.1856

PubMed Abstract | CrossRef Full Text | Google Scholar

Orwoll, E. S., Marshall, L. M., Nielson, C. M., Cummings, S. R., Lapidus, J., Cauley, J. A., et al. (2009). Finite element analysis of the proximal femur and hip fracture risk in older men. J. bone Min. Res. 24, 475–483. doi:10.1359/jbmr.081201

PubMed Abstract | CrossRef Full Text | Google Scholar

Pinilla, T. P., Boardman, K. C., Bouxsein, M. L., Myers, E. R., and Hayes, W. C. (1996). Impact direction from a fall influences the failure load of the proximal femur as much as age-related bone loss. Calcif. Tissue Int. 58, 231–235. doi:10.1007/bf02508641

PubMed Abstract | CrossRef Full Text | Google Scholar

Qasim, M., López Picazo, M., Ruiz Wills, C., Noailly, J., Di Gregorio, S., Del Río Barquero, L. M., et al. (2024). 3D-DXA based finite element modelling for femur strength prediction: evaluation against QCT. J. Clin. Densitom. 27, 101471. doi:10.1016/j.jocd.2024.101471

PubMed Abstract | CrossRef Full Text | Google Scholar

Rinaudo, L., Cuttone, S., Messina, C., Magni, V., Capra, D., Sconfienza, L. M., et al. (2024). Generation and validation of normative, age-specific reference curves for bone strain index in women. Diagnostics 14, 1046. doi:10.3390/diagnostics14101046

PubMed Abstract | CrossRef Full Text | Google Scholar

Schileo, E., and Taddei, F. (2021). Finite element assessment of bone fragility from clinical images. Curr. Osteoporos. Rep. 19, 688–698. doi:10.1007/s11914-021-00714-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Schileo, E., Dall’Ara, E., Taddei, F., Malandrino, A., Schotkamp, T., Baleani, M., et al. (2008a). An accurate estimation of bone density improves the accuracy of subject-specific finite element models. J. Biomech. 41, 2483–2491. doi:10.1016/j.jbiomech.2008.05.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Schileo, E., Taddei, F., Cristofolini, L., and Viceconti, M. (2008b). Subject-specific finite element models implementing a maximum principal strain criterion are able to estimate failure risk and fracture location on human femurs tested in vitro. J. Biomech. 41, 356–367. doi:10.1016/j.jbiomech.2007.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Schileo, E., Balistreri, L., Grassi, L., Cristofolini, L., and Taddei, F. (2014). To what extent can linear finite element models of human femora predict failure under stance and fall loading configurations? J. Biomech. 47, 3531–3538. doi:10.1016/j.jbiomech.2014.08.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Sornay-Rendu, E., Duboeuf, F., Ulivieri, F. M., Rinaudo, L., and Chapurlat, R. (2022). The bone strain index predicts fragility fractures. The OFELY study. Bone 157, 116348. doi:10.1016/j.bone.2022.116348

PubMed Abstract | CrossRef Full Text | Google Scholar

Stone, M. (1974). Cross-validatory choice and assessment of statistical predictions. J. R. Stat. Soc. Ser. B 36, 111–133. doi:10.1111/j.2517-6161.1974.tb00994.x

CrossRef Full Text | Google Scholar

Stone, K. L., Seeley, D. G., Lui, L. Y., Cauley, J. A., Ensrud, K., Browner, W. S., et al. (2003). BMD at multiple sites and risk of fracture of multiple types: long-term results from the Study of Osteoporotic Fractures. J. bone Min. Res. 18, 1947–1954. doi:10.1359/jbmr.2003.18.11.1947

PubMed Abstract | CrossRef Full Text | Google Scholar

Terzini, M., Aldieri, A., Rinaudo, L., Osella, G., Audenino, A. L., and Bignardi, C. (2019). Improving the hip fracture risk prediction through 2D finite element models from DXA images: validation against 3D models. Front. Bioeng. Biotechnol. 7, 220. doi:10.3389/fbioe.2019.00220

PubMed Abstract | CrossRef Full Text | Google Scholar

Thevenot, J., Koivumäki, J., Kuhn, V., Eckstein, F., and Jämsä, T. (2014). A novel methodology for generating 3D finite element models of the hip from 2D radiographs. J. Biomech. 47, 438–444. doi:10.1016/j.jbiomech.2013.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Ulivieri, F. M., and Rinaudo, L. (2022). The Bone Strain Index: an innovative dual X-ray absorptiometry bone strength index and its helpfulness in clinical medicine. J. Clin. Med. 11, 2284. doi:10.3390/jcm11092284

PubMed Abstract | CrossRef Full Text | Google Scholar

Ulivieri, F. M., Rinaudo, L., Messina, C., Aliprandi, A., Sconfienza, L. M., Sardanelli, F., et al. (2022). Bone Strain Index: preliminary distributional characteristics in a population of women with normal bone mass, osteopenia and osteoporosis. Radiol. Med. 127, 1151–1158. doi:10.1007/s11547-022-01543-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Väänänen, S. P., Grassi, L., Flivik, G., Jurvelin, J. S., and Isaksson, H. (2015). Generation of 3D shape, density, cortical thickness and finite element mesh of proximal femur from a DXA image. Med. Image Anal. 24, 125–134. doi:10.1016/j.media.2015.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Viceconti, M., Qasim, M., Bhattacharya, P., and Li, X. (2018). Are CT-based finite element model predictions of femoral bone strengthening clinically useful? Curr. Osteoporos. Rep. 16, 216–223. doi:10.1007/s11914-018-0438-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Watts, N. B., Ettinger, B., and LeBoff, M. S. (2009). FRAX facts. J. Bone Min. Res. 24, 975–979. doi:10.1359/jbmr.090402

PubMed Abstract | CrossRef Full Text | Google Scholar

Wills, C. R., Olivares, A. L., Tassani, S., Ceresa, M., Zimmer, V., González Ballester, M. A., et al. (2019). 3D patient-specific finite element models of the proximal femur based on DXA towards the classification of fracture and non-fracture cases. Bone 121, 89–99. doi:10.1016/j.bone.2019.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, L., Palermo, L., Black, D. M., and Eastell, R. (2014). Prediction of incident hip fracture with the estimated femoral strength by finite element analysis of DXA scans in the study of osteoporotic fractures. J. Bone Min. Res. 29, 2594–2600. doi:10.1002/jbmr.2291

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, S., Leslie, W. D., Luo, Y., Goertzen, A. L., Ahmed, S., Ward, L. M., et al. (2018). Automated DXA-based finite element analysis for hip fracture risk stratification: a cross-sectional study. Osteoporos. Int. 29, 191–200. doi:10.1007/s00198-017-4232-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Zysset, P., Pahr, D., Engelke, K., Genant, H. K., McClung, M. R., Kendler, D. L., et al. (2015). Comparison of proximal femur and vertebral body strength improvements in the FREEDOM trial using an alternative finite element methodology. Bone 81, 122–130. doi:10.1016/j.bone.2015.06.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: biomechanics, osteoporosis, hip fracture discrimination, 3D-DXA, finite element

Citation: Alizadeh E, Ruiz Wills C, del Rio L, Di Gregorio S, Humbert L and Noailly J (2025) Inelastic mechanical descriptors for osteoporotic hip fracture discrimination with 3D-DXA-based nonlinear finite element models. Front. Bioeng. Biotechnol. 13:1673339. doi: 10.3389/fbioe.2025.1673339

Received: 25 July 2025; Accepted: 30 September 2025;
Published: 23 October 2025.

Edited by:

Sandipan Roy, SRM Institute of Science and Technology, India

Reviewed by:

Jayanta Kumar Biswas, National Institute of Technology Patna, India
P.Susai Manickam, SRM University, India

Copyright © 2025 Alizadeh, Ruiz Wills, del Rio, Di Gregorio, Humbert and Noailly. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Elham Alizadeh, ZWxoYW0uYWxpemFkZWhAdXBmLmVkdQ==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.