A crack model of a bone cement interface

This paper is concerned with the fracture mechanics of a bone-cement interface that includes a cohesive zone effect on the crack faces. This accounts for the experimentally observed strengthening mechanism due to the mechanical interlock between the crack faces. Edge crack models are developed where the cohesive zone is simulated by a continuous or a discrete distribution of linear or nonlinear springs. It is shown that the solution obtained by assuming a homogeneous material is fairly close to the exact solution for the bimaterial interface edge crack problem. On the basis of that approximation, the analysis is conducted for the problem of two interacting edge cracks, one at the interface, and the other one in the cement. The small crack that was observed to initiate in the cement, close to the bone-cement interface, does not affect much the mode I stress-intensity factor at the tip of the interface crack. However it may grow, leading to a catastrophic breakdown of the cement. The analysis and following discussion point out an interdependency between bone-cement interface strength and cement strength not previously appreciated. The suggested crack models provide a framework for quantifying the fracture mechanisms at the bone-cement interface.     

Cracks emanating from circular voids or elastic inclusions in PMMA near a bone-implant interface

Mechanical fracture is believed to be a primary reason for loss of fixation at the bone-cement-implant interface. In addition to the expected cracks at the bone-cement interface, cracks are also observed to be formed at voids and inclusions within the cement. An analytical solution is presented for cracks emanating from circular voids or elastic inclusions under uniaxial tension using the solution for a single dislocation as a Green's function. Stress intensity factors are calculated for arbitrary orientations of the cracks, and for varying relative stiffnesses of the inclusion and the matrix, to determine the most favorable combination of parameters for crack growth.     

The mechanical effects of different levels of cement penetration at the cement–bone interface

The mechanical effects of varying the depth of cement penetration in the cement–bone interface were investigated using finite element analysis (FEA) and validated using companion experimental data. Two FEA models of the cement–bone interface were created from micro-computed tomography data and the penetration of cement into the bone was varied over six levels each. The FEA models, consisting of the interdigitated cement–bone constructs with friction between cement and bone, were loaded to failure in tension and in shear. The cement and bone elements had provision for crack formation due to excessive stress. The interfacial strength showed a strong relationship with the average interdigitation (r²=0.97 and r²=0.93 in tension and shear, respectively). Also, the interface strength was strongly related with the contact area (r²=0.98 and r²=0.95 in tension and shear, respectively). The FEA results compared favorably to the stiffness–strength relationships determined experimentally. Overall, the cement–bone interface was 2.5 times stronger in shear than in tension and 1.15 times stiffer in tension than in shear, independent of the average interdigitation. More cracks occurred in the cement than in the bone, independent of the average interdigitation, consistent with the experimental results. In addition, more cracks were generated in shear than in tension. In conclusion, achieving and maintaining maximal infiltration of cement into the bone to obtain large interdigitation and contact area is key to optimizing the interfacial strength.     

Mixed-mode loading of the cement-bone interface: a finite element study

While including the cement-bone interface of complete cemented hip reconstructions is crucial to correctly capture their response, its modelling is often overly simplified. In this study, the mechanical mixed-mode response of the cement-bone interface is investigated, taking into account the effects of the well-defined microstructure that characterises the interface. Computed tomography-based plain strain finite element analyses models of the cement-bone interface are built and loaded in multiple directions. Periodic boundaries are considered and the failure of the cement and bone fractions by cracking of the bulk components are included. The results compare favourably with experimental observations. Surprisingly, the analyses reveal that under shear loading no failure occurs and considerable normal compression is generated to prevent interface dilation. Reaction forces, crack patterns and stress fields provide more insight into the mixed-mode failure process. Moreover, the cement-bone interface analyses provide details which can serve as a basis for the development of a cohesive law.     

Micro-mechanical damage of trabecular bone–cement interface under selected loading conditions: a finite element study

In this study, two micro finite element models of trabecular bone–cement interface developed from high resolution computed tomography (CT) images were loaded under compression and validated using the in situ experimental data. The models were then used under tension and shear to examine the load transfer between the bone and cement and the micro damage development at the bone–cement interface. In addition, one models was further modified to investigate the effect of cement penetration on the bone–cement interfacial behaviour. The simulated results show that the load transfer at the bone–cement interface occurred mainly in the bone cement partially interdigitated region, while the fully interdigitated region seemed to contribute little to the mechanical response. Consequently, cement penetration beyond a certain value would seem to be ineffective in improving the mechanical strength of trabecular bone–cement interface. Under tension and shear loading conditions, more cement failures were found in denser bones, while the cement damage is generally low under compression.     

The effect of cement creep and cement fatigue damage on the micromechanics of the cement–bone interface

The cement–bone interface provides fixation for the cement mantle within the bone. The cement–bone interface is affected by fatigue loading in terms of fatigue damage or microcracks and creep, both mostly in the cement. This study investigates how fatigue damage and cement creep separately affect the mechanical response of the cement–bone interface at various load levels in terms of plastic displacement and crack formation. Two FEA models were created, which were based on micro-computed tomography data of two physical cement–bone interface specimens. These models were subjected to tensile fatigue loads with four different magnitudes. Three deformation modes of the cement were considered: ‘only creep’, ‘only damage’ or ‘creep and damage’. The interfacial plastic deformation, the crack reduction as a result of creep and the interfacial stresses in the bone were monitored. The results demonstrate that, although some models failed early, the majority of plastic displacement was caused by fatigue damage, rather than cement creep. However, cement creep does decrease the crack formation in the cement up to 20%. Finally, while cement creep hardly influences the stress levels in the bone, fatigue damage of the cement considerably increases the stress levels in the bone. We conclude that at low load levels the plastic displacement is mainly caused by creep. At moderate to high load levels, however, the plastic displacement is dominated by fatigue damage and is hardly affected by creep, although creep reduced the number of cracks in moderate to high load region.     

The behavior of the micro-mechanical cement-bone interface affects the cement failure in total hip replacement

In the current study, the effects of different ways to implement the complex micro-mechanical behavior of the cement-bone interface on the fatigue failure of the cement mantle were investigated. In an FEA-model of a cemented hip reconstruction the cement-bone interface was modeled and numerically implemented in four different ways: (I) as infinitely stiff, (II) as infinitely strong with a constant stiffness, (III) a mixed-mode failure response with failure in tension and shear, and (IV) realistic mixed mode behavior obtained from micro-FEA models. Case II, III, and IV were analyzed using data from a stiff and a compliant micro-FEA model and their effects on cement failure were analyzed. The data used for Case IV was derived from experimental specimens that were tested previously. Although the total number of cement cracks was low for all cases, the compliant Case II resulted in twice as many cracks as Case I. All cases caused similar stress distributions at the interface. In all cases, the interface did not display interfacial softening; all stayed the elastic zone. Fatigue failure of the cement mantle resulted in a more favorable stress distribution at the cement-bone interface in terms of less tension and lower shear tractions. We conclude that immediate cement-bone interface failure is not likely to occur, but its local compliancy does affect the formation of cement cracks. This means that at a macro-level the cement-bone interface should be modeled as a compliant layer. However, implementation of interfacial post-yield softening does seems to be necessary.     

Surface pretreatment for prolonged survival of cemented tibial prosthesis components: full- vs. surface-cementation technique]

Aseptic loosening of tibial components due to degradation of the interface between bone cement and metallic tibial shaft component is still a persistent problem, particularly for surface-cemented tibial components. The surface cementation technique has important clinical meaning in case of revision and for avoidance of stress shielding. This study was done to prove crack formation in the bone cement near the metallic surface when this is not coated. We propose a newly developed coating process by SiOx-PVD layering to avoid crack formation. A biomechanical model for a vibration fatigue test was done to prove that crack formation can be significantly reduced in the case of coated surfaces. It was found that coated tibial components showed a highly significant reduction of cement cracking near the metal/bone cement interface (p < 0.01) and a significant reduction of gap formation in the metal-to-bone cement interface (p < 0.05). Coating dramatically reduces hydrolytic- and stress-related crack formation at the prosthesis metal/bone cement interface. This leads to a more homogenous load transfer into the cement mantle which should reduce the frequency of loosening in the metal/bone cement/bone interfaces. With surface coating of the tibial component it should become possible that surface-cemented TKAs reveal similar loosening rates as TKAs both surface- and stem-cemented. This would be an important clinical advantage since it is believed that surface cementing reduces metaphyseal bone loss in the case of revision and stress shielding for a better bone health.     

Fatigue creep damage at the cement–bone interface: An experimental and a micro-mechanical finite element study

The goal of this study was to quantify the micromechanics of the cement–bone interface under tensile fatigue loading using finite element analysis (FEA) and to understand the underlying mechanisms that play a role in the fatigue behavior of this interface. Laboratory cement–bone specimens were subjected to a tensile fatigue load, while local displacements and crack growth on the specimen's surface were monitored. FEA models were created from these specimens based upon micro-computed tomography data. To accurately model interfacial gaps at the interface between the bone and cement, a custom-written erosion algorithm was applied to the bone model. A fatigue load was simulated in the FEA models while monitoring the local displacements and crack propagation. The results showed the FEA models were able to capture the general experimental creep damage behavior and creep stages of the interface. Consistent with the experiments, the majority of the deformation took place at the contact interface. Additionally, the FEA models predicted fatigue crack patterns similar to experimental findings. Experimental surface cracks correlated moderately with FEA surface cracks (r²=0.43), but did not correlate with the simulated crack volume fraction (r²=0.06). Although there was no relationship between experimental surface cracks and experimental creep damage displacement (r²=0.07), there was a strong relationship between the FEA crack volume fraction and the FEA creep damage displacement (r²=0.76). This study shows the additional value of FEA of the cement–bone interface relative to experimental studies and can therefore be used to optimize its mechanical properties.     

Material changes in osteoporotic human cancellous bone following infiltration with acrylic bone cement for a vertebral cement augmentation

Bone cement infiltration can be effective at mechanically augmenting osteoporotic vertebrae. While most published literature describes the gain in mechanical strength of augmented vertebrae, we report the first measurements of viscoelastic material changes of cancellous bone due to cement infiltration. We infiltrated cancellous core specimen harvested from osteoporotic cadaveric spines with acrylic bone cement. Bone specimen before and after cement infiltration were subjected to identical quasi-static and relaxation loading in confined and free compression. Testing data were fitted to a linear viscoelastic model of compressible material and the model parameters for cement, native cancellous bone, and cancellous bone infiltrated (composite) with cement were identified. The fitting demonstrated that the linear viscoelastic model presented in this paper accurately describes the mechanical behaviour of cement and bone, before and after infiltration. Although the composite specimen did not completely adopt the properties of bulk bone cement, the stiffening of cancellous bone due to cement infiltration is considerable. The composite was, for example, 8.5 times stiffer than native bone. The local stiffening of cancellous bone in patients may alter the load transfer of the augmented motion segment and may be the cause of subsequent fractures in the vertebrae adjacent to the ones infiltrated with cement. The material model and parameters in this paper, together with an adequate finite-element model, can be helpful to investigate the load shift, the mechanism for subsequent fractures, and filling patterns for ideal cement infiltration.     

Cracks emanating from slipping inclusions near a bone-implant interface.

It has been hypothesized that mechanical fracture at the bone-cement-implant interface is the initial cause for loosening of orthopedic implants. Previous investigators have observed cracks to emanate from methacrylate beads, apparently acting as inclusions within the cement. It is believed that the bond between these inclusions and the surrounding matrix breaks prior to emanation of radial cracks from the inclusion. An analytical model is developed for radial cracks emanating from circular inclusions that allow slip along their interface. The solution to the interaction of a single dislocation and a slipping inclusion is used as a Green's Function to model the crack. The Mode I stress intensity factors are calculated for arbitrary orientations of the crack and for varying relative stiffness of the matrix and the inclusion to test feasibility of crack growth.     

Material Changes in Osteoporotic Human Cancellous Bone Following Infiltration with Acrylic Bone Cement for a Vertebral Cement Augmentation

Bone cement infiltration can be effective at mechanically augmenting osteoporotic vertebrae. While most published literature describes the gain in mechanical strength of augmented vertebrae, we report the first measurements of viscoelastic material changes of cancellous bone due to cement infiltration. We infiltrated cancellous core specimen harvested from osteoporotic cadaveric spines with acrylic bone cement. Bone specimen before and after cement infiltration were subjected to identical quasi-static and relaxation loading in confined and free compression. Testing data were fitted to a linear viscoelastic model of compressible material and the model parameters for cement, native cancellous bone, and cancellous bone infiltrated (composite) with cement were identified. The fitting demonstrated that the linear viscoelastic model presented in this paper accurately describes the mechanical behaviour of cement and bone, before and after infiltration. Although the composite specimen did not completely adopt the properties of bulk bone cement, the stiffening of cancellous bone due to cement infiltration is considerable. The composite was, for example, 8.5 times stiffer than native bone. The local stiffening of cancellous bone in patients may alter the load transfer of the augmented motion segment and may be the cause of subsequent fractures in the vertebrae adjacent to the ones infiltrated with cement. The material model and parameters in this paper, together with an adequate finite-element model, can be helpful to investigate the load shift, the mechanism for subsequent fractures, and filling patterns for ideal cement infiltration.     

Residual stress due to curing can initiate damage in porous bone cement: experimental and theoretical evidence

Residual stress due to shrinkage of polymethylmethacrylate bone cement after polymerisation is possibly one factor capable of initiating cracks in the mantle of cemented hip replacements. No relationship between residual stress and observed cracking of cement has yet been demonstrated. To investigate if any relationship exists, a physical model has been developed which allows direct observation of damage in the cement layer on the femoral side of total hip replacement. The model contains medial and lateral cement layers between a bony surface and a metal stem; the tubular nature of the cement mantle is ignored. Five specimens were prepared and examined for cracking using manual tracing of stained cracks, observed by transmission microscopy; cracks were located and measured using image analysis. A mathematical approach for the prediction of residual stress due to shrinkage was developed which uses the thermal history of the material to predict when stress-locking occurs, and estimates subsequent thermal stress. The residual stress distribution of the cement layer in the physical model was then calculated using finite element analysis. Results show maximum tensile stresses normal to the observed crack directions, suggesting a link between residual stress and pre-load cracking. The residual stress predicted depends strongly on the definition of the reference temperature for stress-locking. The highest residual stresses (4-7 MPa) are predicted for shrinkage from maximum temperature; in this case, magnitudes are sufficiently high to initiate cracks when the influence of stress raisers such as pores or interdigitation at the bone/cement interface are taken into account (up to 24 MPa when calculating stress around a pore according to the method of Harrigan and Harris (J. Biomech. 24(11) (1991) 1047-1058). We conclude that the damage accumulation failure scenario begins before weight-bearing due to cracking induced by residual stress around pores or stress raisers.     

The influence of cement mantle thickness and stem geometry on fatigue damage in two different cemented hip femoral prostheses

Experimental models can be used for pre-clinical testing of cemented and other type of hip replacements. Total hip replacement (THR) failure scenarios include, among others, cement damage accumulation and the assessment of accurate stress and strain magnitudes at the cement mantle interfaces (stem-cement and cement-bone) can be used to predict mechanical failure. The aseptic loosening scenario in cemented hip replacements is currently not fully understood, and methods of evaluating medical devices must be developed to improve clinical performance. Different results and conclusions concerning the cement micro-cracking mechanism have been reported.The aim of this study was to verify the in vitro behavior of two cemented femoral stems with respect to fatigue crack formation. Fatigue crack damage was assessed at the medial, lateral, anterior and posterior sides of the Lubinus SPII and Charnley stems. All stems were loaded and tested in stair climbing fatigue loading during one million cycles at 2 Hz. After the experiments each implanted synthetic femur was sectioned and analyzed. We observed more damage (cracks per area) for the Lubinus SPII stem, mainly on the proximal part of the cement mantle. The micro-cracking formation initiated in the stem–cement interface and grew towards the direction of cortical bone of the femur.Overall, the cement–bone interface seems to be crucial for the success of the hip replacement. The Charnley stem provoked more damage on the cement–bone interface. A failure index (maximum length of crack/maximum thickness of cement) considered was higher for the cement–stem interface of the Lubinus SPII stem. For a cement mantle thickness higher than 5 mm, cracking initiated at the cement–bone interface and depended on the opening canal process (reaming procedure and instrumentation). The analysis also showed that fatigue-induced damage on the cement mantle, increasing proximally, and depended on the axial position of the stem. The cement thickness is an important factor for the success of THR and this study evidenced that cement thickness higher than 2 mm apparently does not affect the mechanical behavior of the cement mantel and induce more crack formation on the cement–bone interface.     

Cohesive finite element modeling of age-related toughness loss in human cortical bone

Although the age-related loss of bone quality has been implicated in bone fragility, a mechanistic understanding of the relationship is necessary for developing diagnostic and treatment modalities in the elderly population at risk of fracture. In this study, a finite element based cohesive zone model is developed and applied to human cortical bone in order to capture the experimentally shown rising crack growth behavior and age-related loss of bone toughness. The cohesive model developed here is based on a traction–crack opening displacement relationship representing the fracture processes in the vicinity of a propagating crack. The traction–displacement curve, defining the cohesive model, is composed of ascending and descending branches that incorporate material softening and nonlinearity. The results obtained indicate that, in contrast to initiation toughness, the finite element simulations of crack growth in compact tension (CT) specimens successfully capture the rising R-curve (propagation toughness) behavior and the age-related loss of bone toughness. In close correspondence with the experimentally observed decrease of 14–15% per decade, the finite element simulation results show a decrease of 13% in the R-curve slope per decade. The success of the simulations is a result of the ability of cohesive models to capture and predict the parameters related to bone fracture by representing the physical processes occurring in the vicinity of a propagating crack. These results illustrate that fracture mechanisms in the process zone control bone toughness and any modification to these would cause age-related toughness loss.     

Fracture analysis for biological materials with an expanded cohesive zone model

In this study, a theoretical framework for simulation of fracture of bone and bone-like materials is provided. An expanded cohesive zone model with thermodynamically consistent framework has been proposed and used to investigate the crack growth resistance of bone and bone-like materials. The reversible elastic deformation, irreversible plastic deformation caused by large deformation of soft protein matrix, and damage evidenced by the material separation and crack nucleation in the cohesive zone, were all taken into account in the model. Furthermore, the key mechanisms in deformation of biocomposites consisting of mineral platelets and protein interfacial layers were incorporated in the fracture process zone in this model, thereby overcoming the limitations of previous cohesive zone modeling of bone fracture. Finally, applications to fracture of cortical bone and human dentin were presented, which showed good agreement between numerical simulation and reported experiments and substantiated the effectiveness of the model in investigating the fracture behavior of bone-like materials.     

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.     

Crack propagation in cortical bone is affected by the characteristics of the cement line: a parameter study using an XFEM interface damage model

Bulk properties of cortical bone have been well characterized experimentally, and potent toughening mechanisms, e.g., crack deflections, have been identified at the microscale. However, it is currently difficult to experimentally measure local damage properties and isolate their effect on the tissue fracture resistance. Instead, computer models can be used to analyze the impact of local characteristics and structures, but material parameters required in computer models are not well established. The aim of this study was therefore to identify the material parameters that are important for crack propagation in cortical bone and to elucidate what parameters need to be better defined experimentally. A comprehensive material parameter study was performed using an XFEM interface damage model in 2D to simulate crack propagation around an osteon at the microscale. The importance of 14 factors (material parameters) on four different outcome criteria (maximum force, fracture energy, crack length and crack trajectory) was evaluated using ANOVA for three different osteon orientations. The results identified factors related to the cement line to influence the crack propagation, where the interface strength was important for the ability to deflect cracks. Crack deflection was also favored by low interface stiffness. However, the cement line properties are not well determined experimentally and need to be better characterized. The matrix and osteon stiffness had no or low impact on the crack pattern. Furthermore, the results illustrated how reduced matrix toughness promoted crack penetration of the cement line. This effect is highly relevant for the understanding of the influence of aging on crack propagation and fracture resistance in cortical bone.

     

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.     

Anisotropic post-yield response of cancellous bone simulated by stress–strain curves of bulk equivalent structures

During the last decade, finite element (FE) modelling has become ubiquitous in understanding complex mechanobiological phenomena, e.g. bone–implant interactions. The extensive computational effort required to achieve biorealistic results when modelling the post-yield behaviour of microstructures like cancellous bone is a major limitation of these techniques. This study describes the anisotropic biomechanical response of cancellous bone through stress–strain curves of equivalent bulk geometries. A cancellous bone segment, reverse engineered by micro computed tomography, was subjected to uniaxial compression. The material's constitutive law, obtained by nano-indentations, was considered during the simulation of the experimental process. A homodimensionally bulk geometry was employed to determine equivalent properties, resulting in a similar anisotropic response to the trabecular structure. The experimental verification of our model sustained that the obtained stress–strain curves can adequately reflect the post-yield behaviour of the sample. The introduced approach facilitates the consideration of nonlinearity and anisotropy of the tissue, while reducing the geometrical complexity of the model to a minimum.