Modeling cortical bone adaptation using strain gradients

Cyclic and low-magnitude loading promotes osteogenesis (i.e. new bone formation). Normal strain, strain energy density and fatigue damage accumulation are typically considered as osteogenic stimuli in computer models to predict site-specific new bone formation. These models however had limited success in explaining osteogenesis near the sites of minimal normal strain, for example, neutral axis of bending. Other stimuli such as fluid motion or strain gradient also stimulate bone formation. In silico studies modeled the new bone formation as a function of fluid motion, however, computation of fluid motion involves complex mathematical calculations. Strain gradients drive fluid flow and thus can also be established as the stimulus. Osteogenic potential of strain gradients is however not well established. The present study establishes strain gradients as osteogenic stimuli. Bending-induced strain gradients are computed at cortical bone cross-sections reported in animal loading in vivo studies. Correlation analysis between strain gradients and site of osteogenesis is analyzed. In silico model is also developed to test the osteogenic potential of strain gradients. The model closely predicts in vivo new bone distribution as a function of strain gradients. The outcome establishes strain gradient as computationally easy and robust stimuli to predict site-specific osteogenesis. The present study may be useful in the development of biomechanical approaches to mitigate bone loss.