Determination of the heterogeneous anisotropic elastic properties of human femoral bone: from nanoscopic to organ scale.

Cortical bone is a multiscale composite material. Its elastic properties are anisotropic and heterogeneous across its cross-section, due to endosteal bone resorption which might affect bone strength. The aim of this paper was to describe a homogenization method leading to the estimation of the variation of the elastic coefficients across the bone cross-section and along the bone longitudinal axis. The method uses the spatial variations of bone porosity and of the degree of mineralization of the bone matrix (DMB) obtained from the analysis of 3-D synchrotron micro-computed tomography images. For all three scales considered (the foam (100 nm), the ultrastructure (5 microm) and the mesoscale (500 microm)), the elastic coefficients were determined using the Eshelby's inclusion problem. DMB values were used at the scale of the foam. Collagen was introduced at the scale of the ultrastructure and bone porosity was introduced at the mesoscale. The pores were considered as parallel cylinders oriented along the bone axis. Each elastic coefficient was computed for different regions of interest, allowing an estimation of its variations across the bone cross-section and along the bone longitudinal axis. The method was applied to a human femoral neck bone specimen, which is a site of osteoporotic fracture. The computed elastic coefficients for cortical bone were in good agreement with experimental results, but some discrepancies were obtained in the endosteal part (trabecular bone). These results highlight the importance of accounting for the heterogeneity of cortical bone properties across bone cross-section and along bone longitudinal axis. 2010 Elsevier Ltd. All rights reserved.