Influence of fracture gap size on the pattern of long bone healing: a computational study

Following fractures, bones restore their original structural integrity through a complex process in which several cellular events are involved. Among other factors, this process is highly influenced by the mechanical environment of the fracture site. In this study, we present a mathematical model to simulate the effect of mechanical stimuli on most of the cellular processes that occur during fracture healing, namely proliferation, migration and differentiation. On the basis of these three processes, the model then simulates the evolution of geometry, distributions of cell types and elastic properties inside a healing fracture. The three processes were implemented in a Finite Element code as a combination of three coupled analysis stages: a biphasic, a diffusion and a thermoelastic step. We tested the mechano-biological regulatory model thus created by simulating the healing patterns of fractures with different gap sizes and different mechanical stimuli. The callus geometry, tissue differentiation patterns and fracture stiffness predicted by the model were similar to experimental observations for every analysed situation.