Three-dimensional Simulation of Fracture Repair in the Human Tibia

Finite element models of bones can be generated based on images obtained non-invasively in the clinic. One area where such models may prove useful is in the assessment of fracture healing of long bones. To establish the feasibility of such a proposal, a three dimensional finite element model of a fractured tibia was generated, and a model of tissue differentiation and bone regeneration was used to simulate the progress of healing under two different loading magnitudes. Healing is successful under the lower load and unsuccessful under the higher load--this proves that the model has the potential to identify loads that would cause healing to fail. Following a proposal by Richardson et al. [J. Bone Jt Surg. Vol. 76B (1994) pp. 389-394] that the bending stiffness can be used to assess the extent of healing, the bending stiffness was computed during healing--it was shown that the stiffness changed in a similar manner that observed clinically. In conclusion, the paper establishes that 3D computer simulation could be a tool for assessment of the fracture healing under different orthopedic treatments.     

Rigidity and stress analyses of external fracture fixation devices—A theoretical approach

The rigidity and stresses in external fracture fixation devices were studied by means of the finite element method. Different geometries and material parameters were simulated using a beam element model. Axial, bending and torsional loads were applied through the bone ends and the displacement obtained at the fracture sites was used to calculate the fracture fixation stiffness. The key parameters which increased fixation rigidity were identified. High pin stresses were predicted under certain application conditions. Possible clinical implications for the use of such apparatus are discussed in the light of bone fracture healing. The present results are expected to have a significant impact on future design modifications and clinical applications of this popular instrument in orthopedic surgery and traumatology.     

In vivo measurement of bending stiffness in fracture healing

Background  

Measurement of the bending stiffness a healing fracture represents a valid variable in the assessment of fracture healing. However, currently available methods typically have high measurement errors, even for mild pin loosening. Furthermore, these methods cannot provide actual values of bending stiffness, which precludes comparisons among individual fractures. Thus, even today, little information is available with regards to the fracture healing pattern with respect to actual values of bending stiffness. Our goals were, therefore: to develop a measurement device that would allow accurate and sensitive measurement of bending stiffness, even in the presence of mild pin loosening; to describe the course of healing in individual fractures; and help to evaluate whether the individual pattern of bending stiffness can be predicted at an early stage of healing.     

Prediction of the Time Course of Callus Stiffness as a Function of Mechanical Parameters in Experimental Rat Fracture Healing Studies - A Numerical Study

Numerous experimental fracture healing studies are performed on rats, in which different experimental, mechanical parameters are applied, thereby prohibiting direct comparison between each other. Numerical fracture healing simulation models are able to predict courses of fracture healing and offer support for pre-planning animal experiments and for post-hoc comparison between outcomes of different in vivo studies. The aims of this study are to adapt a pre-existing fracture healing simulation algorithm for sheep and humans to the rat, to corroborate it using the data of numerous different rat experiments, and to provide healing predictions for future rat experiments. First, material properties of different tissue types involved were adjusted by comparing experimentally measured callus stiffness to respective simulated values obtained in three finite element (FE) models. This yielded values for Young’s moduli of cortical bone, woven bone, cartilage, and connective tissue of 15,750 MPa, 1,000 MPa, 5 MPa, and 1 MPa, respectively. Next, thresholds in the underlying mechanoregulatory tissue differentiation rules were calibrated by modifying model parameters so that predicted fracture callus stiffness matched experimental data from a study that used rigid and flexible fixators. This resulted in strain thresholds at higher magnitudes than in models for sheep and humans. The resulting numerical model was then used to simulate numerous fracture healing scenarios from literature, showing a considerable mismatch in only 6 of 21 cases. Based on this corroborated model, a fit curve function was derived which predicts the increase of callus stiffness dependent on bodyweight, fixation stiffness, and fracture gap size. By mathematically predicting the time course of the healing process prior to the animal studies, the data presented in this work provides support for planning new fracture healing experiments in rats. Furthermore, it allows one to transfer and compare new in vivo findings to previously performed studies with differing mechanical parameters.     

Computational simulation of bone fracture healing under inverse dynamisation

Adaptive finite element models have allowed researchers to test hypothetical relationships between the local mechanical environment and the healing of bone fractures. However, their predictive power has not yet been demonstrated by testing hypotheses ahead of experimental testing. In this study, an established mechano-biological scheme was used in an iterative finite element simulation of sheep tibial osteotomy healing under a hypothetical fixation regime, “inverse dynamisation”. Tissue distributions, interfragmentary movement and stiffness across the fracture site were compared between stiff and flexible fixation conditions and scenarios in which fixation stiffness was increased at a discrete time-point. The modelling work was conducted blind to the experimental study to be published subsequently. The simulations predicted the fastest and most direct healing under constant stiff fixation, and the slowest healing under flexible fixation. Although low fixation stiffness promoted more callus formation prior to bridging, this conferred little additional stiffness to the fracture in the first 5 weeks. Thus, while switching to stiffer fixation facilitated rapid subsequent bridging of the fracture, no advantage of inverse dynamisation could be demonstrated. In vivo data remains necessary to conclusively test this treatment protocol and this will, in turn, provide an evaluation of the model’s performance. The publication of both hypotheses and their computational simulation, prior to experimental testing, offers an appealing means to test the predictive power of mechano-biological models.     

Determination of mechanical stiffness of bone by pQCT measurements: correlation with non-destructive mechanical four-point bending test data

Mechanical tests of bone provide valuable information about material and structural properties important for understanding bone pathology in both clinical and research settings, but no previous studies have produced applicable non-invasive, quantitative estimates of bending stiffness. The goal of this study was to evaluate the effectiveness of using peripheral quantitative computed tomography (pQCT) data to accurately compute the bending stiffness of bone. Normal rabbit humeri (N=8) were scanned at their mid-diaphyses using pQCT. The average bone mineral densities and the cross-sectional moments of inertia were computed from the pQCT cross-sections. Bending stiffness was determined as a function of the elastic modulus of compact bone (based on the local bone mineral density), cross-sectional moment of inertia, and simulated quasistatic strain rate. The actual bending stiffness of the bones was determined using four-point bending tests. Comparison of the bending stiffness estimated from the pQCT data and the mechanical bending stiffness revealed excellent correlation (R2=0.96). The bending stiffness from the pQCT data was on average 103% of that obtained from the four-point bending tests. The results indicate that pQCT data can be used to accurately determine the bending stiffness of normal bone. Possible applications include temporal quantification of fracture healing and risk management of osteoporosis or other bone pathologies.     

A new device to quantify regenerate torsional stiffness in distraction osteogenesis.

We present a newly developed torsional stiffness measurement device with the potential to quantitatively assess the in vivo torsional stiffness of bone regenerate during distraction osteogenesis. We describe the form and function of this device and its application in a model of regenerate consolidation. The device was able to produce data to assess stiffness of the regenerate with an accuracy between +/- 3 and +/- 9% for material stiffness ranging between 0.1 and 2.4 Nm/o and with a precision of +/- 3.6%. This method provides advantages over similar methods of bone fracture healing assessment with guaranteed maintenance of bone axis, minimized risk of bone misalignment during the bone healing process and a close relation to the functional loading pattern in torsion of bones such as tibia and femora.     

Virtual structural analysis of tibial fracture healing from low-dose clinical CT scans

Quantitative assessment of bone fracture healing remains a significant challenge in orthopaedic trauma research. Accordingly, we developed a new technique for assessing bone healing using virtual mechano-structural analysis of computed tomography (CT) scans. CT scans from 19 fractured human tibiae at 12 weeks after surgery were segmented and prepared for finite element analysis (FEA). Boundary conditions were applied to the models to simulate a torsion test that is commonly used to access the structural integrity of long bones in animal models of fracture healing. The output of each model was the virtual torsional rigidity (VTR) of the healing zone, normalized to the torsional rigidity of each patient’s virtually reconstructed tibia. This provided a structural measure to track the percentage of healing each patient had undergone. Callus morphometric measurements were also collected from the CT scans. Results showed that at 12 weeks post-op, more than 75% of patients achieved a normalized VTR (torsional rigidity relative to uninjured bone) of 85% or above. The predicted intact torsional rigidities compared well with published cadaveric data. Across all patients, callus volume and density were weakly and non-significantly correlated with normalized VTR and time to clinical union. Conversely, normalized VTR was significantly correlated with time to union (R² = 0.383, p = 0.005). This suggests that fracture scoring methods based on the visual appearance of callus may not accurately predict mechanical integrity. The image-based structural analysis presented here may be a useful technique for assessment of bone healing in orthopaedic trauma research.     

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.     

A 3D computational simulation of fracture callus formation: influence of the stiffness of the external fixator

The stiffness of the external fixation highly influences the fracture healing pattern. In this work we study this aspect by means of a finite element model of a simple transverse mid-diaphyseal fracture of an ovine metatarsus fixed with a bilateral external fixator. In order to simulate the regenerative process, a previously developed mechanobiological model of bone fracture healing was implemented in three dimensions. This model is able to simulate tissue differentiation, bone regeneration, and callus growth. A physiological load of 500 N was applied and three different stiffnesses of the external fixator were simulated (2300, 1725, and 1150 N/mm). The interfragmentary strain and load sharing mechanism between bone and the external fixator were compared to those recorded in previous experimental works. The effects of the stiffness on the callus shape and tissue distributions in the fracture site were also analyzed. We predicted that a lower stiffness of the fixator delays fracture healing and causes a larger callus, in correspondence to well-documented clinical observations.     

Quantitative computed tomography-based finite element models of the human lumbar vertebral body: effect of element size on stiffness,damage, and fracture strength predictions

This study investigated the numerical convergence characteristics of specimen-specific "voxel-based" finite element models of 14 excised human cadaveric lumbar vertebral bodies (age: 37-87; M = 6, F = 8) that were generated automatically from clinical-type CT scans. With eventual clinical applications in mind, the ability of the model stiffness to predict the experimentally measured compressive fracture strength of the vertebral bodies was also assessed. The stiffness of "low"-resolution models (3 x 3 x 3 mm element size) was on average only 4% greater (p = 0.03) than for "high"-resolution models (1 x 1 x 1.5 mm) despite interspecimen variations that varied over four-fold. Damage predictions using low- vs high-resolution models were significantly different (p = 0.01) at loads corresponding to an overall strain of 0.5%. Both the high (r2 = 0.94) and low (r2 = 0.92) resolution model stiffness values were highly correlated with the experimentally measured ultimate strength values. Because vertebral stiffness variations in the population are much greater than those that arise from differences in voxel size, these results indicate that imaging resolution is not critical in cross-sectional studies of this parameter. However, longitudinal studies that seek to track more subtle changes in stiffness over time should account for the small but highly significant effects of voxel size. These results also demonstrate that an automated voxel-based finite element modeling technique may provide an excellent noninvasive assessment of vertebral strength.     

Quantitative measures of femoral fracture repair in rats derived by micro-computed tomography

Although fracture healing is frequently studied in pre-clinical models of long bone fractures using rodents, there is a dearth of objective quantitative techniques to assess successful healing. Biomechanical testing is possibly the most quantitative and relevant to a successful clinical outcome, but it is a destructive technique providing little insight into the cellular mechanisms associated with healing. The advent of X-ray computed tomography (CT) has provided the opportunity to quantitatively and non-destructively assess bone structure and density, but it is unknown how measurements derived using this technology relate to successful healing. To examine possible relationships, we used a pre-clinical model to test for statistically significant correlations between quantitative characteristics of the callus by micro-CT (μCT) and the bending strength, stiffness, and energy-to-failure of the callus as assessed by three-point bending of excised bones. A closed, transverse fracture was generated in the mid-shaft of rat femurs by impact loading. Shortly thereafter, the rats received a one-time, local injection of either the vehicle or one of four doses of lovastatin. Following sacrifice after 4 weeks of healing, fractured femurs were extracted for μCT analysis and then three-point bending. Setting the region of interest to be 3.2 mm above and below the fracture line, we acquired standard and new μCT-derived measurements. The mineralized callus volume and the mineral density of the callus correlated positively with callus strength (r_xy=?0.315, p=0.016 and r_xy=0.444, p<0.0005, respectively) and stiffness (r_xy=?0.271, p=0.040 and r_xy=0.325, p=0.013, respectively), but the fraction of the callus that mineralized and the moment of inertia of the callus did not. This fraction did correlate with energy-to-failure (r_xy=?0.343, p=0.0085). Of the μCT-derived measurements, quantifying defects within the outer bridging cortices of the callus produced the strongest correlation with both callus strength (r_xy=0.557, p<0.0001) and stiffness (r_xy=0.468, p=0.0002). By both reducing structural defects and increasing mineralization, lovastatin appears to increase the callus strength.     

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.     

A finite element inverse analysis to assess functional improvement during the fracture healing process

Assessment of the restoration of load-bearing function is the central goal in the study of fracture healing process. During the fracture healing, two critical aspects affect its analysis: (1) material properties of the callus components, and (2) the spatio-temporal architecture of the callus with respect to cartilage and new bone formation. In this study, an inverse problem methodology is used which takes into account both features and yields material property estimates that can analyze the healing changes. Six stabilized fractured mouse tibias are obtained at two time points during the most active phase of the healing process, respectively 10 days (n=3), and 14 days (n=3) after fracture. Under the same displacement conditions, the inverse procedure estimations of the callus material properties are generated and compared to other fracture healing metrics. The FEA estimated property is the only metric shown to be statistically significant (p=0.0194) in detecting the changes in the stiffness that occur during the healing time points. In addition, simulation studies regarding sensitivity to initial guess and noise are presented; as well as the influence of callus architecture on the FEA estimated material property metric. The finite element model inverse analysis developed can be used to determine the effects of genetics or therapeutic manipulations on fracture healing in rodents.     

Localized damage in vertebral bone is most detrimental in regions of high strain energy density

It was hypothesized that damage to bone tissue would be most detrimental to the structural integrity of the vertebral body if it occurred in regions with high strain energy density, and not necessarily in regions of high or low trabecular bone apparent density, or in a particular anatomic location. The reduction in stiffness due to localized damage was computed in 16 finite element models of 10-mm-thick human vertebral sections. Statistical analyses were performed to determine which characteristic at the damage location--strain energy density, apparent density, or anatomic location--best predicted the corresponding stiffness reduction. There was a strong positive correlation between regional strain energy density and structural stiffness reduction in all 16 vertebral sections for damage in the trabecular centrum (p < 0.05, r2 = 0.43-0.93). By contrast, regional apparent density showed a significant negative correlation to stiffness reduction in only four of the sixteen bones (p < 0.05, r2 = 0.47-0.58). While damage in different anatomic locations did lead to different reductions in stiffness (p < 0.0001, ANOVA), no single location was consistently the most critical location for damage. Thus, knowledge of the characteristics of bone that determine strain energy density distributions can provide an understanding of how damage reduces whole bone mechanical properties. A patient-specific finite element model displaying a map of strain energy density can help optimize surgical planning and reinforcement of bone in individuals with high fracture risk.     

Finite element analysis of acetabular fractures—development and validation with a synthetic pelvis

Acetabular fracture presents a challenging situation to trauma surgeons today due to its complexity. Finite element (FE) models can be of great help as they can improve the surgical planning and post surgery patient management for those with acetabular fractures. We have developed a non-linear finite element model of the pelvis and validated its fracture prediction capability with synthetic polyurethane pelves. A mechanical experiment was performed with the synthetic bones and fracture loads and patterns were observed for two different loading cases. Fracture loads predicted by our FE model were within one standard deviation of the experimental fracture loads for both loading cases. The incipient fracture pattern predicted by the model also resembled the actual pattern from the experiment. Although it is not a complete validation with human cadaver bones, the good agreement between model predictions and experimental results indicate the validity of our approach in using non-linear FE formulation along with contact conditions in predicting bone fractures.     

Subject-specific finite element models implementing a maximum principal strain criterion are able to estimate failure risk and fracture location on human femurs tested in vitro

No agreement on the choice of the failure criterion to adopt for the bone tissue can be found in the literature among the finite element studies aiming at predicting fracture risk of bones. The use of stress-based criteria seems to prevail on strain-based ones, while basic bone biomechanics suggest using strain parameters to describe failure. The aim of the present combined experimental-numerical study was to verify, using subject-specific finite element models able to accurately predict strains, if a strain-based failure criterion could identify the failure patterns of bones. Three cadaver femurs were CT-scanned and subsequently fractured in a clinically relevant single-stance loading scenario. Load-displacement curves and high-speed movies were acquired to define the failure load and the location of fracture onset, respectively. Subject-specific finite element models of the three femurs were built from CT data following a validated procedure. A maximum principal strain criterion was implemented in the finite element models, and two stress-based criteria selected for comparison. The failure loads measured were applied to the models, and the computed risks of fracture were compared to the results of the experimental tests. The proposed principal strain criterion managed to correctly identify the level of failure risk and the location of fracture onset in all the modelled specimens, while Von Mises or maximum principal stress criteria did not give significant information. A maximum principal strain criterion can thus be defined a suitable candidate for the in vivo risk factor assessment on long bones.     

Effects of geometric individualisation of a human spine model on load sharing: neuro-musculoskeletal simulation reveals significant differences in ligament and muscle contribution

In spine research, two possibilities to generate models exist: generic (population-based) models representing the average human and subject-specific representations of individuals. Despite the increasing interest in subject specificity, individualisation of spine models remains challenging. Neuro-musculoskeletal (NMS) models enable the analysis and prediction of dynamic motions by incorporating active muscles attaching to bones that are connected using articulating joints under the assumption of rigid body dynamics. In this study, we used forward-dynamic simulations to compare a generic NMS multibody model of the thoracolumbar spine including fully articulated vertebrae, detailed musculature, passive ligaments and linear intervertebral disc (IVD) models with an individualised model to assess the contribution of individual biological structures. Individualisation was achieved by integrating skeletal geometry from computed tomography and custom-selected muscle and ligament paths. Both models underwent a gravitational settling process and a forward flexion-to-extension movement. The model-specific load distribution in an equilibrated upright position and local stiffness in the L4/5 functional spinal unit (FSU) is compared. Load sharing between occurring internal forces generated by individual biological structures and their contribution to the FSU stiffness was computed. The main finding of our simulations is an apparent shift in load sharing with individualisation from an equally distributed element contribution of IVD, ligaments and muscles in the generic spine model to a predominant muscle contribution in the individualised model depending on the analysed spine level.

     

Mineral heterogeneity affects predictions of intratrabecular stress and strain

Knowledge of the influence of mineral variations (i.e., mineral heterogeneity) on biomechanical bone behavior at the trabecular level is limited. The aim of this study is to investigate how this material property affects the intratrabecular distributions of stress and strain in human adult trabecular bone. Two different sets of finite element (FE) models of trabecular samples were constructed; tissue stiffness was either scaled to the local degree of mineralization of bone as measured with microCT (heterogeneous) or tissue stiffness was assumed to be homogeneous. The influence of intratrabecular mineral heterogeneity was analyzed by comparing both models. Interesting effects were seen regarding intratrabecular stress and strain distributions. In the homogeneous model, the highest stresses were found at the surface with a significant decrease towards the core. Higher superficial stresses could indicate a higher predicted fracture risk in the trabeculae. In the heterogeneous model this pattern was different. A significant increase in stress with increasing distance from the trabecular surface was found followed by a significant decrease towards the core. This suggests trabecular bending during a compression. In both models a decrease in strain values from surface to core was predicted, which is consistent with trabecular bending. When mineral heterogeneity was taken into account, the predicted intratrabecular patterns of stress and strain are more consistent with the expected biomechanical behavior as based on mineral variations in trabeculae. Our findings indicate that mineral heterogeneity should not be neglected when performing biomechanical studies on topics such as the (long-term or dose dependent) effects of antiresorptive treatments.     

Characterizing gait induced normal strains in a murine tibia cortical bone defect model

The critical role that mechanical stimuli serve in mediating bone repair is recognized but incompletely understood. Further, previous attempts to understand this role have utilized application of externally applied mechanical loads to study the tissue’s response. In this project, we have therefore endeavored to capitalize on bone’s own consistently diverse loading environment to develop a novel model that would enable assessment of the influence of physiologically engendered mechanical stimuli on cortical defect repair. We used an inverse dynamics approach with finite element analysis (FEA) to first quantify normal strain distributions generated in mouse tibia during locomotion. The strain environment of the tibia, as previously reported for other long bones, was found to arise primarily due to bending and was consistent in orientation through the stance phase of gait. Based on these data, we identified three regions within a transverse cross-section of the mid-diaphysis as uniform locations of either peak tension, peak compression, or the neutral axis of bending (i.e. minimal strain magnitude). We then used FEA to quantify the altered strain environment that would be produced by a 0.6 mm diameter cylindrical cortical bone defect at each diaphyseal site and, in an in situ study confirmed our ability to accurately place defects at the desired diaphyseal locations. The resulting model will enable the exploration of cortical bone healing within the context of physiologically engendered mechanical strain.