A computed tomography-based limit analysis approach to investigate the mechanical behavior of the human femur prone to fracture

Abstract

The paper proposes a refined CT-based FE modelling strategy that implements a limit analysis numerical procedure, namely the Elastic Compensation Method (ECM), to estimate a lower bound to the collapse load of a human femur. In particular, the model geometry was obtained from CT images by segmentation of a fresh-frozen human cadaveric femur that was discretized with second-order tetrahedral 3D finite elements. A yield criterion of Tsai–Wu-type, expressed in principal stress space, was adopted to model the bone tissues for which the strength limit values in tension, compression and shear are computed locally from the femoral density distribution also derived from CT images. The developed CT-based numerical technique showed the ability to predict, at least for the examined femur for which the experimental collapse load is available, a lower bound to the collapse load. The proposed approach seems a promising and effective tool that could be adopted into clinical practice to predict the fracture risk of human femur starting from patient-specific data given by medical imaging.