Can neck fractures in proximal humeri be predicted by CT-based FEA?

BACKGROUND: Proximal humeri fractures at anatomical and surgical neck (∼5% and ∼50% incidence respectively) are frequent in elderly population. Yet, neither in-vitro experiments nor CT-based finite element analyses (CTFEA) have investigated these in depth. Herein we enhance (Dahan et al., 2019) (addressing anatomical neck fractures) by more experiments and specimens, accounting for surgical neck fractures and explore CTFEA's prediction of humeri mechanical response and yield force. METHODS: Four fresh frozen human humeri were tested in a new experimental configuration inducing surgical neck fractures. Digital image correlation (DIC) provided strains and displacements on humeri surfaces and used to validate CTFEA predictions. CTFEA were enhanced herein to improve the accuracy at the proximal neck: A cortical bone mapping (CBM) algorithm was implemented to overcome insufficient scanning resolution, and a new trabecular material mapping was investigated. RESULTS: The new experimental setting induced impacted surgical neck fractures in all humeri. Excellent DIC to CTFEA correlation in strains was obtained at the shaft (slope 0.984, R2=0.99) and a fair agreement (slope 0.807, R2=0.73) at the neck. CBM algorithm had worsened the correlation, whereas the new material mapping had a negligible influence. Yield loads predictions improved considerably when trabecular yielding (maximum principal strain criterion) was considered instead of surface cortical yielding. DISCUSSION: CTFEA well predicts strains on the shaft and reasonably well on the neck. This enhances former conclusions by past studies conducted using SGs, now also evident by DIC. Yield load prediction for surgical neck fractures (involving crushing of trabecular bone) is predicted better by trabecular failure laws rather than cortex ones. Further FEA studies using trabecular orthotropic constitutive models and failure laws are warrant.

Finite element analyses for predicting anatomical neck fractures in the proximal humerus.

BACKGROUND: Proximal humerus fractures which occur as a result of a fall on an outstretched arm are frequent among the elderly population. The necessity of stabilizing such fractures by surgical procedures is a controversial matter among surgeons. Validating a personalized FE analysis by ex-vivo experiments of humeri and mimicking such fractures by experiments is the first step along the path to determine the necessity of such surgeries. METHODS: Four fresh frozen human humeri were loaded using a new simple experimental setting, so to fracture the humeri at the anatomical neck. Strains on humeri's surfaces predicted by the high order FE analyses (as in Dahan et al., 2016) were compared to the experimental observations to further enhance the validity of the FE analyses. A simplified yield criterion based on a linear elastic analysis and principal strains was used to predict the anatomical neck fracture as observed in the experiment. FINDINGS: An excellent correlation between experimental measured and FE predicted strains was obtained (slope of 0.99 and R2=0.98). All humeri were fractured at the anatomical neck. The predicted yield load was within 10%-20% accuracy. INTERPRETATION: High-order FE analyses reliably predict strains and yield loads in the humeri. Fractures induced by the experimental setting correspond to anatomical neck fractures noticed in practice and classified as AO C1.1-C1.3. Surgical neck fractures, which are most common in clinical practice, could not be realized in the proposed experiments, and a different experimental setting should be sought to obtain them ex-vivo.

A novel phase field method for modeling the fracture of long bones.

A proximal humerus fracture is an injury to the shoulder joint that necessitates medical attention. While it is one of the most common fracture injuries impacting the elder community and those who suffer from traumatic falls or forceful collisions, there are almost no validated computational methods that can accurately model these fractures. This could be due to the complex, inhomogeneous bone microstructure, complex geometries, and the limitations of current fracture mechanics methods. In this paper, we develop a novel phase field method to investigate the proximal humerus fracture. To model the fracture in the inhomogeneous domain, we propose a power-law relationship between bone mineral density and critical energy release rate. The method is validated by an in vitro experiment, in which a human humerus is constrained on both ends while subjected to compressive loads on its head, in the longitudinal direction, that lead to fracture at the anatomical neck. CT scans are employed to acquire the bone geometry and material parameters, from which detailed finite element meshes with inhomogeneous Young modulus distributions are generated. The numerical method, implemented in a high performance computing environment, is used to quantitatively predict the complex 3D brittle fracture of the bone and is shown to be in good agreement with experimental observations. Furthermore, our findings show that the damage is initiated in the trabecular bone-head and propagates outward towards the bone cortex. We conclude that the proposed phase field method is a promising approach to model bone fracture.

Scanner influence on the mechanical response of QCT-based finite element analysis of long bones.

Patient-specific QCT-based finite element (QCTFE) analyses enable highly accurate quantification of bone strength. We evaluated CT scanner influence on QCTFE models of long bones. A femur, humerus, and proximal femur without the head were scanned with K2HPO4 phantoms by seven CT scanners (four models) using typical clinical protocols. QCTFE models were constructed. The geometrical dimensions, as well as the QCT-values expressed in Hounsfield unit (HU) distribution was compared. Principal strains at representative regions of interest (ROIs), and maximum principal strains (associated with fracture risk) were compared. Intraclass correlation coefficients (ICCs) were calculated to evaluate strain prediction reliability for different scanners. Repeatability was examined by scanning the femur twice and comparing resulting QCTFE models. Maximum difference in geometry was 2.3%. HU histograms before phantom calibration showed wide variation between QCT scans; however, bone density histogram variability was reduced after calibration and algorithmic manipulation. Relative standard deviation (RSD) in principal strains at ROIs was <10.7%. ICC estimates between scanners were >0.9. Fracture-associated strain had 6.7%, 8.1%, and 13.3% maximum RSD for the femur, humerus, and proximal femur, respectively. The difference in maximum strain location was <2 mm. The average difference with repeat scans was 2.7%. Quantification of strain differences showed mean RSD bounded by ∼6% in ROIs. Fracture-associated strains in "regular" bones showed a mean RSD bounded by ∼8%. Strains were obtained within a ±10% difference relative to the mean; thus, in a longitudinal study only changes larger than 20% in the principal strains may be significant. ICCs indicated high reliability of QCTFE models derived from different scanners.

Verified and validated finite element analyses of humeri.

BACKGROUND: Although ~200,000 emergency room visits per year in the US alone are associated with fractures of the proximal humerus, only limited studies exist on their mechanical response. We hypothesise that for the proximal humeri (a) the mechanical response can be well predicted by using inhomogeneous isotropic material properties, (b) the relation between bone elastic modulus and ash density (E(ρash)) is similar for the humerus and the femur, and may be general for long bones, and (c) it is possible to replicate a proximal humerus fracture in vitro by applying uniaxial compression on humerus׳ head at a prescribed angle. METHODS: Four fresh frozen proximal humeri were CT-scanned, instrumented by strain-gauges and loaded at three inclination angles. Thereafter head displacement was applied to obtain a fracture. CT-based high order (p-) finite element (FE) and classical (h-) FE analyses were performed that mimic the experiments and predicted strains were compared to the experimental observations. RESULTS: The E(ρash) relationship appropriate for the femur is equally appropriate for the humeri: predicted strains in the elastic range showed an excellent agreement with experimental observations with a linear regression slope of m=1.09 and a coefficient of regression R(2)=0.98. p-FE and h-FE results were similar for the linear elastic response. Although fractures of the proximal humeri were realised in the in vitro experiments, the contact FE analyses (FEA) were unsuccessful in representing properly the experimental boundary conditions. DISCUSSION: The three hypotheses were confirmed and the linear elastic response of the proximal humerus, attributed to a stage at which the cortex bone is intact, was well predicted by the FEA. Due to a large post-elastic behaviour following the cortex fracture, a new non-linear constitutive model for proximal humerus needs to be incorporated into the FEA to well represent proximal humerus fractures. Thereafter, more in vitro experiments are to be performed, under boundary conditions that may be well represented by the FEA, to allow a reliable simulation of the fracture process.

Predicting the stiffness and strength of human femurs with real metastatic tumors.

BACKGROUND: Predicting patient specific risk of fracture in femurs with metastatic tumors and the need for surgical intervention are of major clinical importance. Recent patient-specific high-order finite element methods (p-FEMs) based on CT-scans demonstrated accurate results for healthy femurs, so that their application to metastatic affected femurs is considered herein. METHODS: Radiographs of fresh frozen proximal femur specimens from donors that died of cancer were examined, and seven pairs with metastatic tumor were identified. These were CT-scanned, instrumented by strain-gauges and loaded in stance position at three inclination angles. Finally the femurs were loaded until fracture that usually occurred at the neck. Histopathology was performed to determine whether metastatic tumors are present at fractured surfaces. Following each experiment p-FE models were created based on the CT-scans mimicking the mechanical experiments. The predicted displacements, strains and yield loads were compared to experimental observations. RESULTS: The predicted strains and displacements showed an excellent agreement with the experimental observations with a linear regression slope of 0.95 and a coefficient of regression R(2)=0.967. A good correlation was obtained between the predicted yield load and the experimental observed yield, with a linear regression slope of 0.80 and a coefficient of regression R(2)=0.78. DISCUSSION: CT-based patient-specific p-FE models of femurs with real metastatic tumors were demonstrated to predict the mechanical response very well. A simplified yield criterion based on the computation of principal strains was also demonstrated to predict the yield force in most of the cases, especially for femurs that failed at small loads. In view of the limited capabilities to predict risk of fracture in femurs with metastatic tumors used nowadays, the p-FE methodology validated herein may be very valuable in making clinical decisions.