Prediction of Cortical Bone Elastic Constants by a Two-Level Micromechanical Model Using a Generalized Self-Consistent Method

Applying a homogeneous pressure distribution to the upper vertebral endplate: Validation of a new loading system, pilot application to human vertebral bodies, and finite element predictions of DIC measured displacements and strains

Image-based personalized Finite Element Models (pFEM) could detect alterations in physiological deformation of human vertebral bodies, but their accuracy has been seldom reported. Meaningful validation experiments should allow vertebral endplate deformability and ensure well-controlled boundary conditions. This study aimed to (i) validate a new loading system to apply a homogeneous pressure on the vertebral endplate during vertebral body compression regardless of endplate deformation; (ii) perform a pilot study on human vertebral bodies measuring surface displacements and strains with Digital Image Correlation (DIC); (iii) determine the accuracy of pFEM of the vertebral bodies. Homogeneous pressure application was achieved by pressurizing a fluid silicone encased in a rubber silicone film acting on the cranial endplate. The loading system was validated by comparing DIC-measured longitudinal strains and lower-end contact pressures, measured on three homogeneous pseudovertebrae of constant transversal section at 2.0 kN, against theoretically calculated values. Longitudinal strains and contact pressures were rather homogeneous, and their mean values close to theoretical calculations (5% underestimation). DIC measurements of surface longitudinal and circumferential displacements and strains were obtained on three human vertebral bodies at 2.0 kN. Complete displacement and strain maps were achieved for anterolateral aspects with random errors ≤0.2 μm and ≤30 μstrain, respectively. Venous plexus and double curvatures limited the completeness and accuracy of DIC data in posterior aspects. pFEM of vertebral bodies, including cortical bone mapping, were built from computed tomography images. In anterolateral aspects, pFEM accuracy of the three vertebrae was: (i) comparable to literature in terms of longitudinal displacements (R2>0.8); (ii) extended to circumferential displacements (pooled data: R2>0.9) and longitudinal strains (zero median error, 95% error: <27%). Circumferential strains were overestimated (median error: 39%). The new methods presented may permit to study how physiological and pathologic conditions influence the ability of vertebral endplates/bodies to sustain loads.

Effect of geometrical structure variations on strength and damage onset of cortical bone using multi-scale cohesive zone based finite element method

Three-dimensional multi-scale finite element models were designed to examine the effects of geometrical structure variations on the damage onset in cortical bone at multiple structural scales. A cohesive zone finite element approach, together with anisotropic damage initiation criteria, is used to predict the onset of damage. The finite element models are developed to account for the onset of microdamage from the microscopic length scales consisting of collagen fibres, to the macroscopic level consisting of osteons and the Haversian canals. Numerical results indicated that the yield strain at the initiation of microcracks is independent of variations in the local mineral volume fraction at each structural scale. Further, the yield strain and strength properties of cortical bone are dependent on its structural anisotropy and hierarchical structure. A positive correlation is observed between bone strength and mineral content at each length scale.

The multiscale meso-mechanics model of viscoelastic cortical bone

Cortical bone is a complex hierarchical structure consisting of biological fiber composites with transversely isotropic constituents, whose microstructures deserve extensive study to understand the mechanism of living organisms and explore development of biomimetic materials. Based on this, we establish a three-level hierarchical structure from microscale to macroscale and propose a multiscale micromechanics model of cortical bone, which considers Haversian canal, osteonal lamellae, cement line and interstitial lamellae. In order to study the microstructural effect on the elastic behavior of hierarchical structures, the Mori-Tanaka model and locally exact homogenization theory are introduced for the homogenization of heterogeneous materials of microstructure at each level. Within sub-microscale, Haversian canal and Osteonal lamella are treated as fiber and matrix, whose homogenization is surrounded with cement line matrix in microstructure (or what we called "osteon") for the second homogenization; finally, osteon and interstitial lamella establish the meso-structure for the third homogenization, predicting the effective moduli of cortical bone. The correctness of the model in this paper is verified against the data in literature with good agreement. Finally, the dynamic viscoelastic response of cortical bones is investigated from a multiscale perspective, where the measured data are substituted into the present models to study the hydration and aging effect on bones' stiffness and viscoelasticity. It is demonstrated that the hydration is much more influential in affecting the storage and loss moduli of cortical bone than the aging effect. We also present a few numerical investigations on microstructural material and geometric parameters on the overall mechanical properties of cortical bone.

Maximum effect of the heterogeneity of tissue mineralization on the effective cortical bone elastic properties

The mineralization level is heterogeneous in cortical bone extracellular matrix as a consequence of remodeling. Models of the effective elastic properties at the millimeter scale have been developed based on idealizations of the vascular pore network and matrix properties. Some popular models do not take into account the heterogeneity of the matrix. However, the errors on the predicted elasticity when the difference in elastic properties between osteonal and interstitial tissues is not modeled have not been quantified. This work provides an estimation of the maximum error. We compare the effective elasticity of a representative volume element (RVE) assuming (1) different elastic properties in osteonal and interstitial tissues vs. (2) average matrix properties. In order to account for the variability of bone microstructure, we use a collection of high resolution images of the pore network to build RVEs. In each RVE we assumed a constant osteonal wall thickness and we artificially varied this thickness between 35 and 140 [Formula: see text]m to create RVEs with different amounts of osteonal tissue. The homogenization problem was solved with a fast Fourier transform (FFT)-based numerical scheme. We found that the error depends on pore volume fraction and varies on average from 1 to [Formula: see text] depending on the assumed diameter of the osteons. The results suggest that matrix heterogeneity may be disregarded in cortical bone models in most practical cases.

Connection between elastic and electrical properties of cortical bone

The paper focuses on the connection between elastic and electrical properties of cortical bone. Both these properties are governed by microstructure that consists of several pore systems filled with mechanically soft and electrically conductive tissue. Microstructural changes induced by aging, various diseases, microgravity conditions etc. lead to variation in both properties. The paper address the problem of evaluation of the changes in mechanical performance (decrease in Young׳s moduli) via monitoring electrical conductivity. The theoretical results are verified experimentally.

Elastic Properties of Human Osteon and Osteonal Lamella Computed by a Bidirectional Micromechanical Model and Validated by Nanoindentation

Knowledge of the anisotropic elastic properties of osteon and osteonal lamellae provides a better understanding of various pathophysiological conditions, such as aging, osteoporosis, osteoarthritis, and other degenerative diseases. For this reason, it is important to investigate and understand the elasticity of cortical bone. We created a bidirectional micromechanical model based on inverse homogenization for predicting the elastic properties of osteon and osteonal lamellae of cortical bone. The shape, the dimensions, and the curvature of osteon and osteonal lamellae are described by appropriately chosen curvilinear coordinate systems, so that the model operates close to the real morphology of these bone components. The model was used to calculate nine orthotropic elastic constants of osteonal lamellae. The input values have the elastic properties of a single osteon. We also expressed the dependence of the elastic properties of the lamellae on the angle of orientation. To validate the model, we performed nanoindentation tests on several osteonal lamellae. We compared the experimental results with the calculated results, and there was good agreement between them. The inverted model was used to calculate the elastic properties of a single osteon, where the input values are the elastic constants of osteonal lamellae. These calculations reveal that the model can be used in both directions of homogenization, i.e., direct homogenization and also inverse homogenization. The model described here can provide either the unknown elastic properties of a single lamella from the known elastic properties at the level of a single osteon, or the unknown elastic properties of a single osteon from the known elastic properties at the level of a single lamella.

To what extent can cortical bone millimeter-scale elasticity be predicted by a two-phase composite model with variable porosity?

An evidence gap exists in fully understanding and reliably modeling the variations in elastic anisotropy that are observed at the millimeter scale in human cortical bone. The porosity (pore volume fraction) is known to account for a large part, but not all, of the elasticity variations. This effect may be modeled by a two-phase micromechanical model consisting of a homogeneous matrix pervaded by cylindrical pores. Although this model has been widely used, it lacks experimental validation. The aim of the present work is to revisit experimental data (elastic coefficients, porosity) previously obtained from 21 cortical bone specimens from the femoral mid-diaphysis of 10 donors and test the validity of the model by proposing a detailed discussion of its hypotheses. This includes investigating to what extent the experimental uncertainties, pore network modeling, and matrix elastic properties influence the model's predictions. The results support the validity of the two-phase model of cortical bone which assumes that the essential source of variations of elastic properties at the millimeter-scale is the volume fraction of vascular porosity. We propose that the bulk of the remaining discrepancies between predicted stiffness coefficients and experimental data (RMSE between 6% and 9%) is in part due to experimental errors and part due to small variations of the extravascular matrix properties. More significantly, although most of the models that have been proposed for cortical bone were based on several homogenization steps and a large number of variable parameters, we show that a model with a single parameter, namely the volume fraction of vascular porosity, is a suitable representation for cortical bone. The results could provide a guide to build specimen-specific cortical bone models. This will be of interest to analyze the structure-function relationship in bone and to design bone-mimicking materials.

Micromechanical modeling of elastic properties of cortical bone accounting for anisotropy of dense tissue

The paper analyzes the connection between microstructure of the osteonal cortical bone and its overall elastic properties. The existing models either neglect anisotropy of the dense tissue or simplify cortical bone microstructure (accounting for Haversian canals only). These simplifications (related mostly to insufficient mathematical apparatus) complicate quantitative analysis of the effect of microstructural changes - produced by age, microgravity, or some diseases - on the overall mechanical performance of cortical bone. The present analysis fills this gap; it accounts for anisotropy of the dense tissue and uses realistic model of the porous microstructure. The approach is based on recent results of Sevostianov et al. (2005) and Saadat et al. (2012) on inhomogeneities in a transversely-isotropic material. Bone's microstructure is modeled according to books of Martin and Burr (1989), Currey (2002), and Fung (1993) and includes four main families of pores. The calculated elastic constants for porous cortical bone are in agreement with available experimental data. The influence of each of the pore types on the overall moduli is examined.

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.

Biomechanical determinants of the stability of dental implants: influence of the bone-implant interface properties

Dental implants are now widely used for the replacement of missing teeth in fully or partially edentulous patients and for cranial reconstructions. However, risks of failure, which may have dramatic consequences, are still experienced and remain difficult to anticipate. The stability of biomaterials inserted in bone tissue depends on multiscale phenomena of biomechanical (bone-implant interlocking) and of biological (mechanotransduction) natures. The objective of this review is to provide an overview of the biomechanical behavior of the bone-dental implant interface as a function of its environment by considering in silico, ex vivo and in vivo studies including animal models as well as clinical studies. The biomechanical determinants of osseointegration phenomena are related to bone remodeling in the vicinity of the implants (adaptation of the bone structure to accommodate the presence of a biomaterial). Aspects related to the description of the interface and to its space-time multiscale nature will first be reviewed. Then, the various approaches used in the literature to measure implant stability and the bone-implant interface properties in vitro and in vivo will be described. Quantitative ultrasound methods are promising because they are cheap, non invasive and because of their lower spatial resolution around the implant compared to other biomechanical approaches.

Elastic moduli of untreated, demineralized and deproteinized cortical bone: validation of a theoretical model of bone as an interpenetrating composite material

A theoretical experimentally based multi-scale model of the elastic response of cortical bone is presented. It portrays the hierarchical structure of bone as a composite with interpenetrating biopolymers (collagen and non-collagenous proteins) and minerals (hydroxyapatite), together with void spaces (porosity). The model involves a bottom-up approach and employs micromechanics and classical lamination theories of composite materials. Experiments on cortical bone samples from bovine femur include completely demineralized and deproteinized bones as well as untreated bone samples. Porosity and microstructure are characterized using optical and scanning electron microscopy, and micro-computed tomography. Compression testing is used to measure longitudinal and transverse elastic moduli of all three bone types. The characterization of structure and properties of these three bone states provides a deeper understanding of the contributions of the individual components of bone to its elastic response and allows fine tuning of modeling assumptions. Very good agreement is found between theoretical modeling and compression testing results, confirming the validity of the interpretation of bone as an interpenetrating composite material.

Specimen-specific multi-scale model for the anisotropic elastic constants of human cortical bone

The anisotropic elastic constants of human cortical bone were predicted using a specimen-specific micromechanical model that accounted for structural parameters across multiple length scales. At the nano-scale, the elastic constants of the mineralized collagen fibril were estimated from measured volume fractions of the constituent phases, namely apatite crystals and Type I collagen. The elastic constants of the extracellular matrix (ECM) were predicted using the measured orientation distribution function (ODF) for the apatite crystals to average the contribution of misoriented mineralized collagen fibrils. Finally, the elastic constants of cortical bone tissue were determined by accounting for the measured volume fraction of Haversian porosity within the ECM. Model predictions using the measured apatite crystal ODF were not statistically different from experimental measurements for both the magnitude and anisotropy of elastic constants. In contrast, model predictions using common idealized assumptions of perfectly aligned or randomly oriented apatite crystals were significantly different from the experimental measurements. A sensitivity analysis indicated that the apatite crystal volume fraction and ODF were the most influential structural parameters affecting model predictions of the magnitude and anisotropy, respectively, of elastic constants.

Effect of porosity on effective diagonal stiffness coefficients (cii) and elastic anisotropy of cortical bone at 1 MHz: a finite-difference time domain study

Finite-difference time domain (FDTD) numerical simulations coupled to real experimental data were used to investigate the propagation of 1 MHz pure bulk wave propagation through models of cortical bone microstructures. Bone microstructures were reconstructed from three-dimensional high resolution synchrotron radiation microcomputed tomography (SR-muCT) data sets. Because the bone matrix elastic properties were incompletely documented, several assumptions were made. Four built-in bone matrix models characterized by four different anisotropy ratios but the same Poisson's ratios were tested. Combining them with the reconstructed microstructures in the FDTD computations, effective stiffness coefficients were derived from simulated bulk-wave velocity measurements. For all the models, all the effective compression and shear bulk wave velocities were found to decrease when porosity increases. However, the trend was weaker in the axial direction compared to the transverse directions, contributing to the increase of the effective anisotropy. On the other hand, it was shown that the initial Poisson's ratio value may substantially affect the variations of the effective stiffness coefficients. The present study can be used to elaborate sophisticated macroscopic computational bone models incorporating realistic CT-based macroscopic bone structures and effective elastic properties derived from muCT-based FDTD simulations including the cortical porosity effect.