A simplified homogenization model applied to viscoelastic behavior of cortical bone at ultrasonic frequencies

Cortical bone is a complex multiscale medium and its study is of importance for clinical fracture prevention. In particular, cortical attenuation is known to be linked with shock energy absorption and ability to resist fracture. However, the links between cortical bone absorption and its multiscale structure are still not well understood. This work is about the use of homogenized tensors in order to characterize the viscoelastic behavior of cortical bone at ultrasonic frequencies, i.e., about 0.1 to 10 MHz. Such tensors are derived from the cell problem via two-scale homogenization theory for linear elastic and Kelvin-Voigt viscoelastic descriptions. The elliptic formulations obtained from the cell problems are implemented within the range of medically-observed porosities. Microstructure is assessed considering cubic cells with cylindrical inclusion and transverse isotropic assumption. A simplified model, adding one temporal parameter τ per phase, allows a good agreement with experimental data. The corresponding attenuation is proportional to the square of the frequency, in agreement with Kramer-Kronig relations. This development is proposed in the context of robust clinical inverse problem approaches using a restricted number of parameter. Two main properties for the material filling the pores are adjusted and discussed: absorption and shear contribution. Best agreement with experimental data is observed for material inside the pores being solid and highly attenuating.