Relationships Between Femoral Strength Evaluated by Nonlinear Finite Element Analysis and BMD, Material Distribution and Geometric Morphology

QCT-based 3D finite element modeling to assess patient-specific hip fracture risk and risk factors

Early assessment of hip fracture risk may play a critical role in designing preventive mechanisms to reduce the occurrence of hip fracture in geriatric people. The loading direction, clinical, and morphological variables play a vital role in hip fracture. Analyzing the effects of these variables helps predict fractures risk more accurately; thereby suggesting the critical variable that needs to be considered. Hence, this work considered the fall postures by varying the loading direction on the coronal plane (α) and on the transverse plane (β) along with the clinical variables-age, sex, weight, and bone mineral density, and morphological variables-femoral neck axis length, femoral neck width, femoral neck angle, and true moment arm. The strain distribution obtained via finite element analysis (FEA) shows that the angle of adduction (α) during a fall increases the risk of fracture at the greater trochanter and femoral neck, whereas with an increased angle of rotation (β) during the fall, the FRI increases by ∼1.35 folds. The statistical analysis of clinical, morphological, and loading variables (αandβ) delineates that the consideration of only one variable is not enough to realistically predict the possibility of fracture as the correlation between individual variables and FRI is less than 0.1, even though they are shown to be significant (p<0.01). On the contrary, the correlation (R2=0.48) increases as all variables are considered, suggesting the need for considering different variables fork predicting FRI. However, the effect of each variable is different. While loading, clinical, and morphological variables are considered together, the loading direction on transverse plane (β) has high significance, and the anatomical variabilities have no significance.

Biomechanical analysis analyzing association between bone mineral density and lag screw migration

A proximal femoral nail using a helical blade (HB) is commonly utilized to treat proximal femoral fracture but cut through failure of the lag screws is one of the devastating complications following the surgery. While controversial, one of the potential risk factors for cut through failure is poor bone strength which can be predicted by measuring bone mineral density (BMD). In this study, we performed a biomechanical test on the fractured femoral head to validate whether the indirectly measured BMD from the contralateral hip or that measured directly from the retrieved femoral head can elucidate the structural strength of the fractured femoral head and thereby can be used to predict migration of lag screws. Our result showed that directly measured BMD has a significant correlation with the HB migration on the osteoporotic femoral head. However, while the BMDs measured from the contralateral femoral neck or total hip is the most widely used parameter to predict the bone strength of the fractured femur, this may have limited usability to predict HB migration.

Prediction of femoral strength of elderly men based on quantitative computed tomography images using machine learning

Hip fracture is the most common complication of osteoporosis, and its major contributor is compromised femoral strength. This study aimed to develop practical machine learning models based on clinical quantitative computed tomography (QCT) images for predicting proximal femoral strength. Eighty subjects with entire QCT data of the right hip region were randomly selected from the full MrOS cohorts, and their proximal femoral strengths were calculated by QCT-based finite element analysis (QCT/FEA). A total of 50 parameters of each femur were extracted from QCT images as the candidate predictors of femoral strength, including grayscale distribution, regional cortical bone mapping (CBM) measurements, and geometric parameters. These parameters were simplified by using feature selection and dimensionality reduction. Support vector regression (SVR) was used as the machine learning algorithm to develop the prediction models, and the performance of each SVR model was quantified by the mean squared error (MSE), the coefficient of determination (R² ), the mean bias, and the SD of bias. For feature selection, the best prediction performance of SVR models was achieved by integrating the grayscale value of 30% percentile and specific regional CBM measurements (MSE ≤ 0.016, R² ≥ 0.93); and for dimensionality reduction, the best prediction performance of SVR models was achieved by extracting principal components with eigenvalues greater than 1.0 (MSE ≤ 0.014, R² ≥ 0.93). The femoral strengths predicted from the well-trained SVR models were in good agreement with those derived from QCT/FEA. This study provided effective machine learning models for femoral strength prediction, and they may have great potential in clinical bone health assessments.

Femur geometry and body composition influence femoral neck stresses: A combined fall simulation and beam modelling approach

Metrics of femur geometry and body composition have been linked to clinical hip fracture risk. Mechanistic explanations for these relationships have generally focused on femur strength; however, impact loading also modulates fracture risk. We evaluated the potential effects of femur geometry and body composition on femoral neck stresses during lateral impacts. Fifteen female volunteers completed low-energy sideways falls on to the hip. Additionally, participants completed ultrasound and dual-energy x-ray absorptiometry imaging to characterize trochanteric soft tissue thickness (TSTT) over the hip and six metrics of femur geometry, respectively. Subject-specific beam models were developed and utilized to calculate peak femoral neck stress (σ_Neck), utilizing experimental impact dynamics. Except for femoral neck axis length, all metrics of femur geometry were positively correlated with σ_Neck (all p < 0.05). Larger/more prominent proximal femurs were associated with increased force over the proximal femur, whereas a wider neck-shaft angle was associated with greater stress generation independent of force (all p < 0.05). Body mass index (BMI) and TSTT were negatively correlated with σ_Neck (both p < 0.05). Despite strong correlations, these metrics of body composition appear to influence femoral neck stresses through different mechanisms. Increased TSTT was associated with reduced force over the proximal femur, whereas increased BMI was associated with greater resistance to stress generation (both p < 0.05). This study provided novel insights into the mechanistic pathways through which femur geometry and body composition may modulate hip fracture risk. Our findings complement clinical findings and provide one possible explanation for incongruities in the clinical fracture risk and femur strength literature.

Varus malalignment of the lower limb increases the risk of femoral neck fracture: A biomechanical study using a finite element method

INTRODUCTION

The understanding of the stresses and strains and their dependence on loading direction caused by an axial deformity is very important for understanding the mechanism of femural neck fractures. The hypothesis of this study is that lower limb malalignment is correlated with a substantial stress variation on the upper end of the femur. The purpose of this biomechanical trial using the finite element method is to determine the effect of the loading direction on the proximal femur regarding the malalignment of the lower limb, and also enlighten the relation between the lower limb alignment and the risk of a femoral neck fracture.

METHODS

Ten segmentations of CT scans were considered. An axial compression load was applied to the femoral head to digitally simulate the physiological configuration in neutral position as well as in different axial positions in varus/valgus alignment.

RESULTS

The stress at the proximal femur changes as the varus _valgus angle does. It can be observed the smaller absolute stress at angle 10° (valgus) and the higher absolute stress at angle -10° (varus). The mean maximum von Mises stress value was 14.1 (SD=±3.48) MPa for 0°, while the mean maximum von Mises stress value was 17.96 MPa (SD=4.87) for -10° in varus. The fracture risk indicator of the proximal femoral epiphyses changes inversely with angle direction. The FRI was the highest at -10° and the lowest at 10°.

CONCLUSION

Based on the biomechanical findings and the fracture risk indicator determined in this preliminary study, varus malalignment increases the risk of femoral neck fracture. Consideration of other parameters such as bone mineral density and morphological parameters should also help to plan preventive medical strategy in the elderly.

A deep learning-based approach to automatic proximal femur segmentation in quantitative CT images

Automatic CT segmentation of proximal femur has a great potential for use in orthopedic diseases, especially in the imaging-based assessments of hip fracture risk. In this study, we proposed an approach based on deep learning for the fast and automatic extraction of the periosteal and endosteal contours of proximal femur in order to differentiate cortical and trabecular bone compartments. A three-dimensional (3D) end-to-end fully convolutional neural network (CNN), which can better combine the information among neighbor slices and get more accurate segmentation results by 3D CNN, was developed for our segmentation task. The separation of cortical and trabecular bones derived from the QCT software MIAF-Femur was used as the segmentation reference. Two models with the same network structures were trained, and they achieved a dice similarity coefficient (DSC) of 97.82% and 96.53% for the periosteal and endosteal contours, respectively. Compared with MIAF-Femur, it takes half an hour to segment a case, and our CNN model takes a few minutes. To verify the excellent performance of our model for proximal femoral segmentation, we measured the volumes of different parts of the proximal femur and compared it with the ground truth, and the relative errors of femur volume between predicted result and ground truth are all less than 5%. This approach will be expected helpful to measure the bone mineral densities of cortical and trabecular bones, and to evaluate the bone strength based on FEA.

Biomechanical evaluation of the unilateral crossbite on the asymmetrical development of the craniofacial complex. A mechano-morphological approach

BACKGROUND AND OBJECTIVE

The occlusion effect on the craniofacial development is a controversial topic that has attracted the interest of many researchers but that remains unclear, mainly due to the difficulties on measure its mechanical response experimentally. This mechano-morphological relationship of the craniofacial growth is often explained by the periosteal and capsular matrices of the functional matrix hypothesis (FMH); however, its outcomes have not been analytically demonstrated yet. This computational study aims, therefore, to analytically demonstrate the mechano-morphological relationship in the craniofacial development of children with unilateral crossbite (UXB) using the finite element (FE) method.

METHODS

The craniofacial complex asymmetry of ten children, five of whom exhibit UXB, was 3D-analysed and compared with the biomechanical response computed from a FE analysis of each patient's occlusion. Due to the complexity of the geometry and the multitude of contacts involved, the inherent limitations of the model were evaluated by comparing computed occlusal patterns with those recorded by an occlusal analysis on 3D printed copies.

RESULTS

Comparison's outcomes proved the reliability of our models with just a deviation error below 6% between both approaches. Out of validation process, computational results showed that the significant elongation of mandibular branch in the contralateral side could be related to the mandibular shift and increase of thickness on the crossed side, and particularly of the posterior region. These morphological changes could be associated with periodontal overpressure (>4.7 kPa) and mandibular over deformation (0.002 ε) in that side, in agreement with the periosteal matrix's principles. Furthermore, the maxilla's transversal narrowing and the elevation of the maxillary and zygomatic regions on the crossed side were statistically demonstrated and seem to be related with their respective micro displacements at occlusion, as accounted by their specific capsule matrices. Our results were consistent with those reported clinically and demonstrated analytically the mechano-morphological relationship of children's craniofacial development based on the FMH's functional matrices.

CONCLUSIONS

This study is a first step in the understanding of the occlusion's effect on the craniofacial development by computational methods. Our approach could help future engineers, researchers and clinicians to understand better the aetiology of some dental malocclusions and functional disorders improve the diagnosis or even predict the craniofacial development.

Prediction of Femoral Strength Based on Bone Density and Biochemical Markers in Elderly Men With Type 2 Diabetes Mellitus

Purpose: Effects of bone density, bone turnover and advanced glycation end products (AGEs) on femoral strength (FS) are still unclear in patients with type 2 diabetes mellitus (T2DM). This study aims to assess and predict femoral strength and its influencing factors in elderly men with T2DM.

Methods: T2DM patients (n = 10, mean age, 66.98 years) and age-matched controls (n = 8, mean age, 60.38 years) were recruited. Femoral bone mineral density (BMD) and serum biochemical indices of all subjects were measured. FS was evaluated through finite element analysis based on quantitative computed tomography. Multiple linear regression was performed to obtain the best predictive models of FS and to analyze the ability of predictors of FS in both groups.

Results: FS (p = 0.034), HbA1c (p = 0.000) and fasting blood glucose (p = 0.000) levels of T2DM group were significantly higher than those of control group; however, the P1NP level (p = 0.034) was significantly lower. FS was positively correlated with femoral neck T score (FNTS) (r = 0.794, p &lt; 0.01; r = 0.881, p &lt; 0.01) in both groups. FS was correlated with age (r = -0.750, p &lt; 0.05) and pentosidine (r = -0.673, p &lt; 0.05) in T2DM group. According to multiple linear regression, FNTS and P1NP both contributed to FS in two groups. P1NP significantly improved the prediction of FS in both groups, but significant effect of FNTS on predicting FS was only presented in control group. Furthermore, pentosidine, age and HbA1c all played significant roles in predicting FS of T2DM.

Conclusion: Femoral strength was higher in elderly men with T2DM, which might be caused by higher BMD and lower bone turnover rate. Moreover, besides BMD and bone formation level, AGEs, blood glucose and age might significantly impact the prediction of femoral strength in T2DM.

Biomechanical analysis of the correlation between mid-shaft atypical femoral fracture (AFF) and axial varus deformation

Background

Atypical femoral fractures (AFF) are diaphyseal fractures of the elderly that occur at the end of a minor trauma. The objective of this biomechanical study, using finite element modelling, was to evaluate the variations of the femoral diaphysis fracture indicator according to the variations of the mechanical axis of the lower limb, which can explain all the different atypical fracture types identified in the literature.

Methods

In order to measure variations in stress and risk factors for fracture of the femoral diaphysis, the distal end of the femur was constrained in all degrees of freedom. An axial compression load was applied to the femoral head to digitally simulate the bipodal support configuration in neutral position as well as in different axial positions in varus/valgus (− 10°/10°).

Results

The maximum stress value of Von Mises was twice as high (17.96 ± 4.87 MPa) at a varus angle of − 10^° as in the neutral position. The fracture risk indicator of the femoral diaphysis varies proportionally with the absolute value of the steering angle. However, the largest simulated varus deformation (− 10°) found a higher risk of diaphysis fracture indicator than in valgus (10°).

Conclusions

Variations in the mechanical axis of the lower limb influence the stress distribution at the femur diaphysis and consequently increase the risk of AFF. The axial deformation in varus is particularly at risk of AFF. The combination of axial deformation stresses and bone fragility consequently contribute to the creation of an environment favorable to the development of AFF.

Trial registration: ‘retrospectively registered’.

Study of the significance of parameters and their interaction on assessing femoral fracture risk by quantitative statistical analysis

Early assessment of hip fracture helps develop therapeutic and preventive mechanisms that may reduce the occurrence of hip fracture. An accurate assessment of hip fracture risk requires proper consideration of the loads, the physiological and morphological parameters, and the interactions between these parameters. Hence, this study aims at analyzing the significance of parameters and their interactions by conducting a quantitative statistical analysis. A multiple regression model was developed considering different loading directions during a sideways fall (angle [Formula: see text] and [Formula: see text] on the coronal and transverse planes, respectively), age, gender, patient weight, height, and femur morphology as independent parameters and Fracture Risk Index (FRI) as a dependent parameter. Strain-based criteria were used for the calculation of FRI with the maximum principal strain obtained from quantitative computed tomography-based finite element analysis. The statistical result shows that [Formula: see text] [Formula: see text], age [Formula: see text], true moment length [Formula: see text], gender [Formula: see text], FNA [Formula: see text], height [Formula: see text], and FSL [Formula: see text] significantly affect FRI where [Formula: see text] is the most influential parameter. The significance of two-level interaction [Formula: see text] and three-level interaction [Formula: see text] shows that the effect of parameters is dissimilar and depends on other parameters suggesting the variability of FRI from person to person.

The effect of functionally graded materials on bone remodeling around osseointegrated trans-femoral prostheses

Osseointegrated trans-femoral fixations have been used as alternatives for conventional sockets in recent years. Despite numerous advantages, the dissimilarity of the mechanical properties between bone and implant has led to issues in periprosthetic bone adaptation. This study aims to address these issues by proposing fixations made of functionally graded materials (FGMs). The computational study of bone remodeling was performed by linking a bone remodeling algorithm to the finite element analysis. The 3D model of the femur was created by computerized tomography (CT) scan images, and a Titanium fixture, along with nine Titanium/Hydroxyapatite FGM fixtures, were modeled. The analyses revealed evident advantages for the FGM fixtures over the conventionally used Titanium fixtures. Furthermore, it was shown that the gradation direction considerably affects the bone adaptation procedure. The results showed that using a radial FGM with low-stiffness material in the outer layer and less metal composition significantly improves the bone remodeling behavior.

Biomechanical evaluation of axial-loading simulated experiment in wrist fractures: a finite element analysis

Objective

To assess the biomechanical properties that influence wrist fracture, so as to provide the theoretical basis for simulation experiments to aid the optimal design of wrist protectors.

Methods

Six cadaveric wrists were included as experimental specimens. Wrist specimens wearing wrist protectors formed the experimental group and unprotected wrist specimens formed the control group. The wrist specimens were axially loaded under physiological loads and the stress magnitude and distribution of the experimental and control groups were obtained. A three-dimensional wrist finite element model of a healthy volunteer was developed to verify the rationality and effectiveness of the cadaveric wrist models.

Results

Under normal physiological loads, the stress on the radioulnar palmar unit was high and manifested in the form of pressure, while the stress on the radioulnar dorsal unit was lower and manifested in the form of tension. The stresses on the radial distal palmar, ulnar distal palmar, radial distal dorsal, ulnar distal dorsal, radial proximal palmar and ulnar proximal palmar units in the experimental group were less than those in the control group.

Conclusion

Under physiological loads, wearing a wrist protector can reduce the stress on the radioulnar distal palmar, radioulnar proximal palmar and radioulnar distal dorsal units, while having no obvious effect on the radioulnar proximal dorsal units.

A High-Intensity Exercise Intervention Improves Older Women Lumbar Spine and Distal Tibia Bone Microstructure and Function: A 20-Week Randomized Controlled Trial

Introduction: The effects of ageing on bone can be mitigated with different types of physical training, such as power training. However, stimuli that combine increasing external and internal loads concomitantly may improve bone quality. The goal of this study was to assess the efficacy of a combined power and plyometric training on lumbar spine and distal tibia microstructure and function. Methods: 38 sedentary elderly women between 60 and 70 years were randomly allocated in experimental (N = 21) and control group (N = 17). The effects of the 20-week protocol on lumbar spine microstructure and tibia microstructure and function were assessed by trabecular bone score (TBS), high resolution peripheral quantitative computed tomography (HR-pQCT) and microfinite element analysis. Results: when compared to the effects found in the control group, the experimental group showed significant improvements in lumbar spine TBS (Hedges' g = 0.77); and in distal tibia trabecular thickness (g = 0.82) and trabecular bone mineral density (g=0.63). Conclusion: our findings underscore the effectiveness of the proposed intervention, suggesting it as a new strategy to slow down and even reverse the structural and functional losses in the skeletal system due to ageing.

Prediction of lumbar vertebral strength of elderly men based on quantitative computed tomography images using machine learning

The parameters extracted from quantitative computed tomography (QCT) images were used to predict vertebral strength through machine learning models, and the highly accurate prediction indicated that it may be a promising approach to assess fracture risk in clinics.

INTRODUCTION

Vertebral fracture is common in elderly populations. The main factor contributing to vertebral fracture is the reduced vertebral strength. This study aimed to predict vertebral strength based on clinical QCT images by using machine learning.

METHODS

Eighty subjects with QCT data of lumbar spine were randomly selected from the MrOS cohorts. L1 vertebral strengths were computed by QCT-based finite element analysis. A total of 58 features of each L1 vertebral body were extracted from QCT images, including grayscale distribution, grayscale values of 39 partitioned regions, BMD_QCT, structural rigidity, axial rigidity, and BMD_QCTA_min. Feature selection and dimensionality reduction were used to simplify the 58 features. General regression neural network and support vector regression models were developed to predict vertebral strength. Performance of prediction models was quantified by the mean squared error, the coefficient of determination, the mean bias, and the SD of bias.

RESULTS

The 58 parameters were simplified to five features (grayscale value of the 60% percentile, grayscale values of three specific partitioned regions, and BMD_QCTA_min) and nine principal components (PCs). High accuracy was achieved by using the five features or the nine PCs to predict vertebral strength.

CONCLUSIONS

This study provided an effective approach to predict vertebral strength and showed that it may have great potential in clinical applications for noninvasive assessment of vertebral fracture risk.

Assessment of finite element models for prediction of osteoporotic fracture

With increasing life expectancy and mortality rates, the burden of osteoporotic hip fractures is continually on an upward trend. In terms of prevention, there are several osteoporosis treatment strategies such as anti-resorptive drug treatments, which attempt to retard the rate of bone resorption, while promoting the rate of formation. With respect to prediction, several studies have provided insights into obtaining bone strength by non-invasive means through the application of FE analysis. However, what valuable information can we obtain from FE studies that have focused on osteoporosis research, with respect to the prediction of osteoporotic fractures? This paper aims to fine studies that have used FE analysis to predict fractures in the proximal femur through a systematic search of literature using PUBMED, with the main objective of supporting the diagnosis of osteoporosis. The focus of these FE studies is first discussed, and the methodological aspects are summarized, by mainly comparing and contrasting their meshing properties, material properties, and boundary conditions. The implications of these methodological differences in FE modelling processes and propositions with the aim of consolidating or minimalizing these differences are further discussed. We proved that studies need to start converging in terms of their input parameters to make the FE method applicable to clinical settings. This, in turn, will decrease the time needed for in vitro tests. Current advancements in FE analysis need to be consolidated before any further steps can be taken to implement engineering analysis into the clinical scenario.

Whole body vibration with rest days could improve bone quality of distal femoral metaphysis by regulating trabecular arrangement

Low-magnitude, high-frequency vibration (LMHFV) with rest days (particularly seven rest days) was considerably effective in improving the morphological and mechanical properties of rat proximal femur. However, current knowledge is limited regarding the possible benefit of this mechanical regimen to other bone sites and whether the optimal rest days are the same. This study followed our previous experiment on LMHFV loading with rest days for three-month-old male Wistar rats. The experiment involved seven groups, namely, vibrational loading for X day followed with X day rest (X=1, 3, 5, 7), daily vibrational loading, tail suspension and baseline control. Micro-computed tomography (micro-CT) scanning was used to evaluate the microarchitecture of the distal femoral trabecular bone. Micro-CT image-based microfinite element analysis was performed for each distal femoral metaphysis. LMHFV with rest days substantially changed the trabecular arrangement from remarkably plate-like to rod-like. Vibrational loading with 1 day rest was substantially effective in improving the architecture and apparent- and tissuelevel mechanical properties of the rat distal femoral metaphysis. This study may provide an improved understanding of the sitespecific responses of bone tissue to LMHFV with rest days for a substantially effective therapy of a targeted bone site.

Micro-computed tomography and mechanical evaluation of trabecular bone structure in osteopenic and osteoporotic fractures

PURPOSE

Osteopenia and osteoporosis are the two most common musculoskeletal disorders in the elderly population. We determined whether osteopenic and osteoporotic patients with fractures exhibit differences in trabecular morphology and biomechanical properties of bone.

METHODS

Fourteen osteopenic patients and 28 osteoporotic patients with hip fractures who underwent hemiarthroplasty for proximal femoral fractures caused by low-energy injury were included. Bone mineral density (BMD) measurements were performed. Compression tests and high-resolution micro-computed tomography were used to assess cancellous bone samples obtained from the principal compressive region of the femoral head.

RESULTS

The BMD values were lower in the osteoporotic patients than in the osteopenic patients ( p < 0.05). There was a significant difference in the yield stress values between the groups ( p < 0.05). However, no significant differences in the strain energy density, stiffness and Young's modulus were observed between the groups ( p > 0.05). The mean maximum stress was significantly higher in the osteoporotic patients than in the osteopenic patients ( p < 0.05). Although structural parameters, including bone volume (BV), BV fraction, trabecular thickness, trabecular connectivity density and trabecular number, were higher in the osteopenic patients, the differences were not significant ( p > 0.05). Trabecular separation values were significantly higher in the osteoporotic patients ( p < 0.05).

CONCLUSION

Our results showed that the trabecular morphology and biomechanical properties of bone were not significantly different between osteopenic and osteoporotic patients in terms of some parameters.

Quantitative Computed Tomography (QCT) derived Bone Mineral Density (BMD) in finite element studies: a review of the literature

BACKGROUND

Finite element modeling of human bone provides a powerful tool to evaluate a wide variety of outcomes in a highly repeatable and parametric manner. These models are most often derived from computed tomography data, with mechanical properties related to bone mineral density (BMD) from the x-ray energy attenuation provided from this data. To increase accuracy, many researchers report the use of quantitative computed tomography (QCT), in which a calibration phantom is used during image acquisition to improve the estimation of BMD. Since model accuracy is dependent on the methods used in the calculation of BMD and density-mechanical property relationships, it is important to use relationships developed for the same anatomical location and using the same scanner settings, as these may impact model accuracy. The purpose of this literature review is to report the relationships used in the conversion of QCT equivalent density measures to ash, apparent, and/or tissue densities in recent finite element (FE) studies used in common density-modulus relationships. For studies reporting experimental validation, the validation metrics and results are presented.

RESULTS

Of the studies reviewed, 29% reported the use of a dipotassium phosphate (K2HPO4) phantom, 47% a hydroxyapatite (HA) phantom, 13% did not report phantom type, 7% reported use of both K2HPO4 and HA phantoms, and 4% alternate phantom types. Scanner type and/or settings were omitted or partially reported in 31% of studies. The majority of studies used densitometric and/or density-modulus relationships derived from different anatomical locations scanned in different scanners with different scanner settings. The methods used to derive various densitometric relationships are reported and recommendations are provided toward the standardization of reporting metrics.

CONCLUSIONS

This review assessed the current state of QCT-based FE modeling with use of clinical scanners. It was found that previously developed densitometric relationships vary by anatomical location, scanner type and settings. Reporting of all parameters used when referring to previously developed relationships, or in the development of new relationships, may increase the accuracy and repeatability of future FE models.

Assessment of femur geometrical parameters using EOS™ imaging technology in patients with atypical femur fractures; preliminary results

Atypical femur fractures (AFF) arise in the subtrochanteric and diaphyseal regions. Because of this unique distribution, we hypothesized that patients with AFF demonstrate specific geometrical variations of their lower limb whereby baseline tensile forces applied to the lateral cortex are higher and might favor the appearance of these rare stress fractures, when exposed to bisphosphonates. Using the low irradiation 2D-3D X-ray scanner EOS™ imaging technology we aimed to characterize and compare femur geometric parameters between women who sustained bisphosphonate-associated AFF and those who had experienced similar duration of exposure to bisphosphonates but did not sustain fractures. Conditional logistic regression models were constructed to estimate the association between selected geometric parameters and the occurrence of AFF. We identified 16 Caucasian women with AFF and recruited 16 ethnicity-, sex-, age-, height- and cumulative bisphosphonate exposure-matched controls from local osteoporosis clinics. Compared to controls, those with AFF had more lateral femur bowing (-3.2° SD [3.4] versus -0.8° SD [1.9] p=0.02). In regression analysis, lateral femur bowing was associated with the risk of AFF (aOR 1.54; 95% CI 1.04-2.28, p=0.03). Women who sustained a subtrochanteric AFF demonstrated a lesser femoral neck shaft angle (varus geometry) than those with a fracture at a diaphyseal site (121.9 [3.6]° versus 127.6 [7.2]°, p=0.07), whereas femur bowing was more prominent in those with a diaphyseal fracture compared to those with a subtrochanteric fracture (-4.3 [3.2]° versus -0.9 [2.7]°, p=0.07). Our analyses support that subjects with AFF exhibit femoral geometry parameters that result in higher tensile mechanical load on the lateral femur. This may play a critical role in the pathogenesis of AFF and requires further evaluation in a larger size population.

Comparison of proximal femur and vertebral body strength improvements in the FREEDOM trial using an alternative finite element methodology

Denosumab reduced the incidence of new fractures in postmenopausal women with osteoporosis by 68% at the spine and 40% at the hip over 36 months compared with placebo in the FREEDOM study. This efficacy was supported by improvements from baseline in vertebral (18.2%) strength in axial compression and femoral (8.6%) strength in sideways fall configuration at 36 months, estimated in Newtons by an established voxel-based finite element (FE) methodology. Since FE analyses rely on the choice of meshes, material properties, and boundary conditions, the aim of this study was to independently confirm and compare the effects of denosumab on vertebral and femoral strength during the FREEDOM trial using an alternative smooth FE methodology. Unlike the previous FE study, effects on femoral strength in physiological stance configuration were also examined. QCT data for the proximal femur and two lumbar vertebrae were analyzed by smooth FE methodology at baseline, 12, 24, and 36 months for 51 treated (denosumab) and 47 control (placebo) subjects. QCT images were segmented and converted into smooth FE models to compute bone strength. L1 and L2 vertebral bodies were virtually loaded in axial compression and the proximal femora in both fall and stance configurations. Denosumab increased vertebral body strength by 10.8%, 14.0%, and 17.4% from baseline at 12, 24, and 36 months, respectively (p<0.0001). Denosumab also increased femoral strength in the fall configuration by 4.3%, 5.1%, and 7.2% from baseline at 12, 24, and 36 months, respectively (p<0.0001). Similar improvements were observed in the stance configuration with increases of 4.2%, 5.2%, and 5.2% from baseline (p≤0.0007). Differences between the increasing strengths with denosumab and the decreasing strengths with placebo were significant starting at 12 months (vertebral and femoral fall) or 24 months (femoral stance). Using an alternative smooth FE methodology, we confirmed the significant improvements in vertebral body and proximal femur strength previously observed with denosumab. Estimated increases in strength with denosumab and decreases with placebo were highly consistent between both FE techniques.

Assessment of Hip Fracture Risk Using Cross-Section Strain Energy Determined by QCT-Based Finite Element Modeling

Accurate assessment of hip fracture risk is very important to prevent hip fracture and to monitor the effect of a treatment. A subject-specific QCT-based finite element model was constructed to assess hip fracture risk at the critical locations of femur during the single-leg stance and the sideways fall. The aim of this study was to improve the prediction of hip fracture risk by introducing a novel failure criterion to more accurately describe bone failure mechanism. Hip fracture risk index was defined using cross-section strain energy, which is able to integrate information of stresses, strains, and material properties affecting bone failure. It was found that the femoral neck and the intertrochanteric region have higher fracture risk than other parts of the femur, probably owing to the larger content of cancellous bone in these regions. The study results also suggested that women are more prone to hip fracture than men. The findings in this study have a good agreement with those clinical observations reported in the literature. The proposed hip fracture risk index based on strain energy has the potential of more accurate assessment of hip fracture risk. However, experimental validation should be conducted before its clinical applications.

Association of atypical femoral fractures with bisphosphonate use by patients with varus hip geometry

BACKGROUND

There is increasing evidence associating "atypical" femoral fractures with prolonged exposure to bisphosphonate therapy. The cause of these fractures is unknown and likely multifactorial. This study evaluated the hypothesis that patients with primary osteoporosis who sustain atypical femoral fracture(s) while on chronic bisphosphonate therapy have a more varus proximal femoral geometry than patients who use bisphosphonates for primary osteoporosis but do not sustain a femoral fracture.

METHODS

The femoral neck-shaft angle was measured on the radiographs of 111 patients with atypical femoral shaft fracture(s) and thirty-three asymptomatic patients; both groups were on chronic bisphosphonate therapy. Patients with characteristic lateral cortical thickening, stress lines, and thigh pain were included in the fracture group.

RESULTS

The mean neck-shaft angle of the patients who sustained atypical femoral fracture(s) while taking bisphosphonates (case group) differed significantly from that of the patients on bisphosphonate therapy without a fracture (129.5° versus 133.8°; p < 0.001). Fifty-three (48%) of the patients in the case group had a neck-shaft angle that was lower than the lowest angle in the control group (128°). Side-to-side comparison in patients with a unilateral pathologic involvement and an asymptomatic contralateral lower limb did not demonstrate any significant difference between the neck-shaft angles in the two limbs.

CONCLUSIONS

Patients on chronic bisphosphonate therapy who presented with atypical femoral fracture(s) had more varus proximal femoral geometry than those who took bisphosphonates without sustaining a fracture. Although no causative effect can be determined, a finding of varus geometry may help to better identify patients at risk for fracture after long-term bisphosphonate use.

Are we taking full advantage of the growing number of pharmacological treatment options for osteoporosis?

We are becoming increasingly aware that the manner in which our skeleton ages is not uniform within and between populations. Pharmacological treatment options with the potential to combat age-related reductions in skeletal strength continue to become available on the market, notwithstanding our current inability to fully utilize these treatments by accounting for an individual's unique biomechanical needs. Revealing new molecular mechanisms that improve the targeted delivery of pharmaceuticals is important; however, this only addresses one part of the solution for differential age-related bone loss. To improve current treatment regimes, we must also consider specific biomechanical mechanisms that define how these molecular pathways ultimately impact whole bone fracture resistance. By improving our understanding of the relationship between molecular and biomechanical mechanisms, clinicians will be better equipped to take full advantage of the mounting pharmacological treatments available. Ultimately this will enable us to reduce fracture risk among the elderly more strategically, more effectively, and more economically. In this interest, the following review summarizes the biomechanical basis of current treatment strategies while defining how different biomechanical mechanisms lead to reduced fracture resistance. It is hoped that this may serve as a template for the identification of new targets for pharmacological treatments that will enable clinicians to personalize care so that fracture incidence may be globally reduced.

How accurately can we predict the fracture load of the proximal femur using finite element models?

BACKGROUND

Current clinical methods for fracture prediction rely on two-dimensional imaging methods such as dual-energy X-ray absorptiometry and have limited predictive value. Several researchers have tried to integrate three-dimensional imaging techniques with the finite element (FE) method to improve the accuracy of fracture predictions. Before FE models could be used in clinical settings, a thorough validation of their accuracy is required. In this paper, we try to evaluate the current state of accuracy of subject-specific FE models that are used for prediction of the fracture load of proximal femora.

METHODS

All the studies that have used FE for prediction of fracture load and have compared the predicted fracture load with experimentally measured fracture loads in vitro are identified through a systematic search of the literature. A quantitative analysis of the results of those studies has been carried out to determine the absolute prediction error, percentage error, and linear correlations between predicted and measured fracture loads.

FINDINGS

The reported coefficients of determination (R(2)) vary between 0.773 and 0.96 while the percentage error in prediction of fracture load varies between 5 and 46% with most studies reporting percentage errors between 10 and 20%.

INTERPRETATION

We conclude that FE models, which are currently used only experimentally, are in general more accurate than clinically used fracture risk assessment techniques. However, the accuracy of FE models depends on the details of their modeling methodologies. Therefore, modeling procedures need to be optimized and standardized before FE could be used in clinical settings.

Potential of in vivo MRI-based nonlinear finite-element analysis for the assessment of trabecular bone post-yield properties

PURPOSE

Bone strength is the key factor impacting fracture risk. Assessment of bone strength from high-resolution (HR) images have largely relied on linear micro-finite element analysis (μFEA) even though failure always occurs beyond the yield point, which is outside the linear regime. Nonlinear μFEA may therefore be more informative in predicting failure behavior. However, existing nonlinear models applied to trabecular bone (TB) have largely been confined to micro-computed tomography (μCT) and, more recently, HR peripheral quantitative computed tomography (HR-pQCT) images, and typically have ignored evaluation of the post-yield behavior. The primary purpose of this work was threefold: (1) to provide an improved algorithm and program to assess TB yield as well as post-yield properties; (2) to explore the potential benefits of nonlinear μFEA beyond its linear counterpart; and (3) to assess the feasibility and practicality of performing nonlinear analysis on desktop computers on the basis of micro-magnetic resonance (μMR) images obtained in vivo in patients.

METHODS

A method for nonlinear μFE modeling of TB yield as well as post-yield behavior has been designed where material nonlinearity is captured by adjusting the tissue modulus iteratively according to the tissue-level effective strain obtained from linear analysis using a computationally optimized algorithm. The software allows for images at in vivo μMRI resolution as input with retention of grayscale information. Associations between axial stiffness estimated from linear analysis and yield as well as post-yield parameters from nonlinear analysis were investigated from in vivo μMR images of the distal tibia (N = 20; ages: 58-84) and radius (N = 20; ages: 50-75).

RESULTS

All simulations were completed in 1 h or less for 61 strain levels using a desktop computer (dual quad-core Xeon 3.16 GHz CPUs equipped with 40 GB of RAM). Although yield stress and ultimate stress correlated strongly (R(2) > 0.95, p < 0.001) with axial stiffness, toughness correlated moderately at the distal tibia (R(2) = 0.81, p < 0.001) and only weakly at the distal radius (R(2) = 0.34, p = 0.007). Further, toughness was found to vary by up to 16% for bone of very similar axial stiffness (<2%).

CONCLUSIONS

The work demonstrates the practicality of nonlinear μFE simulations at in vivo μMRI resolution, as well as its potential for providing additional information beyond that obtainable from linear analysis. The data suggest that a direct assessment of toughness may provide information not captured by stiffness.

Relationship between microstructure, material distribution, and mechanical properties of sheep tibia during fracture healing process

The aim of this study was to investigate the relationship between microstructural parameters, material distribution, and mechanical properties of sheep tibia at the apparent and tissue levels during the fracture healing process. Eighteen sheep underwent tibial osteotomy and were sacrificed at 4, 8, and 12 weeks. Radiographs and micro-computed tomography (micro-CT) scanning were taken for microstructural assessment, material distribution evaluation, and micro-finite element analysis. A displacement of 5% compressive strain on the longitudinal direction was applied to the micro-finite element model, and apparent and tissue-level mechanical properties were calculated. Principle component analysis and linear regression were used to establish the relationship between principle components (PCs) and mechanical parameters. Visible bony callus formation was observed throughout the healing process from radiographic assessment. Apparent mechanical property increased at 8 weeks, but tissue-level mechanical property did not increase significantly until 12 weeks. Three PCs were extracted from microstructural parameters and material distribution, which accounted for 87.592% of the total variation. The regression results showed a significant relationship between PCs and mechanical parameters (R>0.8, P<0.05). Results of this study show that microstructure and material distribution based on micro-CT imaging could efficiently predict bone strength and reflect the bone remodeling process during fracture healing, which provides a basis for exploring the fracture healing mechanism and may be used as an approach for fractured bone strength assessment.