Age-related properties at the microscale affect crack propagation in cortical bone

Research on ultrasonic bone cutting mechanism based on extended finite element method

The research on the crack propagation mechanism of bone has important research significance and clinical medical value for the selection of cutting parameters and the development of new surgical tools. In this paper, an extended finite element method (X-FEM) model of ultrasonic bone cutting considering microstructure was developed to further study the ultrasonic bone cutting mechanism and to quantitatively analyze the effects of cutting direction, ultrasonic parameters, and cutting parameters on the mechanism of ultrasonic bone cutting crack propagation. The results show that ultrasonic bone cutting is essentially a controlled crack propagation process, in which brittle crack and fatigue crack are the main crack propagation mechanisms. In order to improve the efficiency of ultrasonic bone cutting, large amplitude and high-frequency ultrasonic vibration are preferred. Compared with the other two cutting directions, the crack propagation deflection angle in the transverse cutting direction is the largest, resulting in the worst cutting surface. Therefore, in the path planning of orthopedic surgical robots, the transverse cutting direction should be avoided as much as possible. Frequency only has a significant effect on the crack propagation rate and has a positive correlation. There is a positive correlation between the deflection angle, propagation length, propagation rate, and amplitude, which provides the possibility to control the direction and length of crack propagation by controlling the amplitude of ultrasonic. The feed speed is much lower than the ultrasonic vibration speed, which makes the influence of ultrasonic vibration speed on the crack propagation characteristics dominant. The X-FEM model of ultrasonic bone cutting provides an effective method for selecting reasonable machining parameters of orthopedic robot and optimize the design of ultrasonic osteotome.

A numerical study of dehydration induced fracture toughness degradation in human cortical bone

A 2D plane strain extended finite element method (XFEM) model was developed to simulate three-point bending fracture toughness tests for human bone conducted in hydrated and dehydrated conditions. Bone microstructures and crack paths observed by micro-CT imaging were simulated using an XFEM damage model. Critical damage strains for the osteons, matrix, and cement lines were deduced for both hydrated and dehydrated conditions and it was found that dehydration decreases the critical damage strains by about 50%. Subsequent parametric studies using the various microstructural models were performed to understand the impact of individual critical damage strain variations on the fracture behavior. The study revealed the significant impact of the cement line critical damage strains on the crack paths and fracture toughness during the early stages of crack growth. Furthermore, a significant sensitivity of crack growth resistance and crack paths on critical strain values of the cement lines was found to exist for the hydrated environments where a small change in critical strain values of the cement lines can alter the crack path to give a significant reduction in fracture resistance. In contrast, in the dehydrated state where toughness is low, the sensitivity to changes in critical strain values of the cement lines is low. Overall, our XFEM model was able to provide new insights into how dehydration affects the micromechanisms of fracture in bone and this approach could be further extended to study the effects of aging, disease, and medical therapies on bone fracture.

An integrated experimental-computational framework to assess the influence of microstructure and material properties on fracture toughness in clinical specimens of human femoral cortical bone

Microstructural and compositional changes that occur due to aging, pathological conditions, or pharmacological treatments alter cortical bone fracture resistance. However, the relative importance of these changes to the fracture resistance of cortical bone has not been quantified in detail. In this technical note, we developed an integrated experimental-computational framework utilizing human femoral cortical bone biopsies to advance the understanding of how fracture resistance of cortical bone is modulated due to modifications in its microstructure and material properties. Four human biopsy samples from individuals with varying fragility fracture history and osteoporosis treatment status were converted to finite element models incorporating specimen-specific material properties and were analyzed using fracture mechanics-based modeling. The results showed that cement line density and osteonal volume had a significant effect on crack volume. The removal of cement lines substantially increased the crack volume in the osteons and interstitial bone, representing straight crack growth, compared to models with cement lines due to the lack of crack deflection in the models without cement lines. Crack volume in the osteons and interstitial bone increased when mean elastic modulus and ultimate strength increased and mean fracture toughness decreased. Crack volume in the osteons and interstitial bone was reduced when material property heterogeneity was incorporated in the models. Although both the microstructure and the heterogeneity of the material properties of the cortical bone independently increased the fracture toughness, the relative contribution of the microstructure was more significant. The integrated experimental-computational framework developed here can identify the most critical microscale features of cortical bone modulated by pathological processes or pharmacological treatments that drive changes in fracture resistance and improve our understanding of the relative influence of microstructure and material properties on fracture resistance of cortical bone.

A review on experimental and numerical investigations of cortical bone fracture

This paper comprehensively reviews the various experimental and numerical techniques, which were considered to determine the fracture characteristics of the cortical bone. This study also provides some recommendations along with the critical review, which would be beneficial for future research of fracture analysis of cortical bone. Cortical bone fractures due to sports activities, climbing, running, and engagement in transport or industrial accidents. Individuals having different diseases are also at high risk of cortical bone fracture. It has been observed that osteon orientation influences cortical bone fracture toughness and fracture mechanisms. Apart from this, recent studies indicate that fracture parameters of cortical bone also depend on many factors such as age, sex, temperature, osteoporosis, orientation, location, loading condition, strain rate, and storage facility, etc. The cortical bone regains its fracture toughness due to various toughening mechanisms. Owing to these factors, several experimental, clinical, and numerical investigations have been carried out to determine the fracture parameters of the cortical bone. Cortical bone is the dense outer surface of the bone and contributes to 80%–82% of the skeleton mass. Cortical bone experiences load far exceeding body weight due to muscle contraction and the dynamics of motion. It is very important to know the fracture pattern, direction of fracture, location of the fracture, and toughening mechanism of cortical bone. A basic understanding of the different factors that affect the fracture parameters and fracture mechanisms of the cortical bone is necessary to prevent the failure and fracture of cortical bone. This review has summarized the advancement considered in the various experimental techniques and numerical methods to get complete information about the fracture mechanisms of cortical bone.

Damage analysis of human cortical bone under compressive and tensile loadings

Developing advanced fracture tools can increase the understanding of crack growth trajectories in human cortical bone. The present study investigates fracture micromechanics of human cortical bone under compressive and tensile loadings utilizing a phase field method. We construct two-dimensional finite element models from cortical microstructure of a human tibia cross section. We apply compression on the cortical bone models to create compressive microcracks. Then, we simulate the fracture of these models under tension to discover influential parameters on microcracks formation and post-yielding behavior. The results show that cement lines are susceptible sites to damage nucleation under compression rather than tension. The findings of this study also indicate a higher accumulation of initial damage (induced by compression) can lead to a lower microscopic stiffness as well as a less resistant material to damage initiation under tension. The simulations further indicate that the post-yielding properties (e.g., toughness) can be dependent on different variables such as morphological information of the osteons, the initial accumulation of microcracks, and the total length of cement lines.

Influence of Porosity on Fracture Toughness and Fracture Behavior of Antibiotic-Loaded PMMA Bone Cement

Aseptic loosening is the most common reason for the long-term revision of cemented arthroplasties with fracture of the cement being a postulated cause or contributing factor. In our previous studies we showed that adding an antibiotic to a polymethylmethacrylate (PMMA) bone cement led to detrimental effects on various mechanical properties of the cement such as bending strength, compressive strength and fracture toughness (KIC). This finding implied that the mechanical failure of antibiotic-loaded PMMA bone cement was influenced by its pore volume fraction. Up to now this aspect has not been studied. Hence the purposes of this study were to determine (1) the influence of antibiotic (telavancin) loading on the KIC of a widely used PMMA bone cement brand (Palacos®R) and (2) the influence of pore size and pore distribution on the fracture behavior of the KIC specimens. For (2) both experimental and numerical methods (extended finite element method [XFEM]) were used allowing a comparison between the two sets of results. We found that: (1) KIC decreased with increased porosity with the drop (relative to the value for the control cement) being significant when the telavancin loading was 4.8 wt/wt % (2 g of telavancin added to 40 g of control cement powder); (2) there was a critical pore size above which there was a significant decrease in KIC and is 1 mm; (3) crack propagation was strongly influenced by pore size and pore locations (pore-pore interactions); and, (4) there was good agreement between the experimental and XFEM results. The implications of these findings for the use of a telavancin-loaded PMMA bone cement in cemented total joint arthroplasties are commented upon.

Computational homogenisation based extraction of transverse tensile cohesive responses of cortical bone tissue

The numerical assessment of fracture properties of cortical bone is important in providing suggestions on patient-specific clinical treatments. We present a generic finite element modelling framework incorporating computational fracture approaches and computational homogenisation techniques. Finite element computations for statistical volume elements (SVEs) at the microscale are performed for different sizes with random osteon packing with a fixed volume fraction. These SVEs are loaded in the transverse direction under tension. The minimal SVE size in terms of ensuring a representative effective cohesive law is suggested to be 0.6 mm. Since cement lines as weak interfaces play a key role in bone fracture, the effects of their fracture properties on the effective fracture strength and toughness are investigated. The extracted effective fracture properties can be used as homogenised inputs to a discrete crack simulation at macroscopic or structural scale. The extrinsic toughening mechanisms observed in the SVE models are discussed with a comparison against experimental observations from the literature, giving beneficial insights to cortical bone failure.

Serum exosomes from young rats improve the reduced osteogenic differentiation of BMSCs in aged rats with osteoporosis after fatigue loading in vivo

BACKGROUND

Osteoporosis is a major public health concern for the elderly population and is characterized by fatigue load resulting in bone microdamage. The ability of bone mesenchymal stem cells (BMSCs) to repair bone microdamage diminishes with age, and the accumulation of bone microdamage increases the risk of osteoporotic fracture. There is a lack of effective means to promote the repair of bone microdamage in aged patients with osteoporosis. Exosomes have been shown to be related to the osteogenic differentiation of BMSCs. Here, we aimed to evaluate the changes in the osteogenic differentiation capacity of BMSCs in aged osteoporotic rats after fatigue loading and the treatment potential of serum exosomes from young rats.

METHODS

The tibias of six aged osteoporotic rats were subjected to fatigue loading in vivo for 2 weeks, and the bone microdamage, microstructures, and mechanical properties were assessed. Subsequently, BMSCs were extracted to evaluate their proliferation and osteogenic differentiation abilities. In addition, the BMSCs of aged osteoporotic rats after fatigue loading were treated with serum exosomes from young rats under osteogenic induction conditions, and the expression of osteogenic-related miRNAs was quantified. The osteogenetic effects of miRNA-19b-3p in exosomes and the possible target protein PTEN was detected.

RESULTS

Obvious bone microdamage at the fatigue load stress point, the bone microstructure and biomechanical properties were not obviously changed. A decreased osteogenic differentiation ability of BMSCs was observed after fatigue loading, while serum exosomes from young rats highly expressing miRNA-19b-3p improved the decreased osteogenic differentiation ability of BMSCs. Transfection with miRNA-19b-3p mimic could promote osteoblastic differentiation of BMSCs and decreased the expression of PTEN. After transfection of miRNA-19b-3p inhibitor, the promotional effect of exosomes on bone differentiation was weakened. Treatment with transfected exosomes increased the expression of PTEN.

CONCLUSION

Serum exosomes derived from young rats can improve the decreased osteogenic differentiation ability of BMSCs in aged rats with osteoporosis after fatigue loading and can provide a new treatment strategy for the repair of bone microdamage and prevention of fractures.

Fracture behavior of human cortical bone: Role of advanced glycation end-products and microstructural features

Diabetes is associated with increased fracture risk in human bone, especially in the elderly population. In the present study, we investigate how simulated advanced glycation end-products (AGEs) and materials heterogeneity affect crack growth trajectory in human cortical bone. We used a phase field fracture framework on 2D models of cortical microstructure created from human tibias to analyze crack propagation. The increased AGEs level results in a higher rate of crack formation. The simulations also indicate that the mismatch between the fracture properties (e.g., critical energy release rate) of osteons and interstitial tissue can alter the post-yielding behavior. The results show that if the critical energy release rate of cement lines is lower than that of osteons and the surrounding interstitial matrix, cracks can be arrested by cement lines. Additionally, activation of toughening mechanisms such as crack merging and branching depends on bone microstructural morphology (i.e., osteons geometrical parameters, canals, and lacunae porosities). In conclusion, the present findings suggest that materials heterogeneity of microstructural features and the crack-microstructure interactions can play important roles in bone fragility.

Effect of ageing on microstructure and fracture behavior of cortical bone as determined by experiment and Extended Finite Element Method (XFEM)

Bone fracture is a severe health concern; therefore, understanding the causes of bone fracture are crucial. This paper investigates the microstructure and fracture behaviour of cadaveric cortical bone of two different groups (Young, n= 6; Aged, n=7). The microstructure is obtained from µ-CT images, and the material parameters are measured with nanoindentation. Fracture behaviour in transverse and longitudinal orientations is investigated experimentally and numerically. The results show that the Haversian canal (HC) size increases and the osteon wall thickness (OWT) decreases significantly in the aged group, whereas a nonsignificant difference is found in tissue properties. The crack initiation (J_ic) and crack growth (J_grow) toughness of the aged group are found to be significantly lower (p<0.01) than the young group in the transverse orientation; however, for the longitudinal orientation, only the value of J_ic in the aged group is found significantly lower. Further, a 4-phase XFEM (based on micro-CT image) model is developed to investigate the crack propagation behaviour in both orientations. For the transverse orientation, results show that in the aged group, the crack initially follows the cementline and then penetrates the osteon, whereas, in the young group, it propagates along the cementline. These results are in agreement with experimental results where the decrease in J_grow is more significant than the J_ic in the aged group. This study suggests that ageing leads to a larger HC and reduced OWT, which weakens the crack deflection ability and causes fragility fracture. Further, the XFEM results indicate that the presence of a small microcrack in the vicinity of a major crack tip causes an increase in the critical stress intensity factor.

The influence of microstructure on crack propagation in cortical bone at the mesoscale

The microstructure of cortical bone is key for the tissue's high toughness and strength and efficient toughening mechanisms have been identified at the microscale, for example when propagating cracks interact with the osteonal microstructure. Finite element models have been proposed as suitable tools for analyzing the complex link between the local tissue structure and the fracture resistance of cortical bone. However, previous models that could capture realistic crack paths in cortical bone were due to the required computational effort limited to idealized osteon geometries and small (<1 mm²) model domains. The objective of this study was therefore to bridge the gap between experimental and numerical analysis of crack propagation in cortical bone by introducing image-based models at the mesoscale. Tissue orientation maps from high-resolution micro-CT images were used to define the distribution and orientation of weak interfaces in the models. Crack propagation was simulated using the extended finite element method in combination with an interface damage model, previously developed to simulate crack propagation in microstructural osteon models. The results showed that image-based mesoscale models can be used to capture interactions between cracks and microstructure. The simulated crack paths predicted the general trends seen in experiments with more irregular patterns for cracks propagating perpendicular compared to parallel to the osteon orientation. In all, the proposed method enabled an efficient description of the tissue level microstructure, which is a necessity to predict realistic crack paths in cortical bone and is an important step towards simulating crack propagation in bone models in 3D.

Advanced Modeling Methods-Applications to Bone Fracture Mechanics

PURPOSE OF REVIEW

The goal of this review is to summarize recent advances in modeling of bone fracture using fracture mechanics-based approaches at multiple length scales spanning nano- to macroscale.

RECENT FINDINGS

Despite the additional information that fracture mechanics-based models provide over strength-based ones, the application of this approach to assessing bone fracture is still somewhat limited. Macroscale fracture models of bone have demonstrated the potential of this approach in uncovering the contributions of geometry, material property variation, as well as loading mode and rate on whole bone fracture response. Cortical and cancellous microscale models of bone have advanced the understanding of individual contributions of microstructure, microarchitecture, local material properties, and material distribution on microscale fracture resistance of bone. Nano/submicroscale models have provided additional insight into the effect of specific changes in mineral, collagen, and non-collagenous proteins as well as their interaction on energy dissipation and fracture resistance at small length scales. Advanced modeling approaches based on fracture mechanics provide unique information about the underlying multiscale fracture mechanisms in bone and how these mechanisms are influenced by the structural and material constituents of bone at different length scales. Fracture mechanics-based modeling provides a powerful approach that complements experimental evaluations and advances the understanding of critical determinants of fracture risk.

Subject-specific FE models of the human femur predict fracture path and bone strength under single-leg-stance loading

Hip fractures are a major health problem with high socio-economic costs. Subject-specific finite element (FE) models have been suggested to improve the fracture risk assessment, as compared to clinical tools based on areal bone mineral density, by adding an estimate of bone strength. Typically, such FE models are limited to estimate bone strength and possibly the fracture onset, but do not model the fracture process itself. The aim of this study was to use a discrete damage approach to simulate the full fracture process in subject-specific femur models under stance loading conditions. A framework based on the partition of unity finite element method (PUFEM), also known as XFEM, was used. An existing PUFEM framework previously used on a homogeneous generic femur model was extended to include a heterogeneous material description together with a strain-based criterion for crack initiation. The model was tested on two femurs, previously mechanically tested in vitro. Our results illustrate the importance of implementing a subject-specific material distribution to capture the experimental fracture pattern under stance loading. Our models accurately predicted the fracture pattern and bone strength (1% and 5% error) in both investigated femurs. This is the first study to simulate complete fracture paths in subject-specific FE femur models and it demonstrated how discrete damage models can provide a more complete picture of fracture risk by considering both bone strength and fracture toughness in a subject-specific fashion.