Failure of Mineralized Collagen Microfibrils Using Finite Element Simulation Coupled to Mechanical Quasi-brittle Damage

Mineralized Collagen Fibrils: An Essential Component in Determining the Mechanical Behavior of Cortical Bone

Bone comprises mechanically different materials in a specific hierarchical structure. Mineralized collagen fibrils (MCFs), represented by tropocollagen molecules and hydroxyapatite nanocrystals, are the fundamental unit of bone. The mechanical characterization of MCFs provides the unique adaptive mechanical competence to bone to withstand mechanical load. The structural and mechanical role of MCFs is critical in the deformation mechanisms of bone and the marvelous strength and toughness possessed by bone. However, the role of MCFs in the mechanical behavior of bone across multiple length scales is not fully understood. In the present study, we shed light upon the latest progress regarding bone deformation at multiple hierarchical levels and emphasize the role of MCFs during bone deformation. We propose the concept of hierarchical deformation of bone to describe the interconnected deformation process across multiple length scales of bone under mechanical loading. Furthermore, how the deterioration of bone caused by aging and diseases impairs the hierarchical deformation process of the cortical bone is discussed. The present work expects to provide insights on the characterization of MCFs in the mechanical properties of bone and lays the framework for the understanding of the multiscale deformation mechanics of bone.

A new biological bone remodeling in silico model combined with advanced discretization methods

Bone remodeling remains a highly researched topic investigated by many strands of science. The main purpose of this work is formulating a new computational framework for biological simulation, extending the version of the bone remodeling model previously proposed by Komarova. Thus, considering only the biological aspect of the remodeling process, the action of osteoclasts and osteoblasts is taken into account as well as its impact on bone mass. It is conducted a spatiotemporal analysis of a remodeling cycle obtaining a dynamic behavior of bone cells very similar to the biological process already described in the literature. The numerical example used is based on bone images obtained with scanning electron microscopy. During simulation, it is possible to observe the variation of bone's architecture through isomaps. These maps are obtained through the combination of biological bone remodeling models with three distinct numerical techniques-finite element method (FEM), radial point interpolation method (RPIM), and natural neighbor radial point interpolation method (NNRPIM). A study combining these numerical techniques allows to compare their performance. Ultimately, this work supports the inclusion of meshless methods due to their smoother results and its easiness to be combined with medical images from CT scans and MRI.

Deproteinization of Cortical Bone: Effects of Different Treatments

Bone is a biological composite material having collagen and mineral as its main constituents. In order to better understand the arrangement of the mineral phase in bone, porcine cortical bone was deproteinized using different chemical treatments. This study aims to determine the best method to remove the protein constituent while preserving the mineral component. Chemicals used were H₂O₂, NaOCl, NaOH, and KOH, and the efficacy of deproteinization treatments was determined by thermogravimetric analysis and Raman spectroscopy. The structure of the residual mineral parts was examined using scanning electron microscopy. X-ray diffraction was used to confirm that the mineral component was not altered by the chemical treatments. NaOCl was found to be the most effective method for deproteinization and the mineral phase was self-standing, supporting the hypothesis that bone is an interpenetrating composite. Thermogravimetric analyses and Raman spectroscopy results showed the preservation of mineral crystallinity and presence of residual organic material after all chemical treatments. A defatting step, which has not previously been used in conjunction with deproteinization to isolate the mineral phase, was also used. Finally, Raman spectroscopy demonstrated that the inclusion of a defatting procedure resulted in the removal of some but not all residual protein in the bone.

CROSS-LINKS MULTISCALE EFFECTS ON BONE ULTRASTRUCTURE BIOMECHANICAL BEHAVIOR

Bone is a multiscale combination of collagen molecules merged with mineral crystals. Its high rigidity and stability stem amply from its polymeric organic matrix and secondly from the connections established between interdifferent and intradifferent scale components through cross-links. Several studies have shown that the cross-links inhibition results in a reduction in strength of bone but they do not quantify the degree to which these connections contribute to the bone rigidity and toughness. This report is classified among the few works that measure the cross-links multiscale impact on the ultrastructure bone mechanical behavior.

This work aims firstly to study the effect of cross-links at the molecule scale and secondly to gather from literature studies results handling with cross-links effects on the other bone ultrastructure scales in order to reveal the multiscale effect of cross-links. This study proves that cross-links increasing number improves the mechanical performance of each scale of bone ultrastructure. On the other hand, cross-links have a multiscale contribution that depends on its rank related to existing cross-links connecting the same geometries and it depends on mechanical characteristics of geometries connected.

A multiscale modelling of bone ultrastructure elastic proprieties using finite elements simulation and neural network method

Bone is a living material with a complex hierarchical structure which entails exceptional mechanical properties, including high fracture toughness, specific stiffness and strength. Bone tissue is essentially composed by two phases distributed in approximately 30-70%: an organic phase (mainly type I collagen and cells) and an inorganic phase (hydroxyapatite-HA-and water). The nanostructure of bone can be represented throughout three scale levels where different repetitive structural units or building blocks are found: at the first level, collagen molecules are arranged in a pentameric structure where mineral crystals grow in specific sites. This primary bone structure constitutes the mineralized collagen microfibril. A structural organization of inter-digitating microfibrils forms the mineralized collagen fibril which represents the second scale level. The third scale level corresponds to the mineralized collagen fibre which is composed by the binding of fibrils. The hierarchical nature of the bone tissue is largely responsible of their significant mechanical properties; consequently, this is a current outstanding research topic. Scarce works in literature correlates the elastic properties in the three scale levels at the bone nanoscale. The main goal of this work is to estimate the elastic properties of the bone tissue in a multiscale approach including a sensitivity analysis of the elastic behaviour at each length scale. This proposal is achieved by means of a novel hybrid multiscale modelling that involves neural network (NN) computations and finite elements method (FEM) analysis. The elastic properties are estimated using a neural network simulation that previously has been trained with the database results of the finite element models. In the results of this work, parametric analysis and averaged elastic constants for each length scale are provided. Likewise, the influence of the elastic constants of the tissue constituents is also depicted. Results highlight that intelligent numerical methods are powerful and accurate procedures to deal with the complex multiscale problem in the bone tissue with results in agreement with values found in literature for specific scale levels.

Modelling of bone fracture and strength at different length scales: a review

In this paper, we review analytical and computational models of bone fracture and strength. Bone fracture is a complex phenomenon due to the composite, inhomogeneous and hierarchical structure of bone. First, we briefly summarize the hierarchical structure of bone, spanning from the nanoscale, sub-microscale, microscale, mesoscale to the macroscale, and discuss experimental observations on failure mechanisms in bone at these scales. Then, we highlight representative analytical and computational models of bone fracture and strength at different length scales and discuss the main findings in the context of experiments. We conclude by summarizing the challenges in modelling of bone fracture and strength and list open topics for scientific exploration. Modelling of bone, accounting for different scales, provides new and needed insights into the fracture and strength of bone, which, in turn, can lead to improved diagnostic tools and treatments of bone diseases such as osteoporosis.

Multiscale approach including microfibril scale to assess elastic constants of cortical bone based on neural network computation and homogenization method

The complexity and heterogeneity of bone tissue require a multiscale modeling to understand its mechanical behavior and its remodeling mechanisms. In this paper, a novel multiscale hierarchical approach including microfibril scale based on hybrid neural network (NN) computation and homogenization equations was developed to link nanoscopic and macroscopic scales to estimate the elastic properties of human cortical bone. The multiscale model is divided into three main phases: (i) in step 0, the elastic constants of collagen-water and mineral-water composites are calculated by averaging the upper and lower Hill bounds; (ii) in step 1, the elastic properties of the collagen microfibril are computed using a trained NN simulation. Finite element calculation is performed at nanoscopic levels to provide a database to train an in-house NN program; and (iii) in steps 2-10 from fibril to continuum cortical bone tissue, homogenization equations are used to perform the computation at the higher scales. The NN outputs (elastic properties of the microfibril) are used as inputs for the homogenization computation to determine the properties of mineralized collagen fibril. The mechanical and geometrical properties of bone constituents (mineral, collagen, and cross-links) as well as the porosity were taken in consideration. This paper aims to predict analytically the effective elastic constants of cortical bone by modeling its elastic response at these different scales, ranging from the nanostructural to mesostructural levels. Our findings of the lowest scale's output were well integrated with the other higher levels and serve as inputs for the next higher scale modeling. Good agreement was obtained between our predicted results and literature data.

Finite element prediction with experimental validation of damage distribution in single trabeculae during three-point bending tests

There is growing evidence that information on trabecular microarchitecture can improve the assessment of fracture risk. One current strategy is to exploit finite element (FE) analysis applied to experimental data of mechanically loaded single trabecular bone tissue obtained from non-invasive imaging techniques for the investigation of the damage initiation and growth of bone tissue. FE analysis of this type of bone has mainly focused on linear and non-linear analysis to evaluate the bone's failure properties. However, there is a lack of experimentally validated FE damage models at trabecular bone tissue level allowing for the simulation of the progressive damage process (initiation and growth) till complete fracture. Such models are needed to perform enhanced prediction of the apparent failure mechanical properties needed to assess the fracture risk of bone organs. In the current study, we develop a FE model based on a continuum damage mechanics (CDM) approach to simulate the damage initiation and propagation of a single trabecula till complete facture in quasi-static regime. Three-point bending experiments were performed on single bovine trabeculae and compared to FE results. In order to validate the proposed FE mode, (i) the force displacement curve was compared to the experimental one and (ii) the damage distribution was correlated to the measured one obtained by digital image correlation based on stress whitening in bone, reported to be correlated to microdamage. A very good agreement was obtained between the FE and experimental results, indicating that the proposed damage investigation protocol based on FE analysis and testing is reliable to assess the damage behavior of bone tissue and that the current damage model is able to accurately simulate the damaging and fracturing process of single trabeculae under quasi static load.

Nanomechanical properties of mineralised collagen microfibrils based on finite elements method: biomechanical role of cross-links

Hierarchical structures in bio-composites such as bone tissue have many scales or levels and synergic interactions between the different levels. They also have a highly complex architecture in order to fulfil their biological and mechanical functions. In this study, a new three-dimensional (3D) model based on the finite elements (FEs) method was used to model the relationship between the hierarchical structure and the properties of the constituents at the sub-structure scale (mineralised collagen microfibrils) and to investigate their apparent nanomechanical properties. The results of the proposed FE simulations show that the elastic properties of microfibrils depend on different factors such as the number of cross-links, the mechanical properties and the volume fraction of phases. The results obtained under compression loading at a small deformation < 2% show that the microfibrils have a Young's modulus (Ef) ranging from 0.4 to 1.16 GPa and a Poisson's ratio ranging from 0.26 to 0.3. These results are in excellent agreement with experimental data (X-ray, AFM and MEMS) and molecular simulations.