The elastic properties of trabecular and cortical bone tissues are similar: results from two microscopic measurement techniques

Acoustic microscopy (30-60 microm resolution) and nanoindentation (1-5 microm resolution) are techniques that can be used to evaluate the elastic properties of human bone at a microstructural level. The goals of the current study were (1) to measure and compare the Young's moduli of trabecular and cortical bone tissues from a common human donor, and (2) to compare the Young's moduli of bone tissue measured using acoustic microscopy to those measured using nanoindentation. The Young's modulus of cortical bone in the longitudinal direction was about 40% greater than (p<0.01) the Young's modulus in the transverse direction. The Young's modulus of trabecular bone tissue was slightly higher than the transverse Young's modulus of cortical bone, but substantially lower than the longitudinal Young's modulus of cortical bone. These findings were consistent for both measurement methods and suggest that elasticity of trabecular tissue is within the range of that of cortical bone tissue. The calculation of Young's modulus using nanoindentation assumes that the material is elastically isotropic. The current results, i.e., the average anisotropy ratio (E(L)/E(T)) for cortical bone determined by nanoindentation was similar to that determined by the acoustic microscope, suggest that this assumption does not limit nanoindentation as a technique for measurement of Young's modulus in anisotropic bone.     

Errors induced by off-axis measurement of the elastic properties of bone

Misalignment between the axes of measurement and the material symmetry axes of bone causes error in anisotropic elastic property measurements. Measurements of Poisson's ratio were strongly affected by misalignment errors. The mean errors in the measured Young's moduli were 9.5 and 1.3 percent for cancellous and cortical bone, respectively, at a misalignment angle of 10 degrees. Mean errors of 1.1 and 5.0 percent in the measured shear moduli for cancellous and cortical bone, respectively, were found at a misalignment angle of 10 degrees. Although, cancellous bone tissue was assumed to have orthotropic elastic symmetry, the possibility of the greater symmetry of transverse isotropy was investigated. When the nine orthotropic elastic constants were forced to approximate the five transverse isotropic elastic constants, errors of over 60 percent were introduced. Therefore, it was concluded that cancellous bone is truly orthotropic and not transversely isotropic. A similar but less strong result for cortical bone tissue was obtained.     

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.     

The role of an effective isotropic tissue modulus in the elastic properties of cancellous bone.

Conceptually, the elastic characteristics of cancellous bone could be predicted directly from the trabecular morphology--or architecture--and by the elastic properties of the tissue itself. Although hardly any experimental evidence exists, it is often implicitly assumed that tissue anisotropy has a negligible effect on the apparent elastic properties of cancellous bone. The question addressed in this paper is whether this is actually true. If it is, then micromechanical finite element analysis (micro-FEA) models, representing trabecular architecture, using an 'effective isotropic tissue modulus' should be able to predict apparent elastic properties of cancellous bone. To test this, accurate multi-axial compressive mechanical tests of 29 whale bone specimens were simulated with specimen-specific micro-FEA computer models built from true three-dimensional reconstructions. By scaling the micro-FEA predictions by a constant tissue modulus, 92% of the variation of Young's moduli determined experimentally could be explained. The correlation even increased to 95% when the micro-FEA moduli were scaled to the isotropic tissue moduli of individual specimens. Excellent agreement was also found in the elastic symmetry axes and anisotropy ratios. The prediction of Poisson's ratios was somewhat less precise at 85% correlation. The results support the hypothesis; for practical purposes, the concept of an 'effective isotropic tissue modulus' concept is a viable one. They also suggest that the value of such a modulus for individual cases might be inferred from the average tissue density, hence the degree of mineralization. Future studies must clarify how specific the tissue modulus should be for different types of bone if adequate predictions of elastic behavior are to be made in this way.     

On the dependence of the elasticity and strength of cancellous bone on apparent density

This paper presents a statistical analysis of the pooled data from a number of previous experiments concerning the dependence of the Young's moduli and strength of cancellous bone tissue upon apparent density. The results show that both the Young's moduli and the strength are proportional to the square of apparent density of the tissue and are therefore proportional to one another. It is shown that the coefficient of proportionality is different for human and bovine tissue. It is concluded that the suggestion of Wolff (Das Gesetz der Transformation der Knochen, Hirschwald, Berlin, 1892) that compact bone tissue is simply more dense cancellous bone tissue is not an accurate statement when only the mechanical properties of these two tissues are considered. It is noted that estimates for the elastic modulus of the individual trabecula of human cancellous bone vary from 1 to 20 GPa and it is suggested that this question needs further study.     

Comparison of the elastic and yield properties of human femoral trabecular and cortical bone tissue

The ability to determine trabecular bone tissue elastic and failure properties has biological and clinical importance. To date, trabecular tissue yield strains remain unknown due to experimental difficulties, and elastic moduli studies have reported controversial results. We hypothesized that the elastic and tensile and compressive yield properties of trabecular tissue are similar to those of cortical tissue. Effective tissue modulus and yield strains were calibrated for cadaveric human femoral neck specimens taken from 11 donors, using a combination of apparent-level mechanical testing and specimen-specific, high-resolution, nonlinear finite element modeling. The trabecular tissue properties were then compared to measured elastic modulus and tensile yield strain of human femoral diaphyseal cortical bone specimens obtained from a similar cohort of 34 donors. Cortical tissue properties were obtained by statistically eliminating the effects of vascular porosity. Results indicated that mean elastic modulus was 10% lower (p<0.05) for the trabecular tissue (18.0+/-2.8 GPa) than for the cortical tissue (19.9+/-1.8 GPa), and the 0.2% offset tensile yield strain was 15% lower for the trabecular tissue (0.62+/-0.04% vs. 0.73+/-0.05%, p<0.001). The tensile-compressive yield strength asymmetry for the trabecular tissue, 0.62 on average, was similar to values reported in the literature for cortical bone. We conclude that while the elastic modulus and yield strains for trabecular tissue are just slightly lower than those of cortical tissue, because of the cumulative effect of these differences, tissue strength is about 25% greater for cortical bone.     

Elastic constants of composites formed from PMMA bone cement and anisotropic bovine tibial cancellous bone

An ultrasonic pulse-transit time technique is used to determine the nine orthotropic engineering constants of 32 cement-cancellous bone composites as a function of volume fractions of bone ranging from 0.0 to 0.4. The composites are manufactured using well-aligned bovine cancellous bone from the proximal end of the tibia and low viscosity bone cement. Selected composites are also subjected to mechanical compression tests to compare with the ultrasonic results. There is excellent correlation between the dynamic or ultrasonically determined moduli and the static or mechanically determined moduli; the dynamic moduli are approximately twice the static moduli and this difference is thought to be due to the effect of strain rate. An orthotropic model is assumed requiring nine independent elastic constants to be determined. The dynamic Young's modulus in the direction of major trabecular alignment, E1, increases linearly from 4.9 to 10.4 GPa as bone volume fraction increases from 0 to 0.4; dynamic E2 and E3 values increase from 4.9 to 7 GPa as bone volume fractions increase from 0 to 0.4, with E2 being slightly higher than E3. The dynamic shear modulus, G12, increases from 1.8 to 3.0 GPa, and G31 and G23 increase slightly from 1.8 to 2.2 GPa as bone volume fractions increase from 0 to 0.4. The Poisson's ratios are more sensitive than the Young's moduli and shear moduli to experimental error in the velocity measurements. The mechanically tested modulus (static modulus) in the direction of major trabecular alignment, E1, increases with volume fraction of bone from 2.4 to 4.4 GPa as the bone volume fraction increases from 0 to 0.25; static E2 and E3 values are either equal to or lower than that of pure PMMA.     

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.     

An application of nanoindentation technique to measure bone tissue Lamellae properties

Measuring the microscopic mechanical properties of bone tissue is important in support of understanding the etiology and pathogenesis of many bone diseases. Knowledge about these properties provides a context for estimating the local mechanical environment of bone related cells thait coordinate the adaptation to loads experienced at the whole organ level. The objective of this study was to determine the effects of experimental testing parameters on nanoindentation measures of lamellar-level bone mechanical properties. Specifically, we examined the effect of specimen preparation condition, indentation depth, repetitive loading, time delay, and displacement rate. The nanoindentation experiments produced measures of lamellar elastic moduli for human cortical bone (average value of 17.7 +/- 4.0 GPa for osteons and 19.3 +/- 4.7 GPa for interstitial bone tissue). In addition, the hardness measurements produced results consistent with data in the literature (average 0.52 +/- 0.15 GPa for osteons and 0.59 +/- 0.20 GPa for interstitial bone tissue). Consistent modulus values can be obtained from a 500-nm-deep indent. The results also indicated that the moduli and hardnesses of the dry specimens are significantly greater (22.6% and 56.9%, respectively) than those of the wet and wet and embedded specimens. The latter two groups were not different. The moduli obtained at a 5-nm/s loading rate were significantly lower than the values at the 10- and 20-nm/s loading rates while the 10- and 20-nm/s rates were not significantly different. The hardness measurements showed similar rate-dependent results. The preliminary results indicated that interstitial bone tissue has significantly higher modulus and hardness than osteonal bone tissue. In addition, a significant correlation between hardness and elastic modulus was observed.     

Reliable simulations of the human proximal femur by high-order finite element analysis validated by experimental observations

Background: The mechanical response of patient-specific bone to various load conditions is of major clinical importance in orthopedics. Herein we enhance the methods presented in Yosibash et al. [2007. A CT-based high-order finite element analysis of the human proximal femur compared to in-vitro experiments. ASME Journal of Biomechanical Engineering 129(3), 297–309.] for the reliable simulations of the human proximal femur by high-order finite elements (FEs) and validate the simulations by experimental observations.

Method of approach: A fresh-frozen human femur was scanned by quantitative computed tomography (QCT) and thereafter loaded (in vitro experiments) by a quasi-static force of up to 1250 N. QCT scans were manipulated to generate a high-order FE bone model with distinct cortical and trabecular regions having inhomogeneous isotropic elastic properties with Young's modulus represented by continuous spatial functions. Sensitivity analyses were performed to quantify parameters that mostly influence the mechanical response. FE results were compared to displacements and strains measured in the experiments.

Results: Young moduli correlated to QCT Hounsfield Units by relations in Keyak and Falkinstein [2003. Comparison of in situ and in vitro CT scan-based finite element model predictions of proximal femoral fracture load. Medical Engineering and Physics 25, 781–787.] were found to provide predictions that match the experimental results closely. Excellent agreement was found for both the displacements and strains. The presented study demonstrates that reliable and validated high-order patient-specific FE simulations of human femurs based on QCT data are achievable for clinical computer-aided decision making.     

How is the indentation modulus of bone tissue related to its macroscopic elastic response? A validation study

This work consists of the validation of a novel approach to estimate the local anisotropic elastic constants of the bone extracellular matrix using nanoindentation. For this purpose, nanoindentation on two planes of material symmetry were analyzed and the resulting longitudinal elastic moduli compared to the moduli measured with a macroscopic tensile test. A combined lathe and tensile system was designed that allows machining and testing of cylindrical microspecimens of approximately 4mm in length and 300 microm in diameter. Three bovine specimens were tested in tension and their outer geometry and porosity assessed by synchrotron radiation microtomography. Based on the results of the traction test and the precise outer geometry, an apparent longitudinal Young's modulus was calculated. Results between 20.3 and 27.6 GPa were found that match with previously reported values for bovine compact bone. The same specimens were then characterized by nanoindentation on a transverse and longitudinal plane. A longitudinal Young's modulus for the bone matrix was then derived using the numerical scheme proposed by Swadener and Pharr and the fabric-elasticity relationship by Zysset and Curnier. Based on the matrix modulus and a power law effective volume fraction, an apparent longitudinal Young's modulus was predicted for each microspecimen. This alternative approach provided values between 19.9 and 30.0 GPa, demonstrating differences between 2% and 13% to the values provided by the initial tensile test. This study therefore raises confidence in our nanoindentation protocol of the bone extracellular matrix and supports the underlying hypotheses used to extract the anisotropic elastic constants.     

Analytical methods to determine the effective mesoscopic and macroscopic elastic properties of cortical bone

We compare theoretical predictions of the effective elastic moduli of cortical bone at both the meso- and macroscales. We consider the efficacy of three alternative approaches: the method of asymptotic homogenization, the Mori-Tanaka scheme and the Hashin-Rosen bounds. The methods concur for specific engineering moduli such as the axial Young's modulus but can vary for others. In a past study, the effect of porosity alone on mesoscopic properties of cortical bone was considered, taking the matrix to be isotropic. Here, we consider the additional influence of the transverse isotropy of the matrix. We make the point that micromechanical approaches can be used in two alternative ways to predict either the macroscopic (size of cortical bone sample) or mesoscopic (in between micro- and macroscales) effective moduli, depending upon the choice of representative volume element size. It is widely accepted that the mesoscale behaviour is an important aspect of the mechanical behaviour of bone but models incorporating its effect have started to appear only relatively recently. Before this only macroscopic behaviour was addressed. Comparisons are drawn with experimental data and simulations from the literature for macroscale predictions with particularly good agreement in the case of dry bone. Finally, we show how predictions of the effective mesoscopic elastic moduli can be made which retain dependence on the well-known porosity gradient across the thickness of cortical bone.     

Apparent Young's modulus of human radius using inverse finite-element method

The ability to assess the elastic and failure properties of cortical bone at the radial diaphysis has a clinical importance. A new generation of quantitative ultrasound (QUS) devices and peripheral quantitative computed tomography (p-QCT) has been developed to assess non-invasively bone material and structural properties at the distal radius. This anatomical site is characterized by a thin cortical thickness that complicates traditional mechanical testing methods on specimens. Until now, mechanical properties of cortical bone at distal radius (e.g., elastic modulus, yield stress and strain) remain rarely studied probably due to experimental difficulties. The present study introduces an inverse finite-element method strategy to measure the elastic modulus and yield properties of human cortical specimens of the radial diaphysis. Twenty millimeter-thick portions of diaphysis were cut from 40 human radii (ages 45-90) for biomechanical test. Subsequently the same portion was modeled in order to obtain a specimen-specific three dimensional finite-element model (3D-FEM). Longitudinal elastic constants at the apparent level and stress characterizations were performed by coupling mechanical parameters with isotropic linear-elastic simulations. The results indicated that the mean apparent Young's modulus for radial cortical bone was 16 GPa (SD 1.8) and the yield stress was 153 MPa (SD 33). Breaking load was 12,946 N (SD 3644), cortical thickness 2.9 mm (SD 0.6), structural effective strain at the yield (epsilon(y)=0.0097) and failure (epsilon(u)=0.0154) load were also calculated. The 3D-FEM strategy described here may help to investigate bone mechanical properties when some difficulties arise from machining mechanical sample.     

Elastic moduli, yield stress, and ultimate stress of cancellous bone in the canine proximal femur

Elastic moduli, yield stress and ultimate compressive stress were determined for cancellous bone from the femoral head and neck regions of the canine femur. Unconfined compression tests were performed on 5 mm cubic samples which were cut from two femurs. Elastic moduli were measured in three orthogonal directions, and the yield stress and ultimate stress were measured along the proximal-distal axis. The results from this investigation support previous assumptions that the mechanical behavior of canine cancellous bone is qualitatively similar to human cancellous bone. The canine cancellous bone was observed to be anisotropic in elastic modulus. For two thirds of the cubic specimens tested, the elastic modulus was largest in the load-bearing, proximal-distal direction. A linear relationship between yield stress and elastic modulus was observed for canine bone, as is typical of human bone. A similar linear relationship between ultimate stress and elastic modulus was observed. Thus, for canine bone as well as for human bone, failure appears to be governed by a strain level which is position independent. The yield strain of 0.0259 and ultimate strain of 0.0288 for canine bone were both less than the yield strain of 0.0395 reported for human bone.     

The elastic moduli of human subchondral, trabecular, and cortical bone tissue and the size-dependency of cortical bone modulus

The elastic moduli of human subchondral, trabecular, and cortical bone tissue from a proximal tibia were experimentally determined using three-point bending tests on a microstructural level. The mean modulus of subchondral specimens was 1.15 GPa, and those of trabecular and cortical specimens was 4.59 GPa and 5.44 GPa respectively. Significant differences were found in the modulus values between bone tissues, which may have mainly resulted from the differences in the microstructures of each bone tissue rather than in the mineral density. Furthermore, the size-dependency of the modulus was examined using eight different sizes of cortical specimens (heights h = 100-1000 microns). While the modulus values for relatively large specimens (h greater than 500 microns) remained fairly constant (approximately 15 GPa), the values decreased as the specimens became smaller. A significant correlation was found between the modulus and specimen size. The surface area to volume ratio proved to be a key variable to explain the size-dependency.     

An alternative ultrasonic method for measuring the elastic properties of cortical bone

We studied the elastic properties of bone to analyze its mechanical behavior. The basic principles of ultrasonic methods are now well established for varying isotropic media, particularly in the field of biomedical engineering. However, little progress has been made in its application to anisotropic materials. This is largely due to the complex nature of wave propagation in these media. In the present study, the theory of elastic waves is essential because it relates the elastic moduli of a material to the velocity of propagation of these waves along arbitrary directions in a solid. Transducers are generally placed in contact with the samples which are often cubes with parallel faces that are difficult to prepare. The ultrasonic method used here is original, a rough preparation of the bone is sufficient and the sample is rotated. Moreover, to analyze heterogeneity of the structure we measure velocities in different points on the sample. The aim of the present study was to determine in vitro the anisotropic elastic properties of cortical bones. For this purpose, our method allowed measurement of longitudinal and transverse velocities (C(L) and C(T)) in longitudinal (fiber direction) and the radial directions (orthogonal to the fiber direction) of compact bones. Young's modulus E and Poisson's ratio nu, were then deduced from the velocities measured considering the compact bone as transversely isotropic or orthotropic. The results are in line with those of other methods.     

A cellular model of lung elasticity

The mechanics of the lung parenchyma is studied using models comprised of line members interconnected to form 3-D cellular structures. The mechanical properties are represented as elastic constants of a continuum. These are determined by perturbing each individual cell from a reference state by an increment in stress which is superimposed upon the uniform stretching forces initially present in the members due to the transpulmonary pressure. A force balance on the distorted structure, together with a force-deformation law for the members, leads to a calculation of the strain increments of the members. Predictions based on the analysis of the 3-D isotropic dodecahedron are in good agreement with experimental values for the Young's, shear, and bulk moduli reported in the literature. The model provides an explanation for the dependence of the elastic moduli on transpulmonary pressure, the geometrical details of the structure, and the stress-strain law of the tissue.     

Biomechanical effect of mineral heterogeneity in trabecular bone

Due to daily loading, trabecular bone is subjected to deformations (i.e., strain), which lead to stress in the bone tissue. When stress and/or strain deviate from the normal range, the remodeling process leads to adaptation of the bone architecture and its degree of mineralization to effectively withstand the sustained altered loading. As the apparent mechanical properties of bone are assumed to depend on the degree and distribution of mineralization, the goal of the present study was examine the influences of mineral heterogeneity on the biomechanical properties of trabecular bone in the human mandibular condyle. For this purpose nine right condyles from human dentate mandibles were scanned and evaluated with a microCT system. Cubic regional volumes of interest were defined, and each was transformed into two different types of finite element (FE) models, one homogeneous and one heterogeneous. In the heterogeneous models the element tissue moduli were scaled to the local degree of mineralization, which was determined using microCT. Compression and shear tests were simulated to determine the apparent elastic moduli in both model types. The incorporation of mineralization variation decreased the apparent Young's and shear moduli by maximally 21% in comparison to the homogeneous models. The heterogeneous model apparent moduli correlated significantly with bone volume fraction and degree of mineralization. It was concluded that disregarding mineral heterogeneity may lead to considerable overestimation of apparent elastic moduli in FE models.     

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