Prediction of fracture location in the proximal femur using finite element models

Magnetic resonance imaging based finite element modelling of the proximal femur: a short-term in vivo precision study

Proximal femoral fractures are a serious life-threatening injury with high morbidity and mortality. Magnetic resonance (MR) imaging has potential to non-invasively assess proximal femoral bone strength in vivo through usage of finite element (FE) modelling (a technique referred to as MR-FE). To precisely assess bone strength, knowledge of measurement error associated with different MR-FE outcomes is needed. The objective of this study was to characterize the short-term in vivo precision errors of MR-FE outcomes (e.g., stress, strain, failure loads) of the proximal femur for fall and stance loading configurations using 13 participants (5 males and 8 females; median age: 27 years, range: 21-68), each scanned 3 times. MR-FE models were generated, and mean von Mises stress and strain as well as principal stress and strain were calculated for 3 regions of interest. Similarly, we calculated the failure loads to cause 5% of contiguous elements to fail according to the von Mises yield, Brittle Coulomb-Mohr, normal principal, and Hoffman stress and strain criteria. Precision (root-mean squared coefficient of variation) of the MR-FE outcomes ranged from 3.3% to 11.8% for stress and strain-based mechanical outcomes, and 5.8% to 9.0% for failure loads. These results provide evidence that MR-FE outcomes are a promising non-invasive technique for monitoring femoral strength in vivo.

Scaphoid numerical simulation of the critical loading until fracture

The numerical study of the scaphoid fracture, although it is relatively unexplored, can be of great clinical interest since it is highly common and can result in temporary or persistent disability. In this manuscript, seven combinations of boundary conditions and contacts between adjacent bones, together with four different loads, simulating real hand movements, are assessed. Three different fracture criteria for bones are employed to study the failure of the scaphoid with the aforementioned combination of interaction conditions. The results offer an interesting view of the accuracy of the possible interaction between adjacent bones. For future calculation, it would be possible to choose a combination of the balance between precision and computational cost savings. This study provides a comprehensive assessment into the modeling of the scaphoid bone and its interactions with adjacent bones. The findings reveal that various choices of interactions can yield similar results, allowing for flexibility in selecting interaction models based on desired accuracy or computational efficiency. Ultimately, this study establishes a foundational understanding for future research on modeling scaphoid motion.

Application of Optical Photothermal Infrared (O-PTIR) Spectroscopy for Assessment of Bone Composition at the Submicron Scale

The molecular basis of bone structure and strength is mineralized collagen fibrils at the submicron scale (∼500 nm). Recent advances in optical photothermal infrared (O-PTIR) spectroscopy allow the investigation of bone composition with unprecedented submicron spatial resolution, which may provide new insights into factors contributing to underlying bone function. Here, we investigated (i) whether O-PTIR-derived spectral parameters correlated to standard attenuated total reflection (ATR) Fourier transform infrared spectroscopy spectral data and (ii) whether O-PTIR-derived spectral parameters, including heterogeneity of tissue, contribute to the prediction of proximal femoral bone stiffness. Analysis of serially demineralized bone powders showed a significant correlation (r = 0.96) between mineral content quantified using ATR and O-PTIR spectroscopy, indicating the validity of this technique in assessing bone mineralization. Using femoral neck sections, the principal component analysis showed that differences between O-PTIR and ATR spectra were primarily attributable to the phosphate ion (PO4) absorbance band, which was typically shifter toward higher wavenumbers in O-PTIR spectra. Additionally, significant correlations were found between hydrogen phosphate (HPO4) content (r = 0.75) and carbonate (CO3) content (r = 0.66) quantified using ATR and O-PTIR spectroscopy, strengthening the validity of this method to assess bone mineral composition. O-PTIR imaging of individual trabeculae at 500 nm pixel resolution illustrated differences in submicron composition in the femoral neck from bones with different stiffness. O-PTIR analysis showed a significant negative correlation (r = -0.71) between bone stiffness and mineral maturity, reflective of newly formed bone being an important contributor to bone function. Finally, partial least squares regression analysis showed that combining multiple O-PTIR parameters (HPO4 content and heterogeneity, collagen integrity, and CO3 content) could significantly predict proximal femoral stiffness (R2 = 0.74, error = 9.7%) more accurately than using ATR parameters. Additionally, we describe new findings in the effects of bone tissue orientation in the O-PTIR spectra. Overall, this study highlights a new application of O-PTIR spectroscopy that may provide new insights into molecular-level factors underlying bone mechanical competence.

The Influence of Static Load and Sideways Impact Fall on Extramedullary Bone Plates Used to Treat Intertrochanteric Femoral Fracture: A Preclinical Strength Assessment

Hip fracture accounts for a large number of hospitalizations, thereby causing substantial economic burden. Majority (> 90%) of all hip fractures are associated to sideways fall. Studies on sideways fall usually involve loading at quasi-static or at constant displacement rate, which neglects the physics of actual fall. Understanding femur resonance frequency and associated mode shapes excited by dynamic loads is also critical. Two commercial extramedullary implants, proximal femoral locking plate (PFLP) and variable angle dynamic hip screw (VA-DHS), were chosen to carry out the preclinical assessments on a simulated Evans-I type intertrochanteric fracture. In this study, we hypothesized that the behavior of the implant depends on the loading types-axial static and transverse impact-and a rigid implanted construct will absorb less impact energy for sideways fall. The in silico models were validated using experimental measurements of full-field strain data obtained from a 2D digital image correlation (DIC) study. Under peak axial load of 3 kN, PFLP construct predicted greater axial stiffness (1.07 kN/mm) as opposed to VA-DHS (0.85 kN/mm), although the former predicted slightly higher proximal stress shielding. Further, with greater mode 2 frequency, PFLP predicted improved performance in resisting bending due to sideways fall as compared to the other implant. Overall, the PFLP implanted femur predicted the least propensity to adverse stress intensities, suggesting better structural rigidity and higher capacity in protecting the fractured femur against fall.

BHLHE40 promotes osteoclastogenesis and abnormal bone resorption via c-Fos/NFATc1

Background

Dysregulated osteoclast activity due to altered osteoclast differentiation causes multiple bone diseases. Osteoclasts are multinucleated giant cells derived from hematopoietic stem cells and play a major role in bone absorption. However, the mechanisms underlying the tight regulation of osteoclast differentiation in multiple pathophysiological status remain unknown.

Results

We showed that Bhlhe40 upregulation is tightly associated with osteoclast differentiation and osteoporosis. Functionally, Bhlhe40 promoted osteoclast differentiation in vitro, and Bhlhe40 deficiency led to increased bone mass and decreased osteoclast differentiation in vivo. Moreover, Bhlhe40 deficient mice resisted estrogen deficiency and aging-induced osteoporosis. Mechanism study showed that the increase in bone mass due to Bhlhe40 deficiency was a cell intrinsic defect in osteoclast differentiation in these mice. BHLHE40 upregulated the gene expression of Fos and Nfatc1 by directly binding to their promoter regions. Notably, inhibition of Fos/Nfatc1 abrogated the enhanced osteoclast differentiation induced by BHLHE40 overexpression.

Conclusions

Our research reveals a novel Bhlhe40/c-Fos/Nfatc1 axis involved in regulating osteoclastogenesis and shows that osteoporosis caused by estrogen deficiency and aging can be rescued by regulating Bhlhe40 in mice. This may help in the development of a new strategy for the treatment of osteoporosis.

Comparison of Cortical Bone Fracture Patterns Under Compression Loading Using Finite Element–Discrete Element Numerical Modeling Approach and Destructive Testing

Finite element analysis may not be the only method by which bone fracture initiation and propagation may be analyzed. This study compares fracture patterns generated from compression testing of bone to fracture patterns generated using a combination of both the finite element method (FEM) and discrete element method (DEM) as defined by the finite discrete element method (FDEM). Before testing, a three-dimensional bone model was developed using CT. Force and displacement data were collected during testing. The tested specimen was reimaged using CT. The solid model was discretized and material properties adjusted such that finite element-discrete element macro behavior matched the force-displacement data. A qualitative comparison of the fracture patterns demonstrates that FDEM can successfully be used to simulate and predict fracturing in bone, with this study representing the first time this has been done and reported.

Determinants of fracture type in the proximal femur: Biomechanical study of fresh frozen cadavers and finite element models

BACKGROUND

Proximal femur fractures are usually categorized as either a cervical or trochanteric fracture, but the relationship between fracture type and fall direction is not clear. By cadaveric mechanical testing and finite element analysis (FEA), the aims of this research were to verify the factors that define the proximal femur fracture type and to clarify the change in stress distribution based on fall direction.

METHODS

From fresh frozen cadavers, we obtained 26 proximal femora including ten pairs of 20 femora. We conducted quasi-static compression tests in two fall patterns (lateral and posterolateral), and identified the fracture type. We then examined the relationship between fracture type and the following explanatory variables: age, sex, neck shaft angle, femoral neck length, bone mineral density (cervical and trochanteric), and fall direction. In addition, for the ten pairs of femurs, the effect of fall direction on fracture type was examined by comparing the left and right sides. In addition, we generated the proximal femur finite element (FE) models from computed tomography data to simulate and verify the change of external force in different fall directions.

RESULTS

In mechanical tests, only fall direction was found to have a significant relationship with fracture type (p = 0.0227). The posterolateral fall group had a significantly higher incidence of trochanteric fractures than lateral fall group (p = 0.0325). According to FEA, the equivalent stress in the lateral fall was found to be more concentrated in the cervical area than in the posterolateral fall.

CONCLUSION

In proximal femur fractures, fall direction was significantly associated with fracture type; in particular, trochanteric fractures were more likely to occur following a posterolateral fall than a lateral fall.

Finite element analysis of trabecular bone microstructure using CT imaging and continuum mechanical modelling

Purpose

Osteoporosis is a bone disease associated with enhanced bone loss, microstructural degeneration, and fracture‐risk. Finite element (FE) modeling is used to estimate trabecular bone (Tb) modulus from high‐resolution three‐dimensional (3‐D) imaging modalities including micro‐computed tomography (CT), magnetic resonance imaging (MRI), and high‐resolution peripheral quantitative CT (HR‐pQCT). This paper validates an application of voxel‐based continuum finite element analysis (FEA) to predict Tb modulus from clinical CT imaging under a condition similar to in vivo imaging by comparing with measures derived by micro‐CT and experimental approaches.

Method

Voxel‐based continuum FEA methods for CT imaging were implemented using linear and nonlinear models and applied on distal tibial scans under a condition similar to in vivo imaging. First, tibial axis in a CT scan was aligned with the coordinate z‐axis at 150 μm isotropic voxels. FEA was applied on an upright cylindrical volume of interests (VOI) with its axis coinciding with the tibial bone axis. Voxel volume, edge, and vertex elements and their connectivity were defined as per the isotropic image grid. A calibration phantom was used to calibrate CT numbers in Hounsfield unit to bone mineral density (BMD) values, which was then converted into calcium hydroxyapatite (CHA) density. Mechanical properties at each voxel volume element was defined using its ash‐density defined on CT‐derived CHA density. For FEA, the bottom surface of the cylindrical VOI was fixed and a constant displacement was applied along the z‐direction at each vertex element on the top surface to simulate a physical axial compressive loading condition. Finally, a Poisson's ratio of 0.3 was applied, and Tb modulus (MPa) was computed as the ratio of average von Mises stress (MPa) of volume elements on the top surface and the applied displacement. FEA parameters including mesh element size, substep number, and different tolerance values were optimized.

Results

CT‐derived Tb modulus values using continuum FEA showed high linear correlation with the micro‐CT‐derived reference values (r ∈ [0.87 0.90]) as well as experimentally measured values (r ∈ [0.80 0.87]). Linear correlation of computed modulus with their reference values using continuum FEA with linear modeling was comparable with that obtained by nonlinear modeling. Nonlinear continuum FEA‐based modulus values (mean of 1087.2 MPa) showed greater difference from their reference values (mean of 1498.9 MPa using micro‐CT‐based FEA) as compared with linear continuum methods. High repeat CT scan reproducibility (intra‐class correlation [ICC] = 0.98) was observed for computed modulus values using both linear and nonlinear continuum FEA. It was observed that high stress regions coincide with Tb microstructure as fuzzily characterized by BMD values. Distributions of von Mises stress over Tb microstructure and marrow regions were significantly different (p < 10^–8).

Conclusion

Voxel‐based continuum FEA offers surrogate measures of Tb modulus from CT imaging under a condition similar to in vivo imaging that alleviates the need for segmentation of Tb and marrow regions, while accounting for bone distribution at the microstructural level. This relaxation of binary segmentation will extend the scope of FEA application to assess mechanical properties of bone microstructure at relatively low‐resolution imaging.

Hinge fractures reaching the tibial plateau can be caused by forcible opening of insufficient posterior osteotomy during open-wedge high tibial osteotomy

PURPOSE

The purpose of this study was to use the finite element method (FEM) to reproduce fracture lines that reach the lateral tibial plateau during open-wedge high tibial osteotomy (OWHTO) in patients with Type III lateral hinge fracture (LHF). It was hypothesized that the FEM could clarify biomechanical causes of Type III LHF, enabling prevention of adverse complications.

METHODS

This study used the nonlinear FEM to analyze the data of eight knees in eight patients (two males and six females) with Type III LHF among 82 patients who underwent OWHTO, as well as the data of eight individuals with no LHF. To predict the onset of Type III LHF, simulation models were also developed in which posterior osteotomy sufficiency varied from 50% to perfect, the latter defined as osteotomy reaching the hinge point.

RESULTS

Real-life instances of Type III LHF caused by insufficient posterior osteotomy were reproduced in all patient-specific FEM models, and these models accurately predicted fracture types and locations. During opening of the osteotomy gap, the fracture line reached the lateral tibial plateau, and extended vertically from the end of the insufficient posterior osteotomy, avoiding the rigid proximal tibiofibular joint. In contrast, sufficient posterior osteotomy resulted in a lack of LHF. Posterior osteotomy extension ≥ 70% of the width of the osteotomy plane was the cut-off value to prevent Type III LHF.

CONCLUSION

Forced opening of insufficient posterior osteotomy was found to be a biomechanical cause of Type III LHF that extended perpendicularly to the lateral tibial plateau, avoiding the proximal tibiofibular joint. The clinical significance of this study is that sufficient posterior osteotomy during OWHTO, defined as at least 70% of the width of the osteotomy plane, can prevent Type III LHF.

Finite element modeling of multiple density materials of bone specimens for biomechanical behavior evaluation

The finite elements method allied with the computerized axial tomography (CT) is a mathematical modeling technique that allows constructing computational models for bone specimens from CT data. The objective of this work was to compare the experimental biomechanical behavior by three-point bending tests of porcine femur specimens with different types of computational models generated through the finite elements’ method and a multiple density materials assignation scheme. Using five femur specimens, 25 scenarios were created with differing quantities of materials. This latter was applied to computational models and in bone specimens subjected to failure. Among the three main highlights found, first, the results evidenced high precision in predicting experimental reaction force versus displacement in the models with larger number of assigned materials, with maximal results being an R2 of 0.99 and a minimum root-mean-square error of 3.29%. Secondly, measured and computed elastic stiffness values follow same trend with regard to specimen mass, and the latter underestimates stiffness values a 6% in average. Third and final highlight, this model can precisely and non-invasively assess bone tissue mechanical resistance based on subject-specific CT data, particularly if specimen deformation values at fracture are considered as part of the assessment procedure.

Prediction of Risk Factors for Pathological Fracture After Bone Tumor Biopsy Using Finite Element Analysis

Purpose

We aimed to determine if finite element analysis (FEA) provides useful thresholds for bone biopsy practice patterns.

Methods

The femoral head compression test was performed on rabbit femurs, using FEA to identify the part of the bone that preferentially fractures (n=15/group). Four types of rectangular biopsy holes were made using finite element (FE) models. These models were divided into control (no defect), defect 1 (10% width), defect 2 (20% width), defect 3 (30% width), and defect 4 (40% width) groups (n=15 each). Three types of rectangular biopsy holes (defect A, 27% length; defect B, 40% length; defect C, 53% length) were also made using FE models (n=15 each). The load to failure was then predicted using FEA.

Results

Almost all femurs with no defect were fractured at the femoral shaft in both the femoral head compression test and FEA. The experimental load to failure in intact femurs was predicted well by the FE models (R²=0.74, p<0.001). There was also a strong linear correlation of stiffness between compression test in femurs with no defect and the FEA (R²=0.68, p<0.001). Therefore, the femoral shaft was targeted for FEA. The median predicted loads by FEA were significantly higher for defect 1 than for the other types when testing the widths of the rectangular defects, but there were no significant differences among the three types when testing for defect length.

Conclusion

The FEA results correlated well with those of the femoral head compression test. A width <10% of the circumference length in bone biopsy holes helps minimize bone strength reduction using FEA. It may be useful for orthopedic doctors to perform FEA to avoid pathological fractures after bone tumor biopsy.

Evaluation of fixation after plating of distal radius fractures - a validation study

We used five fresh-frozen cadavers with 10 upper limbs to evaluate by finite element analysis (FEA) the plate fixation of distal radius fractures. The distal radius of the cadavers was fractured using a comminution fracture model. Plate fixation was performed using Synthes VATCP. Compression tests were performed on these specimens and force displacement curves were obtained. FEA was performed using Mechanical Finder. The Keyak, Keller vertebra, Carter, and Matsuyama conversion equations without contact analysis, and the Matsuyama equations with contact analysis, were used for the boundary conditions. We found strong positive correlations with the Matsuyama conversion equations either with or without contact analysis. The validated FEA model will be used for preoperative simulation of actual fractures.

MRI-derived porosity index is associated with whole-bone stiffness and mineral density in human cadaveric femora

Ultrashort echo time (UTE) magnetic resonance imaging (MRI) measures proton signals in cortical bone from two distinct water pools, bound water, or water that is tightly bound to bone matrix, and pore water, or water that is freely moving in the pore spaces in bone. By isolating the signal contribution from the pore water pool, UTE biomarkers can directly quantify cortical bone porosity in vivo. The Porosity Index (PI) is one non-invasive, clinically viable UTE-derived technique that has shown strong associations in the tibia with μCT porosity and other UTE measures of bone water. However, the efficacy of the PI biomarker has never been examined in the proximal femur, which is the site of the most catastrophic osteoporotic fractures. Additionally, the loads experienced during a sideways fall are complex and the femoral neck is difficult to image with UTE, so the usefulness of the PI in the femur was unknown. Therefore, the aim of this study was to examine the relationships between the PI measure in the proximal cortical shaft of human cadaveric femora specimens compared to (1) QCT-derived bone mineral density (BMD) and (2) whole bone stiffness obtained from mechanical testing mimicking a sideways fall. Fifteen fresh, frozen whole cadaveric femora specimens (age 72.1 ± 15.0 years old, 10 male, 5 female) were scanned on a clinical 3-T MRI using a dual-echo UTE sequence. Specimens were then scanned on a clinical CT scanner to measure volumetric BMD (vBMD) and then non-destructively mechanically tested in a sideways fall configuration. The PI in the cortical shaft demonstrated strong correlations with bone stiffness (r = -0.82, P = 0.0014), CT-derived vBMD (r = -0.64, P = 0.0149), and with average cortical thickness (r = -0.60, P = 0.0180). Furthermore, a hierarchical regression showed that PI was a strong predictor of bone stiffness which was independent of the other parameters. The findings from this study validate the MRI-derived porosity index as a useful measure of whole-bone mechanical integrity and stiffness.

Numerical analysis of hip fracture due to a sideways fall

The primary purpose of this paper is to outline a methodology for evaluating the likelihood of cortical bone fracture in the proximal femur in the event of a sideways fall. The approach includes conducting finite element (FE) analysis in which the cortical bone is treated as an anisotropic material, and the admissibility of the stress field is validated both in tension and compression regime. In assessing the onset of fracture, two methodologies are used, namely the Critical Plane approach and the Microstructure Tensor approach. The former is employed in the tension regime, while the latter governs the conditions at failure in compression. The propagation of localized damage is modeled using a constitutive law with embedded discontinuity (CLED). In this approach, the localized deformation is described by a homogenization procedure in which the average properties of cortical tissue intercepted by a macrocrack are established. The key material properties governing the conditions at failure are specified from a series of independent material tests conducted on cortical bone samples tested at different orientations relative to the loading direction. The numerical analysis deals with simulations of experiments involving the sideways fall, and the results are compared with the experimental data. This includes both the evolution of fracture pattern and the local load-displacement characteristics. The proposed approach is numerically efficient, and the results do not display a pathological mesh-dependency. Also, in contrast to the XFEM approach, the analysis does not require any extra degrees of freedom.

Estimation of mass apparent density and Young's modulus of femoral neck-head region

The purpose of this study is to estimate mass apparent density and Young's modulus to investigate biomechanical properties of the proximal femur bone. In this study eleven specimens of sheep femur bone having age between 1-1.25 years and human femur bone having age between 14 and 81years are used. In the present study, the first technique attempts to estimate the density from the image-based Hounsfield unit which is obtained directly from a computed tomography image. The modulus of elasticity is estimated from density-elasticity relation which is available in the literature. Another technique is used to develop a correlation between computed apparent density and greyscale based coefficient obtained by material mapping method using commercial Simpleware ScanIP software. Estimated mean deviation in apparent mass density and Young's modulus is 4.34% and 4.69% in sheep bone and 4.35% and 4.94% in human bone respectively. It is found that apparent density and Young's modulus obtained shows close agreement with values reported in the literature. Moreover, the study attempts to build up a new material model between human and sheep for orthopaedics clinical trials and research in Indian context. In addition, it is also observed that bone mass density of sheep is 1.60 times human. This method can also be useful to study and analyse biomechanical properties of the human femur bone.

Comparison of femur strain under different loading scenarios: Experimental testing

Bone fracture, formation and adaptation are related to mechanical strains in bone. Assessing bone stiffness and strain distribution under different loading conditions may help predict diseases and improve surgical results by determining the best conditions for long-term functioning of bone-implant systems. In this study, an experimentally wide range of loading conditions (56) was used to cover the directional range spanned by the hip joint force. Loads for different stance configurations were applied to composite femurs and assessed in a material testing machine. The experimental analysis provides a better understanding of the influence of the bone inclination angle in the frontal and sagittal planes on strain distribution and stiffness. The results show that the surface strain magnitude and stiffness vary significantly under different loading conditions. For the axial compression, maximal bending is observed at the mid-shaft, and bone stiffness is also maximal. The increased inclination leads to decreased stiffness and increased magnitude of maximum strain at the distal end of the femur. For comparative analysis of results, a three-dimensional, finite element model of the femur was used. To validate the finite element model, strain gauges and digital image correlation system were employed. During validation of the model, regression analysis indicated robust agreement between the measured and predicted strains, with high correlation coefficient and low root-mean-square error of the estimate. The results of stiffnesses obtained from multi-loading conditions experiments were qualitatively compared with results obtained from a finite element analysis of the validated model of femur with the same multi-loading conditions. When the obtained numerical results are qualitatively compared with experimental ones, similarities can be noted. The developed finite element model of femur may be used as a promising tool to estimate proximal femur strength and identify the best conditions for long-term functioning of the bone-implant system in future study.

Hip load capacity and yield load in men and women of all ages

Quantitative computed tomography (QCT) based finite element (FE) models can compute subject-specific proximal femoral strengths, or fracture loads, that are associated with hip fracture risk. These fracture loads are more strongly associated with measured fracture loads than are DXA and QCT measures and are predictive of hip fracture independently of DXA bone mineral density (BMD). However, interpreting FE-computed fracture loads of younger subjects for the purpose of evaluating hip fracture risk in old age is challenging due to limited reference data. The goal of this study was to address this issue by providing reference data for male and female adult subjects of all ages. QCT-based FE models of the left proximal femur of 216 women and 181 men, age 27 to 90 years, from a cohort of Rochester, MN residents were used to compute proximal femoral load capacities, i.e. the maximum loads that can be supported, in single-limb stance and posterolateral fall loading (Stance_LC and Fall_LC, respectively) [US Patent No. 9,245,069] and yield load under fall loading (Fall_yield). To relate these measures to information about hip fracture, the CT scanner and calibration phantom were cross-calibrated with those from our previous prospective study of hip fracture in older fracture and control subjects, the Age Gene/Environment Susceptibility (AGES) Reykjavik cohort. We then plotted Stance_LC_, Fall_LC and Fall_yield versus age for the two cohorts on the same graphs. Thus, proximal femoral strengths in individuals above 70 years of age can be assessed through direct comparison with the FE data from the AGES cohort which were analyzed using identical methods. To evaluate younger individuals, reductions in Stance_LC, Fall_LC and Fall_yield from the time of evaluation to age 70 years can be cautiously estimated from the average yearly cross-sectional decreases found in this study (108 N, 19.4 N and 14.4 N, respectively, in men and 120 N, 19.4 N and 21.6 N, respectively, in women), and the projected fracture loads can be compared with data from the AGES cohort. Although we did not set specific thresholds for identifying individuals at risk of hip fracture, these data provide some guidance and may be used to help establish diagnostic criteria in future. Additionally, given that these data were nearly entirely from Caucasian subjects, future research involving subjects of other races/ethnicities is necessary.

Simulated lesions representative of metastatic disease predict proximal femur failure strength more accurately than idealized lesions

Metastatic disease in bone is characterized by highly amorphous and variable lesion geometry, with increased fracture risk. Assumptions of idealized lesion geometry made in previous finite element (FE) studies of metastatic disease in the proximal femur may not sufficiently capture effects of local stress/strain concentrations on predicted failure strength. The goal of this study was to develop and validate a FE failure model of the proximal femur incorporating artificial defects representative of physiologic metastatic disease. Data from 11 cadaveric femur specimens were randomly divided into either a training set (n = 5) or a test set (n = 6). Clinically representative artificial defects were created, and the femurs were loaded to failure under offset torsion. Voxel-based FE models replicating the experimental setup were created from the training set pre-fracture computed tomography data. Failure loads from the linear model with maximum principal strain failure criterion correlated best with the experimental data (R² = 0.86, p = 0.024). The developed model was found to be reliable when applied to the test dataset with a relatively low RMSE of 46.9 N, mean absolute percent error of 12.7 ± 17.1%, and cross-validation R² = 0.88 (p < 0.001). Models simulating realistic lesion geometry explained an additional 26% of the variance in experimental failure load compared to idealized lesion models (R² = 0.62, p = 0.062). Our validated automated FE model representative of physiologic metastatic disease may improve clinical fracture risk prediction and facilitate research studies of fracture risk during functional activities and with treatment interventions.

MRI-based assessment of proximal femur strength compared to mechanical testing

Half of the women who sustain a hip fracture would not qualify for osteoporosis treatment based on current DXA-estimated bone mineral density criteria. Therefore, a better approach is needed to determine if an individual is at risk of hip fracture from a fall. The objective of this study was to determine the association between radiation-free MRI-derived bone strength and strain simulations compared to results from direct mechanical testing of cadaveric femora. Imaging was conducted on a 3-Tesla MRI scanner using two sequences: one balanced steady-state free precession sequence with 300 μm isotropic voxel size and one spoiled gradient echo with anisotropic voxel size of 234 × 234 × 1500 μm. Femora were dissected free of soft-tissue and 4350-ohm strain-gauges were securely applied to surfaces at the femoral shaft, inferior neck, greater trochanter, and superior neck. Cadavers were mechanically tested with a hydraulic universal test frame to simulate loading in a sideways fall orientation. Sideways fall forces were simulated on MRI-based finite element meshes and bone stiffness, failure force, and force for plastic deformation were computed. Simulated bone strength metrics from the 300 μm isotropic sequence showed strong agreement with experimentally obtained values of bone strength, with stiffness (r = 0.88, p = 0.0002), plastic deformation point (r = 0.89, p < 0.0001), and failure force (r = 0.92, p < 0.0001). The anisotropic sequence showed similar trends for stiffness, plastic deformation point, and failure force (r = 0.68, 0.70, 0.84; p = 0.02, 0.01, 0.0006, respectively). Surface strain-gauge measurements showed moderate to strong agreement with simulated magnitude strain values at the greater trochanter, superior neck, and inferior neck (r = -0.97, -0.86, 0.80; p ≤0.0001, 0.003, 0.03, respectively). The findings from this study support the use of MRI-based FE analysis of the hip to reliably predict the mechanical competence of the human femur in clinical settings.

Heterogeneous Spatial and Strength Adaptation of the Proximal Femur to Physical Activity: A Within-Subject Controlled Cross-Sectional Study

Physical activity (PA) enhances proximal femur bone mass, as assessed using projectional imaging techniques. However, these techniques average data over large volumes, obscuring spatially heterogeneous adaptations. The current study used quantitative computed tomography, statistical parameter mapping, and subject-specific finite element (FE) modeling to explore spatial adaptation of the proximal femur to PA. In particular, we were interested in adaptation occurring at the superior femoral neck and improving strength under loading from a fall onto the greater trochanter. High/long jump athletes (n = 16) and baseball pitchers (n = 16) were utilized as within-subject controlled models as they preferentially load their take-off leg and leg contralateral to their throwing arm, respectively. Controls (n = 15) were included but did not show any dominant-to-nondominant (D-to-ND) leg differences. Jumping athletes showed some D-to-ND leg differences but less than pitchers. Pitchers had 5.8% (95% confidence interval [CI] 3.9%-7.6%) D-to-ND leg differences in total hip volumetric bone mineral density (vBMD), with increased vBMD in the cortical compartment of the femoral neck and trochanteric cortical and trabecular compartments. Voxel-based morphometry analyses and cortical bone mapping showed pitchers had D-to-ND leg differences within the regions of the primary compressive trabeculae, inferior femoral neck, and greater trochanter but not the superior femoral neck. FE modeling revealed pitchers had 4.1% (95% CI 1.4%-6.7%) D-to-ND leg differences in ultimate strength under single-leg stance loading but no differences in ultimate strength to a fall onto the greater trochanter. These data indicate the asymmetrical loading associated with baseball pitching induces proximal femur adaptation in regions associated with weight bearing and muscle contractile forces and increases strength under single-leg stance loading. However, there were no benefits evident at the superior femoral neck and no measurable improvement in ultimate strength to common injurious loading during aging (ie, fall onto the greater trochanter), raising questions as to how to better target these variables with PA. © 2019 American Society for Bone and Mineral Research.

Biomechanics of callus in the bone healing process, determined by specimen-specific finite element analysis

As fractures heal, immature callus formed in the hematoma is calcified by osteoblasts and altered to mature bone. Although the bone strength in the fracture-healing process cannot be objectively measured in clinical settings, bone strength can be predicted by specimen-specific finite element modeling (FEM) of quantitative computed tomography (qCT) scans. FEM predictions of callus strength would enable an objective treatment plan. The present study establishes an equation that converts material properties to bone density and proposes a specimen-specific FEM. In 10 male New Zealand white rabbits, a 10-mm long bone defect was created in the center of the femur and fixed by an external fixator. The callus formed in the defect was extracted after 3-6 weeks, and formed into a (5 × 5 × 5 mm³) cube. The bone density measured by qCT was related to the Young's modulus and the yield stress measured with a mechanical tester. For validation, a 10-mm long bone defect was created in the central femurs of another six New Zealand white rabbits, and fixed by an external fixator. At 3, 4, and 5 weeks, the femur was removed and subjected to Computed tomography (CT) scanning and mechanical testing. A specimen-specific finite element model was created from the CT data. Finally, the bone strength was measured and compared with the experimental value. The bone mineral density σ was significantly and nonlinearly correlated with both the Young's modulus E and the yield stress σ. The material-property conversion equations were E = 0.2391e^8.00ρ and ρ = 30.49σ^2.41. Moreover, the experimental bone strength was significantly linearly correlated with the prospective FEM. We demonstrated the Young's moduli and yield stresses for different bone densities, enabling a FEM of the bone-healing process. An FEM based on these material properties is expected to yield objective clinical judgment criteria.

Perspectives on the non-invasive evaluation of femoral strength in the assessment of hip fracture risk

We reviewed the experimental and clinical evidence that hip bone strength estimated by BMD and/or finite element analysis (FEA) reflects the actual strength of the proximal femur and is associated with hip fracture risk and its changes upon treatment.

INTRODUCTION

The risk of hip fractures increases exponentially with age due to a progressive loss of bone mass, deterioration of bone structure, and increased incidence of falls. Areal bone mineral density (aBMD), measured by dual-energy X-ray absorptiometry (DXA), is the most used surrogate marker of bone strength. However, age-related declines in bone strength exceed those of aBMD, and the majority of fractures occur in those who are not identified as osteoporotic by BMD testing. With hip fracture incidence increasing worldwide, the development of accurate methods to estimate bone strength in vivo would be very useful to predict the risk of hip fracture and to monitor the effects of osteoporosis therapies.

METHODS

We reviewed experimental and clinical evidence regarding the association between aBMD and/orCT-finite element analysis (FEA) estimated femoral strength and hip fracture risk as well as their changes with treatment.

RESULTS

Femoral aBMD and bone strength estimates by CT-FEA explain a large proportion of femoral strength ex vivo and predict hip fracture risk in vivo. Changes in femoral aBMD are strongly associated with anti-fracture efficacy of osteoporosis treatments, though comparable data for FEA are currently not available.

CONCLUSIONS

Hip aBMD and estimated femoral strength are good predictors of fracture risk and could potentially be used as surrogate endpoints for fracture in clinical trials. Further improvements of FEA may be achieved by incorporating trabecular orientations, enhanced cortical modeling, effects of aging on bone tissue ductility, and multiple sideway fall loading conditions.

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.

Numerical modelling of hip fracture patterns in human femur

BACKGROUND AND OBJECTIVE

Hip fracture morphology is an important factor determining the ulterior surgical repair and treatment, because of the dependence of the treatment on fracture morphology. Although numerical modelling can be a valuable tool for fracture prediction, the simulation of femur fracture is not simple due to the complexity of bone architecture and the numerical techniques required for simulation of crack propagation. Numerical models assuming homogeneous fracture mechanical properties commonly fail in the prediction of fracture patterns. This paper focuses on the prediction of femur fracture based on the development of a finite element model able to simulate the generation of long crack paths.

METHODS

The finite element model developed in this work demonstrates the capability of predicting fracture patterns under stance loading configuration, allowing the distinction between the main fracture paths: intracapsular and extracapsular fractures. It is worth noting the prediction of different fracture patterns for the same loading conditions, as observed during experimental tests.

RESULTS AND CONCLUSIONS

The internal distribution of bone mineral density and femur geometry strongly influences the femur fracture morphology and fracture load. Experimental fracture paths have been analysed by means of micro-computed tomography allowing the comparison of predicted and experimental crack surfaces, confirming the good accuracy of the numerical model.

Hip load capacity cut-points for Astronaut Skeletal Health NASA Finite Element Strength Task Group Recommendations

Concerns raised at a 2010 Bone Summit held for National Aeronautics and Space Administration Johnson Space Center led experts in finite element (FE) modeling for hip fracture prediction to propose including hip load capacity in the standards for astronaut skeletal health. The current standards for bone are based upon areal bone mineral density (aBMD) measurements by dual X-ray absorptiometry (DXA) and an adaptation of aBMD cut-points for fragility fractures. Task Group members recommended (i) a minimum permissible outcome limit (POL) for post-mission hip bone load capacity, (ii) use of FE hip load capacity to further screen applicants to astronaut corps, (iii) a minimum pre-flight standard for a second long-duration mission, and (iv) a method for assessing which post-mission physical activities might increase an astronaut’s risk for fracture after return. QCT-FE models of eight astronaut were analyzed using nonlinear single-limb stance (NLS) and posterolateral fall (NLF) loading configurations. QCT data from the Age Gene/Environment Susceptibility (AGES) Reykjavik cohort and the Rochester Epidemiology Project were analyzed using identical modeling procedures. The 75^th percentile of NLS hip load capacity for fractured elderly males of the AGES cohort (9537N) was selected as a post-mission POL. The NLF model, in combination with a Probabilistic Risk Assessment tool, was used to assess the likelihood of exceeding the hip load capacity during post-flight activities. There was no recommendation to replace the current DXA-based standards. However, FE estimation of hip load capacity appeared more meaningful for younger, physically active astronauts and was recommended to supplement aBMD cut-points.

Micro-Finite Element Analysis of the Proximal Femur on the Basis of High-Resolution Magnetic Resonance Images

PURPOSE OF REVIEW

Hip fractures have catastrophic consequences. The purpose of this article is to review recent developments in high-resolution magnetic resonance imaging (MRI)-guided finite element analysis (FEA) of the hip as a means to determine subject-specific bone strength.

RECENT FINDINGS

Despite the ability of DXA to predict hip fracture, the majority of fractures occur in patients who do not have BMD T scores less than - 2.5. Therefore, without other detection methods, these individuals go undetected and untreated. Of methods available to image the hip, MRI is currently the only one capable of depicting bone microstructure in vivo. Availability of microstructural MRI allows generation of patient-specific micro-finite element models that can be used to simulate real-life loading conditions and determine bone strength. MRI-based FEA enables radiation-free approach to assess hip fracture strength. With further validation, this technique could become a potential clinical tool in managing hip fracture risk.

Factors associated with proximal femur fracture determined in a large cadaveric cohort

Many researchers have used cadaveric fracture tests to determine the relationship between proximal femur (hip) fracture strength and a multitude of possible explanatory variables, typically considered one or two at a time. These variables include subject-specific proximal femur variables such as femoral neck areal bone mineral density (aBMD), sex, age, and geometry, as well as physiological hip fracture event variables such as fall speed and angle of impact. However, to our knowledge, no study has included all of these variables simultaneously in the same experimental dataset. To address this gap, the present study simultaneously included all of these subject-specific and fracture event variables in multivariate models to understand their contributions to femoral strength and fracture type. The primary aim of this study was to determine not only whether each of these variables contributed to the prediction of femoral strength, but also to determine the relative importance of each variable in strength prediction. A secondary aim was to similarly characterize the importance of these variables for the prediction of fracture type. To accomplish these aims, we characterized 197 proximal femurs (covering a wide range of subject-specific variables) with DXA and CT scans, and then tested the femurs to fracture in a sideways fall on the hip configuration. Each femur was tested using one of three fall speed conditions and one of four angles of impact (bone orientations). During each test, we acquired measurements of relevant force and displacement data. We then reduced the test data to determine femoral strength, and we used post-fracture CT scans to classify the fracture type (e.g., trochanteric, cervical). Using these results, the explanatory variables were analyzed with mixed statistical models to explain variations in hip fracture strength and fracture type, respectively. Five explanatory variables were statistically significant in explaining the variability in femoral strength: aBMD, sex, age, fall speed, and neck-shaft angle (P ≤ 0.0135). These five variables, including significant interactions, explained 80% of the variability in hip fracture strength. Additionally, when only aBMD, sex, and age (P < 0.0001) were considered in the model, again including significant interactions, these three variables alone explained 79% of the variability in hip fracture strength. So while fall speed (P = 0.0135) and neck-shaft angle (P = 0.0041) were statistically significant, the inclusion of these variables did not appreciably improve the prediction of hip fracture strength compared to the model that considered only aBMD, sex and age. For the variables we included in this study, in the ranges we considered, our findings indicate that the clinically-available information of patient age, sex and aBMD are sufficient for femoral strength assessment. These findings also suggest that there is little value in the extra effort required to characterize the effect of femoral geometry on strength, or to account for the probabilistic nature of fall-related factors such as fall speed and angle of impact. For fracture type, the only explanatory variable found to be significant was aBMD (P ≤ 0.0099). We found that the odds of having intertrochanteric fractures increased by 47% when aBMD decreased by one standard deviation (0.2 g/cm2).

Mechanical Strength of the Proximal Femur After Arthroscopic Osteochondroplasty for Femoroacetabular Impingement: Finite Element Analysis and 3-Dimensional Image Analysis

PURPOSE

To examine the influence of femoral neck resection on the mechanical strength of the proximal femur in actual surgery.

METHODS

Eighteen subjects who received arthroscopic cam resection for cam-type femoroacetabular impingement (FAI) were included. Finite element analyses (FEAs) were performed to calculate changes in simulative fracture load between pre- and postoperative femur models. The finite element femur models were constructed from computed tomographic images; thus, the models represented the shape of the original femur, including the bone resection site. Three-dimensional image analysis of the bone resection site was performed to identify morphometric factors that affect strength in the postoperative femur model. Four oblique sagittal planes running perpendicular to the femoral neck axis were used as reference planes to measure the bone resection site.

RESULTS

At the transcervical reference plane, both the bone resection depth and the cross-sectional area at the resection site correlated strongly with postoperative changes in the simulated fracture load (R² = 0.6, P = .0001). However, only resection depth was significantly correlated with the simulated fracture load at the reference plane for the head-neck junction. The resected bone volume did not correlate with the postoperative changes in the simulated fracture load.

CONCLUSIONS

The results of our FEA suggest that the bone resection depth measured at the head-neck junction and transcervical reference plane correlates with fracture risk after osteochondroplasty. By contrast, bone resection at more proximal areas did not have a significant effect on the postoperative femur model strength in our FEA. The total volume of resected bone was also not significantly correlated with postoperative changes in femur model strength.

CLINICAL RELEVANCE

This biomechanical study using FEA suggest that there is a risk of femoral neck fracture after arthroscopic cam resection, particularly when the resected lesion is located distally.

Finite element analysis of the femoral diaphysis of fresh-frozen cadavers with computed tomography and mechanical testing

Background

The purpose of this study was to validate a diaphyseal femoral fracture model using a finite element analysis (FEA) with mechanical testing in fresh-frozen cadavers.

Methods

We used 18 intact femora (9 right and 9 left) from 9 fresh-frozen cadavers. Specimens were obtained from 5 males and 4 females with a mean age of 85.6 years. We compared a computed tomography (CT)-based FEA model to diaphyseal femoral fracture loads and stiffness obtained by three-point bending. Four material characteristic conversion equations (the Keyak, Carter, and Keller equations plus Keller’s equation for the vertebra) with different shell thicknesses (0.3, 0.4, and 0.5 mm) were compared with the mechanical testing.

Results

The average fracture load was 4582.8 N and the mean stiffness was 942.0 N/mm from actual mechanical testing. FEA prediction using Keller’s equation for the vertebra with a 0.4-mm shell thickness showed the best correlations with the fracture load (R² = 0.76) and stiffness (R² = 0.54). Shell thicknesses of 0.3 and 0.5 mm in Keller’s equation for the vertebra also showed a strong correlation with fracture load (R² = 0.66 for both) and stiffness (R² = 0.50 and 0.52, respectively). There were no significant correlations with the other equations.

Conclusion

We validated femoral diaphyseal fracture loads and stiffness using an FEA in a cadaveric study.

Development of Subject-Specific Proximal Femur Finite Element Models Of Older Adults with Obesity to Evaluate the Effects of Weight Loss on Bone Strength

Study background

Recommendation of intentional weight loss in older adults remains controversial, due in part to the loss of bone mineral density (BMD) known to accompany weight loss. While finite element (FE) models have been used to assess bone strength, these methods have not been used to study the effects of weight loss. The purpose of this study is to develop subject-specific FE models of the proximal femur and study the effect of intentional weight loss on bone strength.

Methods

Computed tomography (CT) scans of the proximal femur of 25 overweight and obese (mean BMI=29.7 ± 4.0 kg/m²), older adults (mean age=65.6 ± 4.1 years) undergoing an 18-month intentional weight loss intervention were obtained at baseline and post-intervention. Measures of volumetric BMD (vBMD) and variable cortical thickness were derived from each subject CT scan and directly mapped to baseline and post-intervention models. Subject-specific FE models were developed using morphing techniques. Bone strength was estimated through simulation of a single-limb stance and sideways fall configuration.

Results

After weight loss intervention, there were significant decreases from baseline to 18 months in vBMD (total hip: -0.024 ± 0.013 g/cm³; femoral neck: -0.012 ± 0.014 g/cm³), cortical thickness (total hip: -0.044 ± 0.032 mm; femoral neck: -0.026 ± 0.039 mm), and estimated strength (stance: -0.15 ± 0.12 kN; fall: -0.04 ± 0.06 kN). Adjusting for baseline bone measures, body mass, and gender, correlations were found between weight change and change in total hip and femoral neck cortical thickness (all p<0.05). For every 1 kilogram of body mass lost cortical thickness in the total hip and femoral neck decreased by 0.003 mm and 0.004 mm, respectively. No significant correlations were present for the vBMD or strength data.

Conclusion

The developed subject-specific FE models could be used to better understand the effects of intentional weight loss on bone health.

Clinical hip fracture is accompanied by compression induced failure in the superior cortex of the femoral neck

Hip fractures pose a major health problem throughout the world due to their devastating impact. Current theories for why these injuries are so prevalent in the elderly point to an increased propensity to fall and decreases in bone mass with ageing. However, the fracture mechanisms, particularly the stress and strain conditions leading to bone failure at the hip remain unclear. Here, we directly examined the cortical bone from clinical intra-capsular hip fractures at a microscopic level, and found strong evidence of compression induced failure in the superior cortex. A total of 143 sections obtained from 24 femoral neck samples that were retrieved from 24 fracturing patients at surgery were examined using laser scanning confocal microscopy (LSCM) after fluorescein staining. The stained microcracks showed significantly higher density in the superior cortex than in the inferior cortex, indicating a greater magnitude of strain in the superior femoral neck during the failure-associated deformation and fracture process. The predominant stress state for each section was reconstructed based on the unique correlation between the microcrack pattern and the stress state. Specifically, we found clear evidence of longitudinal compression and buckling as the primary failure mechanisms in the superior cortex. These findings demonstrate the importance of microcrack analysis in studying clinical hip fractures, and point to the central role of the superior cortex failure as an important aspect of the failure initiation in clinical intra-capsular hip fractures.

Femoral fracture type can be predicted from femoral structure: A finite element study validated by digital volume correlation experiments

Proximal femoral fractures can be categorized into two main types: Neck and intertrochanteric fractures accounting for 53% and 43% of all proximal femoral fractures, respectively. The possibility to predict the type of fracture a specific patient is predisposed to would allow drug and exercise therapies, hip protector design, and prophylactic surgery to be better targeted for this patient rendering fracture preventing strategies more effective. This study hypothesized that the type of fracture is closely related to the patient-specific femoral structure and predictable by finite element (FE) methods. Fourteen femora were DXA scanned, CT scanned, and mechanically tested to fracture. FE-predicted fracture patterns were compared to experimentally observed fracture patterns. Measurements of strain patterns to explain neck and intertrochanteric fracture patterns were performed using a digital volume correlation (DVC) technique and compared to FE-predicted strains and experimentally observed fracture patterns. Although loaded identically, the femora exhibited different fracture types (six neck and eight intertrochanteric fractures). CT-based FE models matched the experimental observations well (86%) demonstrating that the fracture type can be predicted. DVC-measured and FE-predicted strains showed obvious consistency. Neither DXA-based BMD nor any morphologic characteristics such as neck diameter, femoral neck length, or neck shaft angle were associated with fracture type. In conclusion, patient-specific femoral structure correlates with fracture type and FE analyses were able to predict these fracture types. Also, the demonstration of FE and DVC as metrics of the strains in bones may be of substantial clinical value, informing treatment strategies and device selection and design. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:993-1001, 2018.

Prediction of fracture load and stiffness of the proximal femur by CT-based specimen specific finite element analysis: cadaveric validation study

Background

Finite element analysis (FEA) of the proximal femur has been previously validated with large mesh size, but these were insufficient to simulate the model with small implants in recent studies. This study aimed to validate the proximal femoral computed tomography (CT)-based specimen-specific FEA model with smaller mesh size using fresh frozen cadavers.

Methods

Twenty proximal femora from 10 cadavers (mean age, 87.1 years) were examined. CT was performed on all specimens with a calibration phantom. Nonlinear FEA prediction with stance configuration was performed using Mechanical Finder (mesh,1.5 mm tetrahedral elements; shell thickness, 0.2 mm; Poisson’s coefficient, 0.3), in comparison with mechanical testing. Force was applied at a fixed vertical displacement rate, and the magnitude of the applied load and displacement were continuously recorded. The fracture load and stiffness were calculated from force–displacement curve, and the correlation between mechanical testing and FEA prediction was examined.

Results

A pilot study with one femur revealed that the equations proposed by Keller for vertebra were the most reproducible for calculating Young’s modulus and the yield stress of elements of the proximal femur. There was a good linear correlation between fracture loads of mechanical testing and FEA prediction (R² = 0.6187) and between the stiffness of mechanical testing and FEA prediction (R² = 0.5499). There was a good linear correlation between fracture load and stiffness (R² = 0.6345) in mechanical testing and an excellent correlation between these (R² = 0.9240) in FEA prediction.

Conclusions

CT-based specimen-specific FEA model of the proximal femur with small element size was validated using fresh frozen cadavers. The equations proposed by Keller for vertebra were found to be the most reproducible for the proximal femur in elderly people.

Identification of the critical level of implantation of an osseointegrated prosthesis for above-knee amputees

The aim of our study was to identify potential critical levels of implantation of an osseointegrated prosthesis for above-knee amputees. The implant used was the OPRA system. It was inserted in the femur at four different amputation heights, characterized by their residual limb ratios (0.299, 0.44, 0.58 and 0.73). The stress and strain distribution was evaluated in the bone-implant system during walking, considering a body mass of 100 kg. Considerably high stimulus (11,489 με) in the tissue near the tip was found at the highest implantation level. All models presented small non-physiologic stress values in the tissue around the implant. The results revealed that the implantation level has a decisive effect on bone-implant performance. Mainly, the analysis indicates adverse biomechanical conditions for implantations in very short residual limbs.

Are DXA/aBMD and QCT/FEA Stiffness and Strength Estimates Sensitive to Sex and Age?

Dual X-ray absorptiometry (DXA) measures areal bone mineral density (aBMD) by simplifying a complex 3D bone structure to a 2D projection and is not equally effective for explaining fracture strength in women and men. Unlike DXA, subject-specific quantitative computed tomography-based finite element analysis (QCT/FEA) estimates fracture strength using 3D bone mineral distribution and geometry. By using experimentally-measured femoral stiffness and strength from a one hundred sample cadaveric cohort that included variations in sex and age, we wanted to determine if QCT/FEA estimates were able to better predict the experimental variations than DXA/aBMD. For each femur, DXA/aBMD was assessed and a QCT/FEA model was developed to estimate femoral stiffness and strength. Then, the femur was mechanically tested to fracture in a sideways fall on the hip position to measure stiffness and strength. DXA/aBMD and QCT/FEA estimates were compared for their sensitivity to sex and age with multivariate statistical analyses. When comparing the measured data with DXA/aBMD predictions, both age and sex were significant (p ≤ 0.0398) for both femoral stiffness and strength. However, QCT/FEA predictions of stiffness and strength showed sex was insignificant (p ≥ 0.23). Age was still significant (p ≤ 0.0072). These results indicate that QCT/FEA, unlike DXA/aBMD, accounted for bone differences due to sex.

Method and Instrumented Fixture for Femoral Fracture Testing in a Sideways Fall-on-the-Hip Position

Mechanical testing of femora brings valuable insights into understanding the contribution of clinically-measureable variables such as bone mineral density distribution and geometry on the femoral mechanical properties. Currently, there is no standard protocol for mechanical testing of such geometrically complex bones to measure strength, and stiffness. To address this gap we have developed a protocol to test cadaveric femora to fracture and to measure their biomechanical parameters. This protocol describes a set of adaptable fixtures to accommodate the various load magnitudes and directions accounting for possible bone orientations in a fall on the hip configuration, test speed, bone size, and left leg-right leg variations. The femora were prepared for testing by cleaning, cutting, scanning, and potting the distal end and greater trochanter contact surfaces in poly(methyl methacrylate) (PMMA) as presented in a different protocol. The prepared specimens were placed in the testing fixture in a position mimicking a sideways fall on the hip and loaded to fracture. During testing, two load cells measured vertical forces applied to the femoral head and greater trochanter, a six-axis load cell measured forces and moments at the distal femoral shaft, and a displacement sensor measured differential displacement between the femoral head and trochanter contact supports. High speed video cameras were used to synchronously record the sequence of fracture events during testing. The reduction of this data allowed us to characterize the strength, stiffness, and fracture energy for nearly 200 osteoporotic, osteopenic, and normal cadaveric femora for further development of engineering-based diagnostic tools for osteoporosis research.

Influence of Regional Difference in Bone Mineral Density on Hip Fracture Site in Elderly Females by Finite Element Analysis

Hip fracture is a kind of osteoporotic fractures in elderly patients. Its important monitoring indicator is to measure bone mineral density (BMD) using DXA. The stress characteristics and material distribution in different parts of the bones can be well simulated by three-dimensional finite element analysis. Our previous studies have demonstrated a linear positive correlation between clinical BMD and the density of three-dimensional finite element model of the femur. However, the correlation between the density variation between intertrochanteric region and collum femoris region of the model and the fracture site has not been studied yet. The present study intends to investigate whether the regional difference in the density of three-dimensional finite element model of the femur can be used to predict hip fracture site in elderly females. The CT data of both hip joints were collected from 16 cases of elderly female patients with hip fractures. Mimics 15.01 software was used to reconstruct the model of proximal femur on the healthy side. Ten kinds of material properties were assigned. In Abaqus 6.12 software, the collum femoris region and intertrochanteric region were, respectively, drawn for calculating the corresponding regional density of the model, followed by prediction of hip fracture site and final comparison with factual fracture site. The intertrochanteric region/collum femoris region density was [(1.20 ± 0.02) × 10(6)] on the fracture site and [(1.22 ± 0.03) × 10(6)] on the non-fracture site, and the difference was statistically significant (P = 0.03). Among 16 established models of proximal femur on the healthy side, 14 models were consistent with the actual fracture sites, one model was inconsistent, and one model was unpredictable, with the coincidence rate of 87.5 %. The intertrochanteric region or collum femoris region with lower BMD is more prone to hip fracture of the type on the corresponding site.

Exercise loading history and femoral neck strength in a sideways fall: A three-dimensional finite element modeling study

Over 90% of hip fractures are caused by falls. Due to a fall-induced impact on the greater trochanter, the posterior part of the thin superolateral cortex of the femoral neck is known to experience the highest stress, making it a fracture-prone region. Cortical geometry of the proximal femur, in turn, reflects a mechanically appropriate form with respect to habitual exercise loading. In this finite element (FE) modeling study, we investigated whether specific exercise loading history is associated with femoral neck structural strength and estimated fall-induced stresses along the femoral neck. One hundred and eleven three-dimensional (3D) proximal femur FE models for a sideways falling situation were constructed from magnetic resonance (MR) images of 91 female athletes (aged 24.7±6.1years, >8years competitive career) and 20 non-competitive habitually active women (aged 23.7±3.8years) that served as a control group. The athletes were divided into five distinct groups based on the typical loading pattern of their sports: high-impact (H-I: triple-jumpers and high-jumpers), odd-impact (O-I: soccer and squash players), high-magnitude (H-M: power-lifters), repetitive-impact (R-I: endurance runners), and repetitive non-impact (R-NI: swimmers). The von Mises stresses obtained from the FE models were used to estimate mean fall-induced stresses in eight anatomical octants of the cortical bone cross-sections at the proximal, middle, and distal sites along the femoral neck axis. Significantly (p<0.05) lower stresses compared to the control group were observed: the H-I group - in the superoposterior (10%) and posterior (19%) octants at the middle site, and in the superoposterior (13%) and posterior (22%) octants at the distal site; the O-I group - in the superior (16%), superoposterior (16%), and posterior (12%) octants at the middle site, and in the superoposterior (14%) octant at the distal site; the H-M group - in the superior (13%) and superoposterior (15%) octants at the middle site, and a trend (p=0.07, 9%) in the superoposterior octant at the distal site; the R-I group - in the superior (14%), superoposterior (23%) and posterior (22%) octants at the middle site, and in the superoposterior (19%) and posterior (20%) octants at the distal site. The R-NI group did not differ significantly from the control group. These results suggest that exercise loading history comprising various impacts in particular is associated with a stronger femoral neck in a falling situation and may have potential to reduce hip fragility.

Study of stress variations in single-stance and sideways fall using image-based finite element analysis

Image-based finite element analysis (FEA) has been considered an effective computational tool to predict hip fracture risk. The patient specific FEA gives an insight into the inclusive effect of three-dimensional (3D) complex bone geometry, and the distribution of inhomogeneous isotropic material properties in conjunction with loading conditions. The neck region of a femur is primarily the weakest in which fracture is likely to happen, when someone falls. A sideways fall results in the development of greater tensile and compressive stresses, respectively, in the inferior and superior aspects of the femoral neck, whereas the state of stress is reversed in usual gait or stance configuration. Herein, the variations of stresses have been investigated at the femoral neck region considering both single-stance and sideways fall. Finite element models of ten human femora have been generated using Quantitative Computed Tomography (QCT) scan datasets and have been simulated with an equal magnitude of load applied to the aforementioned configurations. Fracture risk indicator, defined as the ratio of the maximum compressive or tensile stress computed at the superior and inferior surfaces to the corresponding yield stress, has been used in this work to measure the variations of fracture risk between single-stance and sideways fall. The average variations of the fracture risk indicators between the fall and stance are at least 24.3% and 8% at the superior and inferior surfaces, respectively. The differences may interpret why sideways fall is more dangerous for the elderly people, causing hip fracture.

ASSOCIATION BETWEEN COMPUTED TOMOGRAPHIC CHARACTERISTICS AND FRACTURES FOLLOWING STEREOTACTIC RADIOSURGERY IN DOGS WITH APPENDICULAR OSTEOSARCOMA

The objective of this observational, descriptive, retrospective study was to report CT characteristics associated with fractures following stereotactic radiosurgery in canine patients with appendicular osteosarcoma. Medical records (1999 and 2012) of dogs that had a diagnosis of appendicular osteosarcoma and undergone stereotactic radiosurgery were reviewed. Dogs were included in the study if they had undergone stereotactic radiosurgery for an aggressive bone lesion with follow-up information regarding fracture status, toxicity, and date and cause of death. Computed tomography details, staging, chemotherapy, toxicity, fracture status and survival data were recorded. Overall median survival time (MST) and fracture rates of treated dogs were calculated. CT characteristics were evaluated for association with time to fracture. Forty-six dogs met inclusion criteria. The median overall survival time was 9.7 months (95% CI: 6.9-14.3 months). The fracture-free rates at 3, 6, and 9 months were 73%, 44%, and 38% (95% CI: 60-86%, 29-60%, and 22-54%), respectively. The region of bone affected was significantly associated with time to fracture. The median time to fracture was 4.2 months in dogs with subchondral bone involvement and 16.3 months in dogs without subchondral bone involvement (P-value = 0.027, log-rank test). Acute and late skin effects were present in 58% and 16% of patients, respectively. Findings demonstrated a need for improved patient selection for this procedure, which can be aided by CT-based prognostic factors to predict the likelihood of fracture.

Subject specific finite element modeling of periprosthetic femoral fracture using element deactivation to simulate bone failure

Subject-specific finite element (FE) modeling methodology could predict peri-prosthetic femoral fracture (PFF) for cementless hip arthoplasty in the early postoperative period. This study develops methodology for subject-specific finite element modeling by using the element deactivation technique to simulate bone failure and validate with experimental testing, thereby predicting peri-prosthetic femoral fracture in the early postoperative period. Material assignments for biphasic and triphasic models were undertaken. Failure modeling with the element deactivation feature available in ABAQUS 6.9 was used to simulate a crack initiation and propagation in the bony tissue based upon a threshold of fracture strain. The crack mode for the biphasic models was very similar to the experimental testing crack mode, with a similar shape and path of the crack. The fracture load is sensitive to the friction coefficient at the implant-bony interface. The development of a novel technique to simulate bone failure by element deactivation of subject-specific finite element models could aid prediction of fracture load in addition to fracture risk characterization for PFF.

Incorporating in vivo fall assessments in the simulation of femoral fractures with finite element models

Femoral fractures are a major health issue. Most experimental and finite element (FE) fracture studies use polymethylmethacrylate cups on the greater trochanter (GT) to simulate fall impact loads. However, in vivo fall studies showed that the femur is loaded distally from the GT. Our objective was to incorporate in vivo fall data in FE models to determine the effects of loading position and direction, and size of simulated impact site on the fracture load and fracture type for a healthy and an osteoporotic femur. Twelve sets of loading position and angles were applied through 'near point loads' on the models. Additional simulations were performed with 'cup loads' on the GT, similar to the literature. The results showed no significant difference between fracture loads from simulations with near point loads distally from the GT and those with cup loads on the GT. However, simulated fracture types differed, as near point loads distally from the GT generally resulted in various neck fractures, whilst cup load simulations predicted superior neck and trochanteric fractures only. This study showed that incorporating in vivo fall assessments in FE models by loading the models distally from the GT results in prediction of realistic fracture loads and fracture types.

Study of DXA-derived lateral-medial cortical bone thickness in assessing hip fracture risk

The currently available clinical tools have limited accuracy in predicting hip fracture risk in individuals. We investigated the possibility of using normalized cortical bone thickness (NCBT) estimated from the patient's hip DXA (dual energy X-ray absorptiometry) as an alternative predictor of hip fracture risk. Hip fracture risk index (HFRI) derived from subject-specific DXA-based finite element model was used as a guideline in constructing the mathematical expression of NCBT. We hypothesized that if NCBT has stronger correlations with HFRI than the single risk factors such as areal BMD (aBMD), then NCBT can be a better predictor. The hypothesis was studied using 210 clinical cases, including 60 hip fracture cases, obtained from the Manitoba Bone Mineral Density Database. The results showed that, in general HFRI has much stronger correlations with NCBT than any of the single risk factors; the strongest correlation was observed at the superior side of the narrowest femoral neck with r² = 0.81 (p < 0.001), which is much higher than the correlation with femoral aBMD, r² = 0.50 (p < 0.001). The capability of aBMD, NCBT, and HFRI in discriminating the hip fracture cases from the non-fracture ones, expressed as the area under the curve with 95% confidence interval, AUC (95% CI), is respectively 0.627 (0.593–0.657), 0.714 (0.644–0.784) and 0.839 (0.787–0.892). The short-term repeatability of aBMD, NCBT, and HFRI, measured by the coefficient of variation (CV, %), was found to be in the range of (0.64–1.22), (1.93–3.41), (3.10–4.16), respectively. We thus concluded that NCBT is potentially a better predictor of hip fracture risk.

•

A new algorithm developed for estimating cortical bone thickness from clinical DXA images.

•

Validation of the algorithm with QCT Pro, a commercial software for processing QCT scans.

•

A more effective risk predictor constructed from finite element studies.

•

The risk predictor can be readily integrated into the current clinical procedure.