In-silico simulation of nanoindentation on bone using a 2D cohesive finite element model

This study proposed and validated a 2D finite element (FE) model for conducting in-silico simulations of in-situ nanoindentation tests on mineralized collagen fibrils (MCF) and the extrafibrillar matrix (EFM) within human cortical bone. Initially, a multiscale cohesive FE model was developed by adapting a previous model of bone lamellae, encompassing both MCF and EFM. Subsequently, nanoindentation tests were simulated in-silico using this model, and the resulting predictions were compared to AFM nanoindentation test data to verify the model's accuracy. The FE model accurately predicted nanoindentation results under wet conditions, closely aligning with outcomes obtained from AFM nanoindentation tests. Specifically, it successfully mirrored the traction/separation curve, nanoindentation modulus, plastic energy dissipation, and plastic energy ratio obtained from AFM nanoindentation tests. Additionally, this in-silico model demonstrated its ability to capture alterations in nanoindentation properties caused by the removal of bound water, by considering corresponding changes in mechanical properties of the collagen phase and the interfaces among bone constituents. Notably, significant changes in the elastic modulus and plastic energy dissipation were observed in both MCF and EFM compartments of bone, consistent with observations in AFM nanoindentation tests. These findings indicate that the proposed in-silico model effectively captures the influence of ultrastructural changes on bone's mechanical properties at sub-lamellar levels. Presently, no experimental methods exist to conduct parametric studies elucidating the ultrastructural origins of bone tissue fragility. The introduction of this in-silico model presents an invaluable tool to bridge this knowledge gap in the future.