Abstract
Structural integrity of bone tissue plays an important role in daily activities of humans. However, traumatic incidents such as sports injuries, collisions and falls can cause bone fracture, servere pain and mobility loss. In addition, ageing and degenerative bone diseases such as osteoporosis can increase the risk of fracture [1]. As a composite-like material, a cortical bone tissue is capable of tolerating moderate fracture/cracks without complete failure. The key to this is its heterogeneously distributed microstructural constituents providing both intrinsic and extrinsic toughening mechanisms. At micro-scale level, cortical bone can be considered as a four-phase composite material consisting of osteons, Haversian canals, cement lines and interstitial matrix. These microstructural constituents can directly affect local distributions of stresses and strains, and, hence, crack initiation and propagation. Therefore, understanding the effect of micromorphology of cortical bone on crack initiation and propagation, especially under dynamic loading regimes is of great importance for fracture risk evaluation. In this study, random microstructures of a cortical bone tissue were modelled with finite elements for four groups: healthy (control), young age, osteoporosis and bisphosphonate-treated, based on osteonal morphometric parameters measured from microscopic images for these groups. The developed models were loaded under the same dynamic loading conditions, representing a direct impact incident, resulting in progressive crack propagation. An extended finite-element method (X-FEM) was implemented to realize solution-dependent crack propagation within the microstructured cortical bone tissues. The obtained simulation results demonstrate significant differences due to micromorphology of cortical bone, in terms of crack propagation characteristics for different groups, with the young group showing highest fracture resistance and the senior group the lowest.