It is generally accepted that trabecular architecture plays a pivotal role in the mechanical behaviour of bone.With age, bone undergoes structural changes, which can result in osteoporosis, leading to lifethreatening fractures, and inevitable decrease in the quality of life.While mathematical laws governing bone remodelling are under continued investigation, the aim of this project was to apply a simple in silico model to simulate changes in the bone architecture due to age, as previously reported in clinical studies.In addition, the effects of the current recommended treatments were investigated.Using high-resolution three-dimensional µCT scans of whole human vertebrae, age-related bone loss and recovery simulation produced realistic simulations of structural change over 30 years.
Modelling and remodelling are the processes by which bone adapts its shape and internal structure to external influences. However, the cellular mechanisms triggering osteoclastic resorption and osteoblastic formation are still unknown. In order to investigate current biological theories, in silico models can be applied. In the past, most of these models were based on the continuum assumption, but some questions related to bone adaptation can be addressed better by models incorporating the trabecular microstructure. In this paper, existing simulation models are reviewed and one of the microstructural models is extended to test the hypothesis that bone adaptation can be simulated without particular knowledge of the local strain distribution in the bone. Validation using an experimental murine loading model showed that this is possible. Furthermore, the experimental model revealed that bone formation cannot be attributed only to an increase in trabecular thickness but also to structural reorganization including the growth of new trabeculae. How these new trabeculae arise is still an unresolved issue and might be better addressed by incorporating other levels of hierarchy, especially the cellular level. The cellular level sheds light on the activity and interplay between the different cell types, leading to the effective change in the whole bone. For this reason, hierarchical multi-scale simulations might help in the future to better understand the biomathematical laws behind bone adaptation.