Sideways Fall Research Articles

Osteoporotic Hip Fracture Prediction: Is T-Score-Based Criterion Enough? A Hip Structural Analysis-Based Model.

At present, the current gold-standard for osteoporosis diagnosis is based on bone mineral density (BMD) measurement, which, however, has been demonstrated to poorly estimate fracture risk. Further parameters in the hands of the clinicians are represented by the hip structural analysis (HSA) variables, which include geometric information of the proximal femur cross section. The purpose of this study was to investigate the suitability of HSA parameters as additional hip fracture risk predictors. With this aim, twenty-eight three-dimensional patient-specific models of the proximal femur were built from computed tomography (CT) images and a sideways fall condition was reproduced by finite element (FE) analyses. A tensile or compressive predominance based on minimum and maximum principal strains was determined at each volume element and a risk factor (RF) was calculated. The power of HSA variables combinations to predict the maximum superficial RF values was assessed by multivariate linear regression analysis. The optimal regression model, identified through the Akaike information criterion (AIC), only comprises two variables: the buckling ratio (BR) and the neck-shaft angle (NSA). In order to validate the study, the model was tested on two additional patients who suffered a hip fracture after a fall. The results classified the patients in the high risk level, confirming the prediction power of the adopted model.

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

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

Ageing Effects on 3-Dimensional Femoral Neck Cross-Sectional Asymmetry: Implications for Age-Related Bone Fragility in Falling

This paper explores the effects of aging on femoral neck (FN) anatomy in a study of women aged 20–90years in relation to implications for FN fracture propensity in buckling. Five hundred and four participants were scanned by Quantitative Computed Tomography and analyzed using Quantitative Computed Tomography Pro BIT (Mindways). FN cross-section was split through geometric center into superior and inferior sectors. Bone mass, structural measurements, and bone mineral density were analyzed. Buckling ratio was calculated as ratio of buckling radius to cortical thickness. Between 2nd decade and 8th decade, age-related integral bone mass reduction in superior sector was substantially larger than in inferior sector (33% compared to 21%), especially in cortical bone superiorly compared to inferiorly (53% vs 21%; p < 0.001), principally due to reduction in cortical thickness, averaged cortical thickness (56%) with little difference in density. Superior and inferior sector trabecular bone mineral density reduction was similar at 41% and 43% respectively. Differential cortical bone loss in superior sector resulted in a 59% inferior displacement (δ) of center-of-mass from geometric center. Differences in δ and averaged cortical thickness with age accounted for a 151% increase in mean superior buckling ratio from 9 to 23. Analysis confirms significant progressive age-related superior cortical bone loss as the major age effect on FN structure with relative preservation of inferior cortex probably related to maintenance of inferior sector by regular loading as a result of standing and walking. Computation of buckling ratio may allow prediction of fracture propensity in a sideways fall.

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).

A novel sideways fall simulator to study hip fractures ex vivo.

Falls to the side are the leading cause of hip fractures in the elderly. The load that a person experiences during a fall cannot be measured with volunteers for ethical reasons. To evaluate injurious loads, while considering relevant energy input and body posture for a sideways fall, a subject-specific cadaveric impact experiment was developed. Full cadaveric femur-pelvis constructs (N = 2) were embedded in surrogate soft tissue material and attached to metallic surrogate lower limbs. The specimens were then subjected to an inverted pendulum motion, simulating a fall to the side with an impact to the greater trochanter. The load at the ground and the deformation of the pelvis were evaluated using a 6-axis force transducer and two high-speed cameras. Post-test, a trauma surgeon (PG) evaluated specimen injuries. Peak ground contact forces were 7132 N and 5641 N for the fractured and non-fractured specimen, respectively. We observed a cervical fracture of the femur in one specimen and no injuries in a second specimen, showing that the developed protocol can be used to differentiate between specimens at high and low fracture risk.

Mechanical Loading of the Femoral Neck in Human Locomotion.

Advancing age and reduced loading are associated with a reduction in bone formation. Conversely, loading increases periosteal apposition and may reduce remodeling imbalance and slow age-related bone loss, an important outcome for the proximal femur, which is a common site of fracture. The ability to take advantage of bone's adaptive response to increase bone strength has been hampered by a lack of knowledge of which exercises and specific leg muscles load the superior femoral neck: a common region of microcrack initiation and progression following a sideways fall. We used an in vivo method of quantifying focal strains within the femoral neck in postmenopausal women during walking, stair ambulation, and jumping. Relative to walking, stair ambulation and jumping induced significantly higher strains in the anterior and superior aspects of the femoral neck, common regions of microcrack initiation and progression following a fall. The gluteus maximus, a hip extensor muscle, induced strains in the femoral neck during stair ambulation and jumping, in contrast to walking which induced strains via the iliopsoas, a hip flexor. The ground reaction force was closely associated with the level of strain during each task, providing a surrogate indicator of the potential for a given exercise to load the femoral neck. The gluteal muscles combined with an increased ground reaction force relative to walking induce high focal strains within the anterosuperior region of the femoral neck and therefore provide a target for exercise regimens designed to slow bone loss and maintain or improve microstructural strength. Model files used for calculating femoral neck strains are available at uitbl.mechse.illinois.edu/downloads © 2018 American Society for Bone and Mineral Research.

Ricci-flow based conformal mapping of the proximal femur to identify exercise loading effects

The causal relationship between habitual loading and adaptive response in bone morphology is commonly explored by analysing the spatial distribution of mechanically relevant features. In this study, 3D distribution of features in the proximal femur of 91 female athletes (5 exercise loading groups representing habitual loading) is contrasted with 20 controls. A femur specific Ricci-flow based conformal mapping procedure was developed for establishing correspondence among the periosteal surfaces. The procedure leverages the invariance of the conformal mapping method to isometric shape differences to align surfaces in the 2D parametric domain, to produce dense correspondences across an isotopological set of surfaces. This is implemented through a multi-parametrisation approach to detect surface features and to overcome the issue of inconsistency in the anatomical extent present in the data. Subsequently, the group-wise distribution of two mechanically relevant features was studied – cortical thickness and surface principal strains (simulation results of a sideways fall). Statistical inferences over the surfaces were made by contrasting the athlete groups with the controls through statistical parametric mapping. With the aid of group-wise and composite-group maps, proximal femur regions affected by specific loading groups were identified with a high degree of spatial localisation.

Comparison of femoral strength and fracture risk index derived from DXA-based finite element analysis for stratifying hip fracture risk: A cross-sectional study.

Dual-energy X-ray absorptiometry (DXA)-based finite element analysis (FEA) has been studied for assessment of hip fracture risk. Femoral strength (FS) is the maximum force that the femur can sustain before its weakest region reaches the yielding limit. Fracture risk index (FRI), which also considers subject-specific impact force, is defined as the ratio of von Mises stress induced by a sideways fall to the bone yield stress over the proximal femur. We compared risk stratification for prior hip fracture using FS and FRI derived from DXA-based FEA. The study cohort included women aged ≥65years undergoing baseline hip DXA, with femoral neck T-scores <-1 and no osteoporosis treatment; 324 cases had prior hip fracture and 655 controls had no prior fracture. Using anonymized DXA hip scans, we measured FS and FRI. Separate multivariable logistic regression models were used to estimate odds ratios (ORs), c-statistics and their 95% confidence intervals (95% CIs) for the association of hip fracture with FS and FRI. Increased hip fracture risk was associated with lower FS (OR per SD 1.36, 95% CI: 1.15, 1.62) and higher FRI (OR per SD 1.99, 95% CI: 1.63, 2.43) after adjusting for Fracture Risk Assessment Tool (FRAX) hip fracture probability computed with bone mineral density (BMD). The c-statistic for the model containing FS (0.69; 95% CI: 0.65, 0.72) was lower than the c-statistic for the model with FRI (0.77; 95% CI: 0.74, 0.80) or femoral neck BMD (0.74; 95% CI: 0.71, 0.77; all P<0.05). FS and FRI were independently associated with hip fracture, but there were differences in performance characteristics.

On the Failure Initiation in the Proximal Human Femur Under Simulated Sideways Fall.

The limitations of areal bone mineral density measurements for identifying at-risk individuals have led to the development of alternative screening methods for hip fracture risk including the use of geometrical measurements from the proximal femur and subject specific finite element analysis (FEA) for predicting femoral strength, based on quantitative CT data (qCT). However, these methods need more development to gain widespread clinical applications. This study had three aims: To investigate whether proximal femur geometrical parameters correlate with obtained femur peak force during the impact testing; to examine whether or not failure of the proximal femur initiates in the cancellous (trabecular) bone; and finally, to examine whether or not surface fracture initiates in the places where holes perforate the cortex of the proximal femur. We found that cortical thickness around the trochanteric-fossa is significantly correlated to the peak force obtained from simulated sideways falling (R 2=0.69) more so than femoral neck cortical thickness (R 2=0.15). Dynamic macro level FE simulations predicted that fracture generally initiates in the cancellous bone compartments. Moreover, our micro level FEA results indicated that surface holes may be involved in primary failure events.

Association of incident hip fracture with the estimated femoral strength by finite element analysis of DXA scans in the Osteoporotic Fractures in Men (MrOS) study.

Bone fractures occur only when it is loaded beyond its ultimate strength. The goal of this study was to determine the association of femoral strength, as estimated by finite element (FE) analysis of DXA scans, with incident hip fracture as a single condition or with femoral neck (FN) and trochanter (TR) fractures separately in older men. This prospective case-cohort study included 91 FN and 64 TR fracture cases and a random sample of 500 men (14 had a hip fracture) from the Osteoporotic Fractures in Men study during a mean ± SD follow-up of 7.7 ± 2.2years. We analysed the baseline DXA scans of the hip using a validated plane-stress, linear-elastic FE model of the proximal femur and estimated the femoral strength during a sideways fall. The estimated strength was significantly (P < 0.05) associated with hip fracture independent of the TR and total hip (TH) BMDs but not FN BMD, and combining the strength with BMD did not improve the hip fracture prediction. The strength estimate was associated with FN fractures independent of the FN, TR and TH BMDs; the age-BMI-BMD adjusted hazard ratio (95% CI) per SD decrease of the strength was 1.68 (1.07-2.64), 2.38 (1.57, 3.61) and 2.04 (1.34, 3.11), respectively. This association with FN fracture was as strong as FN BMD (Harrell's C index for the strength 0.81 vs. FN BMD 0.81) and stronger than TR and TH BMDs (0.8 vs. 0.78 and 0.81 vs. 0.79). The strength's association with TR fracture was not independent of hip BMD. Although the strength estimate provided additional information over the hip BMDs, its improvement in predictive ability over the hip BMDs was confined to FN fracture only and limited.

Material mapping strategy to improve the predicted response of the proximal femur to a sideways fall impact.

Sideways falls are largely responsible for the highly prevalent osteoporotic hip fractures in today's society. These injuries are dynamic events, therefore dynamic FE models validated with dynamic ex vivo experiments provide a more realistic simulation than simple quasi-static analysis. Drop tower experiments using cadaveric specimens were used to identify the material mapping strategy that provided the most realistic mechanical response under impact loading. The present study tested the addition of compression-tension asymmetry, tensile bone damage, and cortical-specific strain rate dependency to the material mapping strategy of fifteen dynamic FE models of the proximal femur, and found improved correlations and reduced error for whole bone stiffness (R2 = 0.54, RSME = 0.87kN/mm) and absolute maximum force (R2 = 0.56, RSME =0.57kN), and a high correlation in impulse response (R2 = 0.82, RSME =12.38kg/s). Simulations using fully bonded nodes between the rigid bottom plate and PMMA cap supporting the femoral head had higher correlations and less error than simulations using a frictionless sliding at this contact surface. Strain rates over 100/s were observed in certain elements in the femoral neck and trochanter, indicating that additional research is required to better quantify the strain rate dependencies of both trabecular and cortical bone at these strain rates. These results represent the current benchmark in dynamic FE modeling of the proximal femur in sideways falls. Future work should also investigate improvements in experimental validation techniques by developing better displacement measurements and by enhancing the biofidelity of the impact loading wherever possible.

Automated DXA-based finite element analysis for hip fracture risk stratification: a cross-sectional study.

The study aims to examine the association between prior hip fracture and FRIs derived from automated finite element analysis (FEA) of DXA hip scans. Femoral neck, intertrochanteric, and subtrochanteric FRIs were calculated as the von Mises stress induced by a sideways fall divided by the bone yield stress over the specified region of interest (ROI). Using the Manitoba Bone Mineral Density Database, we selected women age ≥65years with femoral neck T-scores below -1 and no osteoporosis treatment. From this population, we identified 324 older women with hip fracture before DXA testing and a random sample of 658 non-fracture controls. FRIs were derived from the anonymized DXA scans. Logistic regression models were used to estimate odds ratios (ORs) and 95% confidence intervals (95% CIs) for the associations between FRIs (per SD increase) and hip fracture. After adjusting for FRAX score (hip fracture with BMD), femoral neck FRI (OR 1.36, 95% CI 1.13, 1.64), intertrochanteric FRI (OR 1.81, 95% CI 1.44, 2.27), and subtrochanteric FRI (OR 2.09, 95% CI 1.68, 2.60) were associated with hip fracture. Intertrochanteric and subtrochanteric FRIs gave significantly higher c-statistics (all P≤0.05) than femoral neck BMD. Subgroup analyses showed that all FRIs were more strongly associated with hip fracture in women who were younger and had higher body mass index (BMI) or non-osteoporotic BMD (all P interaction<0.1). FRIs derived from DXA-based FEA were independently associated with prior hip fracture, suggesting that they could potentially improve hip fracture risk assessment.

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.

Comparison of non-invasive assessments of strength of the proximal femur

It is not clear which non-invasive method is most effective for predicting strength of the proximal femur in those at highest risk of fracture. The primary aim of this study was to compare the abilities of dual energy X-ray absorptiometry (DXA)-derived aBMD, quantitative computed tomography (QCT)-derived density and volume measures, and finite element analysis (FEA)-estimated strength to predict femoral failure load. We also evaluated the contribution of cortical and trabecular bone measurements to proximal femur strength. We obtained 76 human cadaveric proximal femurs (50 women and 26 men; age 74±8.8years), performed imaging with DXA and QCT, and mechanically tested the femurs to failure in a sideways fall configuration at a high loading rate. Linear regression analysis was used to construct the predictive model between imaging outcomes and experimentally-measured femoral strength for each method. To compare the performance of each method we used 3-fold cross validation repeated 10 times. The bone strength estimated by QCT-based FEA predicted femoral failure load (R2adj=0.78, 95%CI 0.76–0.80; RMSE=896N, 95%CI 830–961) significantly better than femoral neck aBMD by DXA (R2adj=0.69, 95%CI 0.66–0.72; RMSE=1011N, 95%CI 952–1069) and the QCT-based model (R2adj=0.73, 95%CI 0.71–0.75; RMSE=932N, 95%CI 879–985). Both cortical and trabecular bone contribute to femoral strength, the contribution of cortical bone being higher in femurs with lower trabecular bone density. These findings have implications for optimizing clinical approaches to assess hip fracture risk. In addition, our findings provide new insights that will assist in interpretation of the effects of osteoporosis treatments that preferentially impact cortical versus trabecular bone.

New approaches for cement-based prophylactic augmentation of the osteoporotic proximal femur provide enhanced reinforcement as predicted by non-linear finite element simulations

BackgroundHigh incidence and increased mortality related to secondary, contralateral proximal femoral fractures may justify invasive prophylactic augmentation that reinforces the osteoporotic proximal femur to reduce fracture risk. Bone cement-based approaches (femoroplasty) may deliver the required strengthening effect; however, the significant variation in the results of previous studies calls for a systematic analysis and optimization of this method. Our hypothesis was that efficient generalized augmentation strategies can be identified via computational optimization. MethodsThis study investigated, by means of finite element analysis, the effect of cement location and volume on the biomechanical properties of fifteen proximal femora in sideways fall. Novel cement cloud locations were developed using the principles of bone remodeling and compared to the “single central” location that was previously reported to be optimal. FindingsThe new augmentation strategies provided significantly greater biomechanical benefits compared to the “single central” cement location. Augmenting with approximately 12ml of cement in the newly identified location achieved increases of 11% in stiffness, 64% in yield force, 156% in yield energy and 59% in maximum force, on average, compared to the non-augmented state. The weaker bones experienced a greater biomechanical benefit from augmentation than stronger bones. The effect of cement volume on the biomechanical properties was approximately linear. Results of the “single central” model showed good agreement with previous experimental studies. InterpretationThese findings indicate enhanced potential of cement-based prophylactic augmentation using the newly developed cementing strategy. Future studies should determine the required level of strengthening and confirm these numerical results experimentally.

側面方向の外乱に対するCapturabilityに基づく二足ロボットの転倒回避

側面方向の外乱に対するCapturabilityに基づく二足ロボットの転倒回避

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.

Can CT image deblurring improve finite element predictions at the proximal femur?

The aim of this study was to determine if a CT image deblurring algorithm can improve CT-based FE modelling accuracy at the proximal femur.Experimental data (CT scans of fourteen proximal fresh-frozen cadaveric femurs, non-destructive surface strain measurements in stance and sideways fall loading configurations on all femurs, and failure loads obtained in stance for seven specimens, in sideways fall for the other seven) were taken from a recent study (Schileo et al., 2014). An estimate of the 3D Point Spread Function for each CT scan was used within a deconvolution solver to perform the deblurring. The most proximal regions of three specimens were scanned using an HRpQCT scanner and compared to the original and deblurred CT images to quantify errors in bone contour estimates and determine correlation of intensity values within the bone contours. Subject-specific FE models of the proximal femur were generated. The accuracy of deblurred FE predictions against experimental measurements was compared to the published (non-deblurred) FE results.When compared to HRpQCT, CT deblurring led to lower mean surface distances (0.31 vs. 0.49mm) and higher CT intensity correlations with respect to the original CT. All indicators of strain prediction accuracy were significantly improved in deblurred FE models, more markedly at the femoral neck (peak error reduced by 38%). Failure load prediction, based on a simple elastic limit model, was also improved in deblurred FE models, although differently for stance and sideways fall loading conditions. In stance, correlation was unchanged, but specimen-wise errors were reduced (mean error 10% vs. 15%). In sideways fall, correlation notably increased (R2=0.95 vs. 0.81), despite a general overestimation of failure load.In summary, the proposed CT deblurring technique yielded moderate but significant improvements in FE predictions, and may thus be considered a first step toward the improvement of CT-based FE models of the human femur.

The influence of the modulus–density relationship and the material mapping method on the simulated mechanical response of the proximal femur in side-ways fall loading configuration

Contributing to slow advance of finite element (FE) simulations for hip fracture risk prediction, into clinical practice, could be a lack of consensus in the biomechanics community on how to map properties to the models. Thus, the aim of the present study was first, to systematically quantify the influence of the modulus–density relationship (E–ρ) and the material mapping method (MMM) on the predicted mechanical response of the proximal femur in a side-ways fall (SWF) loading configuration and second, to perform a model-to-model comparison of the predicted mechanical response within the femoral neck for all the specimens tested in the present study, using three different modelling techniques that have yielded good validation outcome in terms of surface strain prediction and whole bone response according to the literature. We found the outcome to be highly dependent on both the E–ρ relationship and the MMM. In addition, we found that the three modelling techniques that have resulted in good validation outcome in the literature yielded different principal strain prediction both on the surface as well as internally in the femoral neck region of the specimens modelled in the present study. We conclude that there exists a need to carry out a more comprehensive validation study for the SWF loading mode to identify which combination of MMMs and E–ρ relationship leads to the best match for whole bone and local mechanical response. The MMMs tested in the present study have been made publicly available at https://simtk.org/home/mitk-gem.

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

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