Numerical simulation of callus healing for optimization of fracture fixation stiffness

The stiffness of fracture fixation devices together with musculoskeletal loading defines the mechanical environment within a long bone fracture, and can be quantified by the interfragmentary movement. In vivo results suggested that this can have acceleratory or inhibitory influences, depending on direction and magnitude of motion, indicating that some complications in fracture treatment could be avoided by optimizing the fixation stiffness. However, general statements are difficult to make due to the limited number of experimental findings. The aim of this study was therefore to numerically investigate healing outcomes under various combinations of shear and axial fixation stiffness, and to detect the optimal configuration. A calibrated and established numerical model was used to predict fracture healing for numerous combinations of axial and shear fixation stiffness under physiological, superimposed, axial compressive and translational shear loading in sheep. Characteristic maps of healing outcome versus fixation stiffness (axial and shear) were created. The results suggest that delayed healing of 3 mm transversal fracture gaps will occur for highly flexible or very rigid axial fixation, which was corroborated by in vivo findings. The optimal fixation stiffness for ovine long bone fractures was predicted to be 1000-2500 N/mm in the axial and >300 N/mm in the shear direction. In summary, an optimized, moderate axial stiffness together with certain shear stiffness enhances fracture healing processes. The negative influence of one improper stiffness can be compensated by adjustment of the stiffness in the other direction.