The numerical simulation of biocomposites consisting of zirconia-based ceramics and cortical bone was performed with the use of a multilevel approach. The mechanical properties of the ceramic biocomposite were determined. The evolution of mesoscopic stress distributions in the biocomposite components during the process of its deformation was investigated, taking into account damage accumulation up to the fulfillment of the macro strength criterion. It is shown that damage accumulation has an impact on the stress distribution laws at the mesoscopic level, which is manifested through the appearance of a threshold for the stress distribution, as well as through a significant decrease in the distribution amplitude.