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.     

Micromechanics fracture in osteonal cortical bone: a study of the interactions between microcrack propagation, microstructure and the material properties

A two-dimensional micromechanical fibre reinforced composite materials model for osteonal cortical bone is presented. The interstitial bone is modelled as a matrix, the osteons are modelled as fibres, and the cement line is presented as interface tissue. The interaction between osteons and microcracks is evaluated by linear elastic fracture mechanics theory, followed by a determination of the stress intensity factor at the vicinity of the microcrack tips. The results indicate that bone microstructural heterogeneity greatly influences fracture parameters. Furthermore, microstructural morphology and loading conditions affect growth trajectories, the microcrack propagation trajectory deviates from the osteon under tensile loading, and osteon penetration is observed under compressive loads.     

The dependence of transversely isotropic elasticity of human femoral cortical bone on porosity

The objective of this study was to examine the dependence of the elastic properties of cortical bone as a transversely isotropic material on its porosity. The longitudinal Young's modulus, transverse Young's modulus, longitudinal shear modulus, transverse shear modulus, and longitudinal Poisson's ratio of cortical bone were determined from eighteen groups of longitudinal and transverse specimens using tensile and torsional tests on a servo-hydraulic material testing system. These cylindrical waisted specimens of cortical bone were harvested from the middle diaphysis of three pairs of human femora. The porosity of these specimens was assessed by means of histology. Our study demonstrated that the longitudinal Young's and shear moduli of human femoral cortical bone were significantly (p<0.01) negatively correlated with the porosity of cortical bone. Conversely, the elastic properties in the transverse direction did not have statistically significant correlations with the porosity of cortical bone. As a result, the transverse elastic properties of cortical bone were less sensitive to changes in porosity than those in the longitudinal direction. Additionally, the anisotropic ratios of cortical bone elasticity were found to be significantly (p<0.01) negatively correlated with its porosity, indicating that cortical bone tended to become more isotropic when its porosity increased. These results may help a number of researchers develop more accurate micromechanics models of cortical bone.     

Micromechanics modeling of Haversian cortical bone properties.

A finite-element micromechanics model for Haversian cortical bone tissue has been developed and studied. The model is an extension of two-dimensional micromechanics techniques for fiber-reinforced composite materials. Haversian systems, or secondary osteons, are considered to be the fiber component, and interstitial lamellar bone the matrix material. The cement line is included as an 'interphase' component along the fiber/matrix interface. The model assumes a regular repeatable spacing of the longitudinally aligned continuous fibers and is, therefore, restricted to approximating Haversian cortical bone in its present form. Haversian porosity is modeled explicitly by incorporating a hollow fiber to represent the Haversian canal. Solutions have been obtained by applying uniform macroscopic stresses to the boundaries of the repeating unit cell model. Macroscopic mechanical property predictions correspond reasonably well with the experimental data for cortical bone, but are necessarily dependent on the input properties for each constituent, which are not well established. The predicted variation in the elastic modulus with porosity is not as sensitive as that observed experimentally. Stresses within the constituents can also be modeled with this method and are demonstrated to deviate from the macroscopic applied stress levels.     

A determination of the minimum sizes of representative volume elements for the prediction of cortical bone elastic properties

At its highest level of microstructural organization—the mesoscale or millimeter scale—cortical bone exhibits a heterogeneous distribution of pores (Haversian canals, resorption cavities). Multi-scale mechanical models rely on the definition of a representative volume element (RVE). Analytical homogenization techniques are usually based on an idealized RVE microstructure, while finite element homogenization using high-resolution images is based on a realistic RVE of finite size. The objective of this paper was to quantify the size and content of possible cortical bone mesoscale RVEs. RVE size was defined as the minimum size: (1) for which the apparent (homogenized) stiffness tensor becomes independent of the applied boundary conditions or (2) for which the variance of elastic properties for a set of microstructure realizations is sufficiently small. The field of elastic coefficients and microstructure in RVEs was derived from one acoustic microscopy image of a human femur cortical bone sample with an overall porosity of 8.5%. The homogenized properties of RVEs were computed with a finite element technique. It was found that the size of the RVE representative of the overall tissue is about 1.5 mm. Smaller RVEs (~0.5 mm) can also be considered to estimate local mesoscopic properties that strongly depend on the local pores volume fraction. This result provides a sound basis for the application of homogenization techniques to model the heterogeneity of cortical microstructures. An application of the findings to estimate elastic properties in the case of a porosity gradient is briefly presented.     

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.     

The relative effects of collagen fiber orientation, porosity, density, and mineralization on bone strength

This investigation determined the relative importance of collagen fiber orientation, porosity, density, and mineralization in determining the tensile strength of bovine cortical bone. Thirty-nine specimens were tested for failure stress and the values of eight histologic and compositional parameters: collagen fiber orientation, wet and dry apparent density, percent mineralization of the bone matrix, and several components of porosity (Haversian canals, Volkmann's canals, and plexiform vascular spaces). Linear regression analysis showed that collagen fiber orientation was consistently the single best predictor of strength. Mineralization of the bone matrix was generally a poor predictor of strength. Density and porosity ranked between these variables in importance. Multiple regression equations containing all significantly correlated variables achieved correlation coefficients of 0.607 for plexiform bone and 0.881 for osteonal bone. Also, separate analysis of plexiform and osteonal specimens showed that the latter group was weaker even though it was less porous, apparently because it had collagen fibers which were less longitudinally oriented. This study suggests it is feasible to develop better empirical formulae for the prediction of cortical bone strength than are currently available if a variety of variables is introduced. Additional data are needed to confirm these results.     

Micromechanical modelling of cortical bone

Cortical bone is a heterogeneous material with a complex hierarchical microstructure. In this work, unit cell finite element models were developed to investigate the effect of microstructural morphology on the macroscopic properties of cortical bone. The effect of lacunar and vascular porosities, percentage of osteonal bone and orientation of the Haversian system on the macroscopic elastic moduli and Poisson's ratios was investigated. The results presented provide relationships for applying more locally accurate material properties to larger scale and whole bone models of varying porosity. Analysis of the effect of the orientation of the Haversian system showed that its effects should not be neglected in larger scale models. This study also provides insight into how microstructural features effect local distributions and cause a strain magnification effect. Limitations in applying the unit cell methodology approach to bone are also discussed.     

Fluid pressure relaxation depends upon osteonal microstructure: modeling an oscillatory bending experiment.

When bone is mechanically loaded, bone fluid flow induces shear stresses on bone cells that have been proposed to be involved in bone's mechanosensory system. To investigate bone fluid flow and strain-generated potentials, several theoretical models have been proposed to mimic oscillatory four-point bending experiments performed on thin bone specimens. While these previous models assume that the bone fluid relaxes across the specimen thickness, we hypothesize that the bone fluid relaxes primarily through the vascular porosity (osteonal canals) instead and develop a new poroelastic model that integrates the microstructural details of the lacunar-canalicular porosity, osteonal canals, and the osteonal cement lines. Local fluid pressure profiles are obtained from the model, and we find two different fluid relaxation behaviors in the bone specimen, depending on its microstructure: one associated with the connected osteonal canal system, through which bone fluid relaxes to the nearby osteonal canals; and one associated with the thickness of a homogeneous porous bone specimen (approximately 1 mm in our model), through which bone fluid relaxes between the external surfaces of the bone specimen at relatively lower loading frequencies. Our results suggest that in osteonal bone specimens the fluid pressure response to cyclic loading is not sensitive to the permeability of the osteonal cement lines, while it is sensitive to the applied loading frequency. Our results also reveal that the fluid pressure gradients near the osteonal canals (and thus the fluid shear stresses acting on the nearby osteocytes) are significantly amplified at higher loading frequencies.     

Simulation of creep in non-homogenous samples of human cortical bone

Characterising the mechanisms causing viscoelastic mechanical properties of human cortical bone, as well as understanding sources of variation, is important in predicting response of the bone to creep and fatigue loads. Any better understanding, when incorporated into simulations including finite element analysis, would assist bioengineers, clinicians and biomedical scientists. In this study, we used an empirically verified model of creep strain accumulation, in a simulation of 10 non-homogeneous samples, which were created from micro-CT scans of human cortical bone of the femur midshaft obtained from a 74-year-old female cadaver. These non-homogeneous samples incorporate the presence of Haversian canals and resorption cavities. The influence of inhomogeneity on the response and variation in the samples in both creep and stress relaxation tests are examined. The relationship between steady-state creep rate, applied loads (stress relaxation and creep tests) and microstructure, that is bone apparent porosity, is obtained. These relations may provide insight into damage accumulation of whole human bones and be relevant to studies on osteoporosis.     

Simulation of creep in non-homogenous samples of human cortical bone

Characterising the mechanisms causing viscoelastic mechanical properties of human cortical bone, as well as understanding sources of variation, is important in predicting response of the bone to creep and fatigue loads. Any better understanding, when incorporated into simulations including finite element analysis, would assist bioengineers, clinicians and biomedical scientists. In this study, we used an empirically verified model of creep strain accumulation, in a simulation of 10 non-homogeneous samples, which were created from micro-CT scans of human cortical bone of the femur midshaft obtained from a 74-year-old female cadaver. These non-homogeneous samples incorporate the presence of Haversian canals and resorption cavities. The influence of inhomogeneity on the response and variation in the samples in both creep and stress relaxation tests are examined. The relationship between steady-state creep rate, applied loads (stress relaxation and creep tests) and microstructure, that is bone apparent porosity, is obtained. These relations may provide insight into damage accumulation of whole human bones and be relevant to studies on osteoporosis.     

Elastic modulus and hardness of cortical and trabecular bone lamellae measured by nanoindentation in the human femur.

The mechanical properties of bone tissue are determined by composition as well as structural, microstructural and nanostructural organization. The aim of this study was to quantify the elastic properties of bone at the lamellar level and compare these properties among osteonal, interstitial and trabecular microstructures from the diaphysis and the neck of the human femur. A nanoindentation technique with a custom irrigation system was used for simultaneously measuring force and displacement of a diamond tip pressed 500 nm into the moist bone tissue. An isotropic elastic modulus was calculated from the unloading curve with an assumed Poisson ratio of 0.3, while hardness was defined as the maximal force divided by the corresponding contact area. The elastic moduli ranged from 6.9 +/- 4.3 GPa in trabecular tissue from the femoral neck of a 74 yr old female up to 25.0 +/- 4.3 GPa in interstitial tissue from the diaphyseal cortex of a 69 yr old female. The mean elastic modulus was found to be significantly influenced by the type of lamella (p < 10(-6)) and by donor (p < 10(-6)). The interaction between the type of lamella and the donor was also highly significant (p < 10(-6)). Hardness followed a similar distribution as elastic modulus among types of lamellae and donor, but with lower statistical contrast. It is concluded that the nanostructure of bone tissue must differ substantially among lamellar types, anatomical sites and individuals and suggests that tissue heterogeneity is of potential importance in bone fragility and adaptation.     

Estimation of the effective transversely isotropic elastic constants of a material from known values of the material's orthotropic elastic constants

A method is illustrated for determining the effective transversely isotropic (or isotropic) elastic constants from measured orthotropic elastic constants. This method consists of constructing upper and lower bounds on the effective transversely isotropic (or isotropic) elastic constants using the known orthotropic values. This method is illustrated using three sets of elastic constants for bone. Fortunately, the upper and lower bounds are very close. Thus very good approximations for the effective transversely isotropic (or isotropic) elastic constants for cortical and cancellous bone are obtained from previously published data on the orthotropic elastic constants for those tissue types. This work is undertaken to build a greater database for the transversely isotropic elastic constants of bone with the intention of employing them in a transversely isotropic model of bone poroelasticity. An interesting aspect of the present result is that the Voigt and Reuss bounds are very tight for these anisotropic materials. This is not always the case for these bounds. Received: 14 November 2001 / Accepted: 25 February 2002     

Elastic modulus variation in mandibular bone: a microindentation study of Macaca fascicularis

We characterized the heterogeneous anisotropic elastic properties of mandibular bone in an adult female specimen of Macaca fascicularis using the technique of microindentation. This approach used an indenter of known mass and geometry to sample bone hardness at a spatial resolution in the order of 100 mum. Hardness values were converted to elastic modulus using empirically derived regression. We determined properties in alveolar, midcorpus, and basal regions of coronal and transverse sections taken from multiple locations along the corpus and ramus. Within sections, we determined properties from endosteal, midcortical, and periosteal regions. We found regional variations in bone structure, including bands of orthotropic circumferential lamellar bone at the endosteal and periosteal corpus base, angular region, and ramus. Transversely isotropic osteonal bone characterizes the midcortices of alveolar and basal regions, with many resorption spaces in alveolar regions restricting sampling opportunities. Regional variations in elasticity include relatively compliant bone in the anterior corpus and ramus. Basal cortical bone is stiffer longitudinally than transversely or superoinferiorly, while the evidence for directional dependence in alveolar bone is equivocal. Alveolar bone appears to be relatively compliant with respect to bone found in midcorpus or basal regions. Considerable variation exists in structure and material properties on a highly localized scale, more so than is discernible through conventional approaches for determining material property variation.     

En bloc staining of articular cartilage and bone

Blocks of canine and porcine articular cartilage were stained en bloc with Weigert's iron hematoxylin or Harris' hematoxylin with or without eosin Y counterstaining and cleared in methyl salicylate. The morphology and three-dimensional relationships of chondrocytes were best demonstrated with Weigert's iron hematoxylin. The morphology of the cartilage and chondrocytes was superior to that in sections of routine hematoxylin and eosin stained, paraffin processed samples. The three-dimensional localization of intracellular lipids in individual and clones of chondrocytes was observed when cartilage samples were stained with oil red O and mounted directly in a water-based medium. Blocks of decalcified bone were stained en bloc with Weigert's iron hematoxylin and cleared with methyl salicylate. The three-dimensional orientation of osteocytes around osteonal canals, in circumferential lamellae, and in interstitial lamellae was demonstrated. The morphology of "cutting cones" in cortical bone also was observed.     

En Bloc Staining of Articular Cartilage and Bone

Blocks of canine and porcine articular cartilage were stained en bloc with Weigert's iron hematoxylin or Harris' hematoxylin with or without eosin Y counterstaining and cleared in methyl salicylate. The morphology and three-dimensional relationships of chondrocytes were best demonstrated with Weigert's iron hematoxylin. The morphology of the cartilage and chondrocytes was superior to that in sections of routine hematoxylin and eosin stained, paraffin processed samples. The three-dimensional localization of intracellular lipids in individual and clones of chondrocytes was observed when cartilage samples were stained with oil red O and mounted directly in a water-based medium. Blocks of decalcified bone were stained en bloc with Weigert's iron hematoxylin and cleared with methyl salicylate. The three-dimensional orientation of osteocytes around osteonal canals, in circumferential lamellae, and in interstitial lamellae was demonstrated. The morphology of “cutting cones” in cortical bone also was observed.     

Can the diverse elastic properties of trabecular and cortical bone be attributed to only a few tissue-independent phase properties and their interactions?

As candidates for tissue-independent phases of cortical and trabecular bone we consider (i) hydroxyapatite, (ii) collagen, (iii) ultrastructural water and non-collagenous organic matter, and (iv) marrow (water) filling the Haversian canals and the intertrabecular space. From experiments reported in the literature, we assign stiffness properties to these phases (experimental set I).On the basis of these phase definitions, we develop, within the framework of continuum micromechanics, a two-step homogenization procedure: (i) at a length scale of 100–200 nm, hydroxyapatite (HA) crystals build up a crystal foam (polycrystal), and water and non-collagenous organic matter fill the intercrystalline space (homogenization step I); (ii) at the ultrastructural scale of mineralized tissues (5–10 microns), collagen assemblies composed of collagen molecules are embedded into the crystal foam, acting mechanically as cylindrical templates. At an enlarged material scale of 5–10 mm, the second homogenization step also accommodates the micropore space as cylindrical pore inclusions (Haversian and Volkmann canals, inter-trabecular space) that are suitable for both trabecular and cortical bone. The inputs for this micromechanical model are the tissue-specific volume fractions of HA, collagen, and of the micropore space. The outputs are the tissue-specific ultrastructural and microstructural (=macroscopic=apparent) elasticity tensors.A second independent experimental set (composition data and experimental stiffness values) is employed to validate the proposed model. We report a small mean prediction error for the macroscopic stiffness values. The validation suggests that hydroxyapatite, collagen, and water are tissue-independent phases, which define, through their mechanical interaction, the elasticity of all bones, whether cortical or trabecular.     

Mechanical properties of porcine femoral cortical bone measured by nanoindentation

This study uses a nanoindentation technique to examine variations in the local mechanical properties of porcine femoral cortical bone under hydrated conditions. Bone specimens from three age groups (6, 12 and 42 months), representing developing bone, ranging from young to mature animals, were tested on the longitudinal and transverse cross-sectional surfaces. Elastic modulus and hardness of individual lamellae within bone's microstructure: laminar bone, interstitial bone, and osteons, were measured. Both the elastic modulus and hardness increased with age. However, the magnitudes of these increases were different for each microstructural component. The longitudinal moduli were higher than the transverse moduli. Dehydrated samples were also tested to allow a comparison with hydrated samples and these resulted in higher moduli and hardness than the hydrated samples. Again, the degree of variation was different for each microstructural component. These results indicate that the developmental changes in bone have different rates of mechanical change within each microstructural component.     

Patterns of femoral bone remodeling dynamics in a medieval Nubian population

The relationship between age, sex and histomorphometry in femoral cortical bone was examined in a skeletal population of late Medieval antiquity (AD 1250–1450) from Kulubnarti, in Sudanese Nubia. These skeletal remains are naturally mummified and in an excellent state of preservation. The study sample consisted of femoral cross sections from 24 females and 19 males ranging in age from 20 to 50+ years. Femoral cross sections were examined using an image analysis system. Numbers of secondary osteons and osteon fragments were counted, osteon area and Haversian canal area were measured, and several variables were calculated to assess differences between sexes and among age groups in bone remodeling variables. The results indicate significant differences between the sexes in osteon number and size. Males had significantly more intact osteons than females, whereas females had significantly larger osteons than males. Haversian canal dimensions were not statistically significant between the sexes. Sex differences in activity patterns in which males were involved in more physically strenuous tasks may have contributed to differences in remodeling variables. Interpopulational comparisons suggest that mechanical strain affects the microstructural features examined in this study. In particular, small Haversian canals in some archaeological skeletal populations are associated with higher bone volume, which may result from high levels of mechanical strain. Am J Phys Anthropol 104:133–146, 1997. © 1997 Wiley-Liss, Inc.     

A two-parameter model of the effective elastic tensor for cortical bone

Multiscale models of cortical bone elasticity require a large number of parameters to describe the organization and composition of the tissue. We hypothesize that the macro-scale anisotropic elastic properties of different bones can be modeled retaining only two variable parameters, and setting the others to universal values identical for all bones. Cortical bone is regarded as a two-phase composite material: a dense mineralized matrix (ultrastructure) and a soft phase (pores). The ultrastructure is assumed to be a homogeneous and transversely isotropic tissue whose elastic properties in different directions are mutually dependent and can be scaled with a single parameter driving the overall rigidity. This parameter is taken to be the volume fraction of mineral f(ha). The pore network is modeled as an ensemble of water-filled cylinders and described only by the porosity p. The effective macroscopic elasticity tensor C(ij)(f(ha),p) is calculated with a multiscale micromechanics approach starting from existing models. The modeled stiffness coefficients compare favorably to four literature datasets which were chosen because they provide the full stiffness tensors of groups of human samples. Since the physical counterparts of f(ha) and p were unknown for the datasets, their values which allow the best fit of experimental tensors by the modeled ones were determined by optimization. Optimum values of f(ha) and p are found to be unique and realistic. These results suggest that a two-parameter model may be sufficient to model the elasticity of different samples of human femora and tibiae. Such a model would in particular be useful in large-scale parametric studies of bone mechanical response.