Do microcracks decrease or increase fatigue resistance in cortical bone?

Fatigue of cortical bone produces microcracks; it has been hypothesized that these cracks are analogous to those occurring in engineered composite materials and constitute a similar mechanism for fatigue resistance. However, the numbers of these linear microcracks increase substantially with age, suggesting that they contribute to increased fracture incidence among the elderly. To test these opposing hypotheses, we fatigued 20 beams of femoral cortical bone from elderly men and women in load-controlled four point bending having initial strain ranges of 3000 or 5000 microstrain. Loading was stopped at fracture or 10(6) cycles, whichever occurred first, and microcrack density and length were measured in the loaded region and in a control region that was not loaded. We studied the dependence of fatigue life and induced microdamage on initial microdamage, cortical region, subject gender and age, and several other variables. When the effect of modulus variability was controlled, longer fatigue life was associated with higher rather than lower initial crack density, particularly in the medial cortex. The increase in crack density following fatigue loading was greater in specimens from older individuals and those initially having longer microcracks. Crack density increased as much in specimens fatigued short of the failure point as in those that fractured, and microcracks were, on average, shorter in specimens with greater numbers of resorption spaces, a measure of remodeling rate.     

Non Destructive Characterization of Cortical Bone Micro-Damage by Nonlinear Resonant Ultrasound Spectroscopy

The objective of the study was to evaluate the ability of a nonlinear ultrasound technique, the so-called nonlinear resonant ultrasound spectroscopy (NRUS) technique, for detecting early microdamage accumulation in cortical bone induced by four-point bending fatigue. Small parallelepiped beam-shaped human cortical bone specimens were subjected to cyclic four-point bending fatigue in several steps. The specimens were prepared to control damage localization during four-point bending fatigue cycling and to unambiguously identify resonant modes for NRUS measurements. NRUS measurements were achieved to follow the evolution of the nonlinear hysteretic elastic behavior during fatigue-induced damage. After each fatigue step, a small number of specimens was removed from the protocol and set apart to quantitatively assess the microcrack number density and length using synchrotron radiation micro-computed tomography (SR-µCT). The results showed a significant effect of damage steps on the nonlinear hysteretic elastic behavior. No significant change in the overall length of microcracks was observed in damaged regions compared to the load-free control regions. Only an increased number of shortest microcracks, those in the lowest quartile, was noticed. This was suggestive of newly formed microcracks during the early phases of damage accumulation. The variation of nonlinear hysteretic elastic behavior was significantly correlated to the variation of the density of short microcracks. Our results suggest that the nonlinear hysteretic elastic behavior is sensitive to early bone microdamage. Therefore NRUS technique can be used to monitor fatigue microdamage progression in in vitro experiments.     

The behaviour of microcracks in compact bone

This paper summarises four separate studies carried out by our group over the past number of years in the area of bone microdamage. The first study investigated the manner by which microcracks accumulate and interact with bone microstructure during fatigue testing of compact bone specimens. In a series of fatigue tests carried out at four different stress ranges between 50 and 80 MPA, crack density increased with loading cycles at a rate determined by the applied stress. Variations in the patterns of microdamage accumulation suggest that that at low stress levels, larger amounts of damage can build up without failure occurring. In a second study using a series of four-pont bending tests carried out on ovine bone samples, it was shown that bone microstructure influenced the ability of microcracks to propagate, with secondary osteons acting as barriers to crack growth. In a third study, the manner by which crack growth disrupts the canalicular processes connecting osteocytes was investigated. Analysis of individual cracks showed that disruption of the canalicular processes connecting osteocytes occurred due to shear displacement at the face of propagating microcracks, suggesting that this may play some role in the mechanism that signals bone remodelling. In a fourth in vivo study, it was shown that altering the mechanical load applied to the long bones of growing rats causes microcrack formation. In vivo microdamage was present in rats subjected to hindlimb suspension with a higher microcrack density found in the humeri than the femora. Microdamage was also found in control animals. This is the first study to demonstrate in vivo microcracks in normally loaded bones in a rat model.     

A fatigue microcrack alters fluid velocities in a computational model of interstitial fluid flow in cortical bone

Targeted remodeling is activated by fatigue microcracks and plays an important role in maintaining bone integrity. It is widely believed that fluid flow-induced shear stress plays a major role in modulating the mechanotransduction process. Therefore, it is likely that fluid flow-induced shear stress plays a major role in the initiation of the repair of fatigue damage. Since no in vivo measurements of fluid flow within bone exist, computational and mathematical models must be employed to investigate the fluid flow field and the shear stress occurring within cortical bone. We developed a computational fluid dynamic model of cortical bone to examine the effect of a fatigue microcrack on the fluid flow field. Our results indicate that there are alterations in the fluid flow field as far as 150 microm away from the crack, and that at distances farther than this, the fluid flow field is similar to the fluid flow field of intact bone. Through the crack and immediately above and below it, the fluid velocity is higher, while at the lateral edges it is lower than that calculated for the intact model, with a maximum change of 29%. Our results suggest that the presence of a fatigue microcrack can alter the shear stress in regions near the crack. These alterations in shear stress have the potential to significantly alter mechanotransduction and may play a role in the initiation of the repair of fatigue microcracks.     

Deformation behaviour and damage accumulation of cortical bone specimens from the equine tibia under cyclic loading

Despite its clinical importance, the fatigue behaviour of cortical bone has not been examined as widely as its static behaviour. In the present study, specimens from the tibiae of horses have been subjected to load-controlled single step tests. The cyclic deformation behaviour was described by the development of stress-strain hysteresis parameters over the lifetime. The fatigue behaviour of bone is characterised by cyclic softening which is most distinctive towards the end of the lifetime. The microstructural damage accumulated during cyclic loading results in a loss of stiffness, asymmetrical deformation of the bone in tension and compression in cyclic creep. As shown by light and scanning electron microscopy, microcrack formation and growth is the main damage mechanism. The crack growth behaviour is strongly influenced by the microstructure, the stress components and the absolute value of the local stresses. Lower local stresses and/or compressive mean stresses lead to a dominant influence of the shear stress components with shear failure at inner interfaces. With increasing crack length, that is, higher local stress amplitudes, or tensile mean stresses, the microstructure is more and more ignored and failure occurs primarily under the influence of the normal stress components. This can be clearly seen on the fracture and specimen surfaces.     

Can deterministic mechanical size effects contribute to fracture and microdamage accumulation in trabecular bone?

Failure of bone under monotonic and cyclic loading is related to the bone mineral density, the quality of the bone matrix, and the evolution of microcracks. The theory of linear elastic fracture mechanics has commonly been applied to describe fracture in bone. Evidence is presented that bone failure can be described through a non-linear theory of fracture. Thereby, deterministic size effects are introduced. Concepts of a non-linear theory are applied to discern how the interaction among bone matrix constituents (collagen and mineral), microcrack characteristics, and trabecular architecture can create distinctively differences in the fracture resistance at the bone tissue level. The non-linear model is applied to interpret pre-clinical data concerning the effects of anti-osteoporotic agents on bone properties. The results show that bisphosphonate (BP) treatments that suppress bone remodeling will change trabecular bone in ways such that the size of the failure process zone relative to the trabecular thickness is reduced. Selective estrogen receptor modulators (SERMs) that suppress bone remodeling will change trabecular bone in ways such that the size of the failure process zone relative to the trabecular thickness is increased. The consequences of these changes are reflected in bone mechanical response and predictions are consistent with experimental observations in the animal model which show that BP treatment is associated with more brittle fracture and microcracks without altering the average length of the cracks, whereas SERM treatments lead to a more ductile fracture and mainly increase crack length with a smaller increase in microcrack density. The model suggests that BPs may be more effective in cases in which bone mass is very low, whereas SERMS may be more effective when milder osteoporotic symptoms are present.     

Experimental validation of a microcracking-based toughening mechanism for cortical bone

It has been proposed that cortical bone derives its toughness by forming microcracks during the process of crack propagation (J. Biomech. 30 (1997) 763; J. Biomech. 33 (2000) 1169). The purpose of this study was to experimentally validate the previously proposed microcrack-based toughening mechanism in cortical bone. Crack initiation and propagation tests were conducted on cortical bone compact tension specimens obtained from the antlers of red deer. For these tests, the main fracture crack was either propagated to a predetermined crack length or was stopped immediately after initiating from the notch. The microcracks produced in both groups of specimens were counted in the same surface area of interest around and below the notch, and crack growth resistance and crack propagation velocity were analyzed. There were more microcracks in the surface area of interest in the propagation than in initiation specimens showing that the formation of microcracks continued after the initiation of a fracture crack. Crack growth resistance increased with crack extension, and crack propagation velocity vs. crack extension curves demonstrated the characteristic jump increase and decrease pattern associated with the formation of microcracks. The scanning electron micrographs of crack initiation and propagation displayed the formation of a frontal process zone and a wake, respectively. These results support the microcrack-based toughening mechanism in cortical bone. Bone toughness is, therefore, determined by its ability to form microcracks during fracture.     

Propagation of surface fatigue cracks in human cortical bone

An understanding of how fatigue cracks grow in bone is of importance as fatigue is thought to be the main cause of clinical stress fractures. This study presents new results on the fatigue-crack growth behavior of small surface cracks (approximately 75-1000 microm in size) in human cortical bone, and compares their growth rates with data from other published studies on the behavior of both surface cracks and many millimeter, through-thickness large cracks. Results are obtained with a cyclically loaded cantilever-beam geometry using optical microscopy to examine for crack growth after every 100-500 cycles. Based on the current and previous results, small fatigue cracks appear to become more resistant to fatigue-crack growth with crack extension, analogous to the way the fracture resistance of cortical bone increases with crack growth. Mechanistically, a theory attributing such behavior to the development of bridges in the wake of the crack with crack growth is presented. The existence of such bridges is directly confirmed using optical microscopy.     

The role of the lamellar interface during torsional yielding of human cortical bone

Fragility fractures are a result of alterations in bone quantity, tissue properties, applied loads, or a combination of these factors. The current study addresses the contribution of cortical bone tissue properties to skeletal fragility by characterizing the shear damage accumulation processes which occur during torsional yielding in normal bone. Samples of human femoral cortical bone were loaded in torsion and damaged at a post-yield twist level. The number of microcracks within osteons, interstitial tissue, and along cement lines were assessed using basic fuchsin staining. Damage density measures (number of cracks/mm2) were correlated with stiffness degradation and changes in relaxation. Damaged samples exhibited a wide variation in total microcrack density, ranging from 1.1 to 43.3 cracks/mm2 with a mean density of 19.7 +/- 9.8 cracks/mm2. Lamellar interface cracks comprised more than 75% of the total damage, indicating that the lamellar interface is weak in shear and is a principal site of shear damage accumulation. Damage density was positively correlated with secant stiffness degradation, but only explained 22% of the variability in degradation. In contrast, damage density was uncorrelated with the changes in relaxation, indicating that a simple crack counting measure such as microcrack density was not an appropriate measure of relaxation degradation. Finally, a nonuniform microcrack density distribution was observed, suggesting that internal shear stresses were redistributed within the torsion samples during post-yield loading. The results suggested that the lamellar interface in human cortical bone plays an important role in torsional yielding by keeping cracks physically isolated from each other and delaying microcrack coalescence in order to postpone the inevitable formation of the fatal crack.     

Rupture of osteocyte processes across microcracks: the effect of crack length and stress

Bone cells are connected to one another in a network, via their dendritic cellular processes. Previously, we hypothesized that these processes could be ruptured by microcracks. We proposed this as a mechanism by which osteoctyes could detect the presence of microcracks. In order for this mechanism to be effective, the number of ruptured processes would have to increase with microcrack length and also with the applied cyclic stress. This paper presents for the first time experimental data, which shows that this is indeed the case. We examined samples of bovine, ovine and murine bone ex vivo and observed processes passing across crack faces: some were still intact whilst others had ruptured. The number of intact processes per unit crack length decreased significantly with increasing crack length and also decreased in samples, which had been tested in vitro at higher stress levels. A theoretical model that we had developed previously was able to predict the overall magnitude and general trends in the experimental data. This work has provided further support for our "scissors" model, which proposes that microcracks can be detected because they disturb the osteocyte network, specifically by rupturing cellular processes where they pass across the crack faces.     

Material heterogeneity,microstructure, and microcracks demonstrate differential influence on crack initiation and propagation in cortical bone

The recent studies have shown that long-term bisphosphonate use may result in a number of mechanical alterations in the bone tissue including a reduction in compositional heterogeneity and an increase in microcrack density. There are limited number of experimental and computational studies in the literature that evaluated how these modifications affect crack initiation and propagation in cortical bone. Therefore, in this study, the entire crack growth process including initiation and propagation was simulated at the microscale by using the cohesive extended finite element method. Models with homogeneous and heterogeneous material properties (represented at the microscale capturing the variability in material property values and their distribution) as well as different microcrack density and microstructure were compared. The results showed that initiation fracture resistance was higher in models with homogeneous material properties compared to heterogeneous ones, whereas an opposite trend was observed in propagation fracture resistance. The increase in material heterogeneity level up to 10 different material property sets increased the propagation fracture resistance beyond which a decrease was observed while still remaining higher than the homogeneous material distribution. The simulation results also showed that the total osteonal area influenced crack propagation and the local osteonal area near the initial crack affected the crack initiation behavior. In addition, the initiation fracture resistance was higher in models representing bisphosphonate treated bone (low material heterogeneity, high microcrack density) compared to untreated bone models (high material heterogeneity, low microcrack density), whereas an opposite trend was observed at later stages of crack growth. In summary, the results demonstrated that tissue material heterogeneity, microstructure, and microcrack density influenced crack initiation and propagation differently. The findings also elucidate how possible modifications in material heterogeneity and microcrack density due to bisphosphonate treatment may influence the initiation and propagation fracture resistance of cortical bone.     

Fatigue microcracks that initiate fracture are located near elevated intracortical porosity but not elevated mineralization

In vivo microcracks in cortical bone are typically observed within more highly mineralized interstitial tissue, but postmortem investigations are inherently limited to cracks that did not lead to fracture which may be misleading with respect to understanding fracture mechanisms. We hypothesized that the one fatigue microcrack which initiates fracture is located spatially adjacent to elevated intracortical porosity but not elevated mineralization. Therefore, the spatial correlation between intracortical porosity, elevated mineralization, and fatigue microdamage was investigated by combining, for the first time, sequential, nondestructive, three-dimensional micro-computed tomography (micro-CT) measurements of each in cortical bone specimens subjected to compressive fatigue loading followed by a tensile overload to fracture. Fatigue loading resulted in significant microdamage accumulation and compromised mechanical properties upon tensile overload compared to control specimens. The microdamage that initiated fracture upon tensile overload was able to be identified in all fatigue-loaded specimens using contrast-enhanced micro-CT and registered images. Two-point (or pair) correlation functions revealed a spatial correlation between microdamage at the fracture initiation site and intracortical porosity, but not highly mineralized tissue, confirming the hypothesis. This difference was unique to the fracture initiation site. Intracortical porosity and highly mineralized tissue exhibited a significantly lower and higher probability, respectively, of being located spatially adjacent to all sites of microdamage compared to the fracture initiation site. Therefore, the results of this study suggest that human cortical bone is tolerant of most microcracks, which are generally compartmentalized within the more highly mineralized interstitial tissue, but a single microcrack of sufficient size located in spatial proximity to intracortical porosity can compromise fracture resistance.     

Contribution, development and morphology of microcracking in cortical bone during crack propagation

A fracture mechanics study of cortical bone is presented to investigate the contribution, development morphology of microcracking in cortical bone during crack propagation. Post-hoc analyses of microcrack orientation, crack propagation velocity and fracture surface roughness were conducted on previously tested human and bovine bone compact tension specimens. It was found that, consistent with its higher toughness, bovine bone formed significantly more longitudinal, transverse and inclined microcracks than human bone. However, in human bone more of the microcracks that formed were longitudinal than transverse or inclined, a feature that would optimise bone's toughness. Crack propagation velocity in human and bovine bone displayed the same characteristic pattern with crack extension, where an increase in velocity is followed by a consequent decrease and vice versa. On the basis of this pattern, a model or crack propagation has been proposed. It provides a detailed account of mocrocrack formation and contribution towards the propagation of a fracture crack. Analyses of fracture surfaces indicated that, consistent with its higher toughness, bovine bone displays a rougher surface than human bone but they both have the same basic fractured element, i.e. a mineralised collagen fibril.     

Microcrack accumulation at different intervals during fatigue testing of compact bone

Fatigue damage in bone occurs in the form of microcracks. This microdamage contributes to the formation of stress fractures and acts as a stimulus for bone remodelling. A technique has been developed, which allows microcrack growth to be monitored during the course of a fatigue test by the application of a series of fluorescent chelating agents. Specimens were taken from bovine tibiae and fatigue tested in cyclic compression at a stress range of 80MPa. The specimens were stained before testing with alizarin and up to three other chelating agents were applied during testing to label microcracks formed at different times. Microcracks initiated in interstitial bone in the early part of a specimen's life. Further accumulation of microcracks is then suppressed until the period late in the specimen's life. Microcracks were found to be longer in the longitudinal than in the transverse direction. Only a small proportion of cracks are actively propagating; these are longer than non-propagating cracks. These results support the concept of a microstructural barrier effect existing in bone, whereby cracks initiate easily but slow down or stop at barriers such as cement lines.     

Microdamage accumulation in the cement layer of hip replacements under flexural loading.

Mechanical fatigue of bone cement leading to damage accumulation is implicated in the loosening of cemented hip components. Even though cracks have been identified in autopsy-retrieved mantles, damage accumulation by continuous growth and increase in number of microcracks has not yet been demonstrated experimentally. To determine just how damage accumulation occurs in the cement layer of a hip replacement, a physical model of the joint was used in an experimental study. The model regenerates the stress pattern found in the cement layers whilst at the same time allowing visualisation of microcrack initiation and growth. In this way the gradual process of damage accumulation can be determined. Six specimens were tested to 5 million cycles and a total of 1373 cracks were observed. It was found that, under the flexural loading allowed by the model, the majority of cracks come from pores in the bulk cement and not from the interfaces. Furthermore, the lateral and medial sides have statistically different damage accumulation behaviours, and pre-load cracks significantly accelerate the damage accumulation process. The experimental results confirm that damage accumulation commences early on in the loading history and that it is continuously increasing with load in the form of crack initiation and crack propagation. The results highlight the importance of replicating the loading and restraint conditions of clinical cement mantles when endeavouring to accurately model the damage accumulation process.     

Synchrotron radiation micro-CT at the micrometer scale for the analysis of the three-dimensional morphology of microcracks in human trabecular bone

Bone quality is an important concept to explain bone fragility in addition to bone mass. Among bone quality factors, microdamage which appears in daily life is thought to have a marked impact on bone strength and plays a major role in the repair process. The starting point for all studies designed to further our understanding of how bone microdamage initiate or dissipate energy, or to investigate the impact of age, gender or disease, remains reliable observation and measurement of microdamage. In this study, 3D Synchrotron Radiation (SR) micro-CT at the micrometric scale was coupled to image analysis for the three-dimensional characterization of bone microdamage in human trabecular bone specimens taken from femoral heads. Specimens were imaged by 3D SR micro-CT with a voxel size of 1.4 μm. A new tailored 3D image analysis technique was developed to segment and quantify microcracks. Microcracks from human trabecular bone were observed in different tomographic sections as well as from 3D renderings. New 3D quantitative measurements on the microcrack density and morphology are reported on five specimens. The 3D microcrack density was found between 3.1 and 9.4/mm3 corresponding to a 2D density between 0.55 and 0.76 /mm2. The microcrack length and width measured in 3D on five selected microcrack ranged respectively from 164 μm to 209 μm and 100 μm to 120 μm. This is the first time that various microcracks in unloaded human trabecular bone--from the simplest linear crack to more complex cross-hatch cracks--have been examined and quantified by 3D imaging at this scale. The suspected complex morphology of microcracks is here considerably more evident than in the 2D observations. In conclusion, this technique opens new perspective for the 3D investigation of microcracks and the impact of age, disease or treatment.     

Micromechanics fracture in osteonal cortical bone: a study of the interactions between microcrack propagation, microstructure and the material properties

A two-dimensional micromechanical fibre reinforced composite materials model for osteonal cortical bone is presented. The interstitial bone is modelled as a matrix, the osteons are modelled as fibres, and the cement line is presented as interface tissue. The interaction between osteons and microcracks is evaluated by linear elastic fracture mechanics theory, followed by a determination of the stress intensity factor at the vicinity of the microcrack tips. The results indicate that bone microstructural heterogeneity greatly influences fracture parameters. Furthermore, microstructural morphology and loading conditions affect growth trajectories, the microcrack propagation trajectory deviates from the osteon under tensile loading, and osteon penetration is observed under compressive loads.     

The fatigue strength of compact bone in torsion

We have conducted a series of fatigue tests on samples of bovine compact bone loaded in cyclic torsion. The fatigue strength (i.e. the range of stress needed to cause failure in a given number of cycles) was found to be lower than the fatigue strength of the same material in compression by more than a factor of two. We also tested intact chicken metatarsals and found a similar reduction in strength compared to compression testing of chicken tibiae. These results were predicted using a theoretical model in which fatigue failure was assumed to be dependent on the growth of microcracks, oriented approximately parallel to the bone's longitudinal axis but having misorientation angles of up to 30 degrees. An effective stress range was derived which is a function of the normal and shear stresses, and thus of the Mode I and Mode II stress intensities experienced by the crack. These results may have important consequences for the understanding of fatigue in bone in vivo; relatively small amounts of longitudinal shear stress, which are often ignored in analysis, may contribute significantly to fatigue failures. This may shed light on the phenomenon of stress fractures and on the need for repair and adaptation in living bone.     

Near-terminal creep damage does not substantially influence fatigue life under physiological loading

Cortical bone specimens were damaged using repeated blocks of tensile creep loading until a near-terminal amount of creep damage was generated (corresponding to a reduction in elastic modulus of 15%). One group of cortical bone specimens was submitted to the near-terminal damage protocol and subsequently underwent fatigue loading in tension with a maximum strain of 2000 με (Damage Fatigue, n=5). A second group was submitted to cyclic fatigue loading but was not pre-damaged (Control Fatigue, n=5). All but one specimen (a damaged specimen) reached run-out (10 million cycles, 7.7 days). No significant differences in microscopic cracks or other tissue damage were observed between the two groups or between either group and additional, completely unloaded specimens. Our results suggest that damage in cortical bone allograft that is not obvious or associated with a stress riser may not substantially affect its fatigue life under physiologic loading.