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

Multiscale and multidisciplinary analysis of aging processes in bone

The world population is increasingly aging, deeply affecting our society by challenging our healthcare systems and presenting an economic burden, thus turning the spotlight on aging-related diseases: exempli gratia, osteoporosis, a silent disease until you suddenly break a bone. The increase in bone fracture risk with age is generally associated with a loss of bone mass and an alteration in the skeletal architecture. However, such changes cannot fully explain increased fragility with age. To successfully tackle age-related bone diseases, it is paramount to comprehensively understand the fundamental mechanisms responsible for tissue degeneration. Aging mechanisms persist at multiple length scales within the complex hierarchical bone structure, raising the need for a multiscale and multidisciplinary approach to resolve them. This paper aims to provide an overarching analysis of aging processes in bone and to review the most prominent outcomes of bone aging. A systematic description of different length scales, highlighting the corresponding techniques adopted at each scale and motivating the need for combining diverse techniques, is provided to get a comprehensive description of the multi-physics phenomena involved.

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.

XFEM for Composites, Biological, and Bioinspired Materials: A Review

The eXtended finite element method (XFEM) is a powerful tool for structural mechanics, assisting engineers and designers in understanding how a material architecture responds to stresses and consequently assisting the creation of mechanically improved structures. The XFEM method has unraveled the extraordinary relationships between material topology and fracture behavior in biological and engineered materials, enhancing peculiar fracture toughening mechanisms, such as crack deflection and arrest. Despite its extensive use, a detailed revision of case studies involving XFEM with a focus on the applications rather than the method of numerical modeling is in great need. In this review, XFEM is introduced and briefly compared to other computational fracture models such as the contour integral method, virtual crack closing technique, cohesive zone model, and phase-field model, highlighting the pros and cons of the methods (e.g., numerical convergence, commercial software implementation, pre-set of crack parameters, and calculation speed). The use of XFEM in material design is demonstrated and discussed, focusing on presenting the current research on composites and biological and bioinspired materials, but also briefly introducing its application to other fields. This review concludes with a discussion of the XFEM drawbacks and provides an overview of the future perspectives of this method in applied material science research, such as the merging of XFEM and artificial intelligence techniques.

Hormonal and non-hormonal oral contraceptives given long-term to pubertal rats differently affect bone mass, quality and metabolism

Introduction

We investigated the effects of hormonal and non-hormonal oral contraceptives (OCs) on bone mass, mineralization, composition, mechanical properties, and metabolites in pubertal female SD rats.

Methods

OCs were given for 3-, and 7 months at human equivalent doses. The combined hormonal contraceptive (CHC) was ethinyl estradiol and progestin, whereas the non-hormonal contraceptive (NHC) was ormeloxifene. MicroCT was used to assess bone microarchitecture and BMD. Bone formation and mineralization were assessed by static and dynamic histomorphometry. The 3-point bending test, nanoindentation, FTIR, and cyclic reference point indentation (cRPI) measured the changes in bone strength and material composition. Bone and serum metabolomes were studied to identify potential biomarkers of drug efficacy and safety and gain insight into the underlying mechanisms of action of the OCs.

Results

NHC increased bone mass in the femur metaphysis after 3 months, but the gain was lost after 7 months. After 7 months, both OCs decreased bone mass and deteriorated trabecular microarchitecture in the femur metaphysis and lumbar spine. Also, both OCs decreased the mineral: matrix ratio and increased the unmineralized matrix after 7 months. After 3 months, the OCs increased carbonate: phosphate and carbonate: amide I ratios, indicating a disordered hydroxyapatite crystal structure susceptible to resorption, but these changes mostly reversed after 7 months, indicating that the early changes contributed to demineralization at the later time. In the femur 3-point bending test, CHC reduced energy storage, resilience, and ultimate stress, indicating increased susceptibility to micro-damage and fracture, while NHC only decreased energy storage. In the cyclic loading test, both OCs decreased creep indentation distance, but CHC increased the average unloading slope, implying decreased microdamage risk and improved deformation resistance by the OCs. Thus, reduced bone mineralization by the OCs appears to affect bone mechanical properties under static loading, but not its cyclic loading ability. When compared to an age-matched control, after 7 months, CHC affected 24 metabolic pathways in bone and 9 in serum, whereas NHC altered 17 in bone and none in serum. 6 metabolites were common between the serum and bone of CHC rats, suggesting their potential as biomarkers of bone health in women taking CHC.

Conclusion

Both OCs have adverse effects on various skeletal parameters, with CHC having a greater negative impact on bone strength.

Biomechanical evaluation of different strain judging criteria on the prediction precision of cortical bone fracture simulation under compression

Introduction: The principal strain or equivalent strain is mainly used in current numerical studies to determine the mechanical state of the element in the cortical bone finite element model and then perform fracture simulation. However, it is unclear which strain is more suitable for judging the element mechanical state under different loading conditions due to the lack of a general strain judging criterion for simulating the cortical bone fracture. Methods: This study aims to explore a suitable strain judging criterion to perform compressive fracture simulation on the rat femoral cortical bone based on continuum damage mechanics. The mechanical state of the element in the cortical bone finite element model was primarily assessed using the principal strain and equivalent strain separately to carry out fracture simulation. The prediction accuracy was then evaluated by comparing the simulated findings with different strain judging criteria to the corresponding experimental data. Results: The results showed that the fracture parameters predicted using the principal strain were closer to the experimental values than those predicted using the equivalent strain. Discussion: Therefore, the fracture simulation under compression was more accurate when the principal strain was applied to control the damage and failure state in the element. This finding has the potential to improve prediction accuracy in the cortical bone fracture simulation.

Bionic design based on micro-nano structure of osteon and its low-velocity impact damage behavior

It is found that the osteon is composed of thin and thick lamellae which are periodic and approximately concentric, every 5 lamellae is a cycle, the periodic helix angle of mineralized collagen fibers in two adjacent sub-lamellae is 30°. Four bionic composite models with different fiber helix angles were established and fabricated according to the microstructure of mineralized collagen fibers in osteon. Based on the impact analysis of four kinds of bionic composite models, the effects of the fiber periodic spiral structure on the impact resistance and energy dissipation of multi-layer bionic composite were investigated. The analysis results show that the fiber helix angle affects the impact damage resistance and energy dissipation of multi-layer fiber reinforced composites. Among the 4 kinds of multi-layer composite models, the composite model with helix angle of 30° has better comprehensive ability to resist impact damage. The test results show that the impact damage area of the specimen with 30° helix angle is smallest among the 4 types of bionic specimens, which is consistent with the results of finite-element impact analysis. Furthermore, in the case of without impact damage, the smaller the fiber helix angle is, the more uniform the stress distribution is and more energy is dissipated in the impact process. The periodic spiral structure of mineralized collagen fibers in osteon are the result of natural selection of biological evolution. This structure can effectively improve the ability of cortical bone to resist external impact. The research results can provide useful guidance for the design and manufacture of high-performance and strong impact resistant bionic composites.

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.