Computational study of the mechanical influence of lacunae and perilacunar zones in cortical bone microcracking

The mechanical behavior of cortical bone is influenced by microstructural components such as osteons, Haversian canals, and osteocyte lacunae that arise from biological remodeling processes. This study takes a computational approach to investigate the role of the perilacunar zones formed by the local remodeling processes of lacunar-dwelling osteocytes by utilizing phase-field finite element models based on histological imaging of human bone. The models simulated the microdamage accumulation that occurs in cortical bone under transverse compression in bone without lacunae, with lacunae, and with a perilacunar zone surrounding lacunae in order to investigate the role of these features. The results of the simulations found that while lacunae create stress concentration which initiate further damage, perilacunar regions can delay or prevent the emergence and growth of microcracks.