Verified and validated finite element analyses of humeri

Gal Dahan, Nir Trabelsi, Ori Safran, Zohar Yosibash

Research output: Contribution to journalArticlepeer-review

25 Scopus citations


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 R2=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.

Original languageEnglish
Pages (from-to)1094-1102
Number of pages9
JournalJournal of Biomechanics
Issue number7
StatePublished - 3 May 2016
Externally publishedYes

Bibliographical note

Publisher Copyright:
© 2016 Elsevier Ltd.


The authors thank Prof. Charles Milgrom, from the Hadassah Hospital for helpful discussions and Mr. Ilan Gilad and Yekutiel Katz, from the Ben Gurion University for their assistance with the experiments. This study was made possible through the generous support of a (Milgrom Foundation for Science) grant.

FundersFunder number
Milgrom Foundation for Science


    • Finite element analysis
    • Humerus
    • Impacted-fracture


    Dive into the research topics of 'Verified and validated finite element analyses of humeri'. Together they form a unique fingerprint.

    Cite this