Determination of anisotropic elastic parameters from morphological parameters of cancellous bone for osteoporotic lumbar spine

In biomechanics, large finite element models with macroscopic representation of several bones or joints are necessary to analyze implant failure mechanisms. In order to handle large simulation models of human bone, it is crucial to homogenize the trabecular structure regarding the mechanical behavior without losing information about the realistic material properties. Accordingly, morphology and fabric measurements of 60 vertebral cancellous bone samples from three osteoporotic lumbar spines were performed on the basis of X-ray microtomography (μCT) images to determine anisotropic elastic parameters as a function of bone density in the area of pedicle screw anchorage. The fabric tensor was mapped in cubic bone volumes by a 3D mean-intercept-length method. Fabric measurements resulted in a high degree of anisotropy (DA = 0.554). For the Young's and shear moduli as a function of bone volume fraction (BV/TV, bone volume/total volume), an individually fit function was determined and high correlations were found (97.3 ≤ R² ≤ 99.1,p < 0.005). The results suggest that the mathematical formulation for the relationship between anisotropic elastic constants and BV/TV is applicable to current μCT data of cancellous bone in the osteoporotic lumbar spine. In combination with the obtained results and findings, the developed routine allows determination of elastic constants of osteoporotic lumbar spine. Based on this, the elastic constants determined using homogenization theory can enable efficient investigation of human bone using finite element analysis (FEA). Graphical Abstract Cancellous Bone with Fabric Tensor Ellipsoid representing anisotropy and principal axis (colored coordinate system) of given trabecular structure.