A microstructural finite element simulation of mechanically induced bone formation

A finite element method to simulate the formation of an interconnected trabectular bone microstructure oriented with respect to applied in vivo mechanical forces is introduced and quantitatively compared to experimental data from a hydraulic bone chamber implant model. Randomly located 45 microm mineralized nodules were used as the initial condition for the model simulations to represent an early stage of intramembranous bone formation. Boundary conditions were applied consistent with the mechanical environment provided by the in vivo bone chamber model. A two-dimensional repair simulation algorithim that incorporated strain energy density (SED), SED gradient, principal strain, or principal strain gradient as the local objective criterion was utilized to simulate the formation of an oriented trabecular bone microstructure. The simulation solutions were convergent, unique, and relatively insensitive to the assumed initial distribution of mineralized nodules. Model predictions of trabecular bone morphology and anisotropy were quantitatively compared to experimental results. All simulations produced structures that qualitatively resembled oriented trabecular bone. However only simulations utilizing a gradient objective criterion yielded results quantitatively similar to in vivo observations. This simulation approach coupled with an experimental model that delivers controlled in vivo mechanical stimuli can be utilized to study the relationship between physical factors and microstructural adaptation during bone repair.