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.     

Validation of a voxel-based FE method for prediction of the uniaxial apparent modulus of human trabecular bone using macroscopic mechanical tests and nanoindentation

Assessment of the mechanical properties of trabecular bone is of major biological and clinical importance for the investigation of bone diseases, fractures and their treatments. Finite element (FE) methods are getting increasingly popular for quantifying the elastic and failure properties of trabecular bone. In particular, voxel-based FE methods have been previously used to calculate the effective elastic properties of trabecular microstructures. However, in most studies, bone tissue moduli were assumed or back-calculated to match the apparent elastic moduli from experiments, which often lead to surprisingly low values when compared to nanoindentation results. In this study, voxel-based FE analysis of trabecular bone is combined with physical measures of volume fraction, micro-CT (microCT) reconstructions, uniaxial mechanical tests and specimen-specific nanoindentation tests for proper validation of the method. Cylindrical specimens of cancellous bone were extracted from human femurs and their volume fraction determined with Archimede's method. Uniaxial apparent modulus of the specimens was measured with an improved tension-compression testing protocol that minimizes boundary artefacts. Their microCT reconstructions were segmented to match the measured bone volume fraction and used to create full-size voxel models with 30-45 microm element size. For each specimen, linear isotropic elastic material properties were defined based on specific nanoindentation measurements of its embedded bone tissue. Linear FE analyses were finally performed to simulate the uniaxial mechanical tests. Additional parametric analyses were performed to evaluate the potential errors on the predicted apparent modulus arising from variations in segmentation threshold, tissue modulus, and the use of 125-mm(3) cubic sub-regions. The results demonstrate an excellent correspondence between experimental measures and FE predictions of uniaxial apparent modulus. In conclusion, the adopted voxel-based FE approach is found to be a robust method to predict the linear elastic properties of human cancellous bone, provided segmentation of the microCT reconstructions is carefully calibrated, tissue modulus is known a priori and the entire region of interest is included in the analysis.     

High-resolution finite element models with tissue strength asymmetry accurately predict failure of trabecular bone

The ability to predict trabecular failure using microstructure-based computational models would greatly facilitate study of trabecular structure–function relations, multiaxial strength, and tissue remodeling. We hypothesized that high-resolution finite element models of trabecular bone that include cortical-like strength asymmetry at the tissue level, could predict apparent level failure of trabecular bone for multiple loading modes. A bilinear constitutive model with asymmetric tissue yield strains in tension and compression was applied to simulate failure in high-resolution finite element models of seven bovine tibial specimens. Tissue modulus was reduced by 95% when tissue principal strains exceeded the tissue yield strains. Linear models were first calibrated for effective tissue modulus against specimen-specific experimental measures of apparent modulus, producing effective tissue moduli of (mean±S.D.) 18.7±3.4 GPa. Next, a parameter study was performed on a single specimen to estimate the tissue level tensile and compressive yield strains. These values, 0.60% strain in tension and 1.01% strain in compression, were then used in non-linear analyses of all seven specimens to predict failure for apparent tensile, compressive, and shear loading. When compared to apparent yield properties previously measured for the same type of bone, the model predictions of both the stresses and strains at failure were not statistically different for any loading case (p>0.15). Use of symmetric tissue strengths could not match the experimental data. These findings establish that, once effective tissue modulus is calibrated and uniform but asymmetric tissue failure strains are used, the resulting models can capture the apparent strength behavior to an outstanding level of accuracy. As such, these computational models have reached a level of fidelity that qualifies them as surrogates for destructive mechanical testing of real specimens.     

Postfailure modulus strongly affects microcracking and mechanical property change in human iliac cancellous bone: a study using a 2D nonlinear finite element method

A two-dimensional (2D) finite element (FE) method was used to estimate the ability of bone tissue to sustain damage as a function of postfailure modulus. Briefly, 2D nonlinear compact-tension FE models were created from quantitative back-scattered electron images taken of human iliac crest bone specimens. The effects of different postfailure moduli on predicted microcrack propagation were examined. The 2D FE models were used as surrogates for real bone tissues. The crack number was larger in models with higher postfailure modulus, while mean crack length and area were smaller in these models. The rate of stiffness reduction was greater in the models with lower postfailure modulus. Hence, the current results supported the hypothesis that hard tissue postfailure properties have strong effects on bone microdamage morphology and the rate of change in apparent mechanical properties.     

Anatomical variation of orthotropic elastic moduli of the proximal human tibia

The anatomical variation of orthotropic elastic moduli of the cancellous bone from three human proximal tibiae was investigated using an ultrasonic technique. With this technique, it was possible to measure three orthogonal elastic moduli and three shear moduli from cubic specimens of cancellous bone as small as 8 mm per side. Correlation with mechanical tensile testing has shown this technique to offer a precise measure of cancellous modulus (Eten = 0.94Eult + 144.6 MPa, r2 = 0.96, n = 34). The cancellous bone of the proximal tibia was found to be very inhomogeneous, with the axial modulus ranging between 340 and 3350 MPa. A course map is presented, showing measured Young's moduli as a function of anatomical position. The anisotropy of the cancellous bone, determined by the relative differences between the three orthogonal moduli, was shown to be relatively constant over the entire range of cancellous densities tested. The relationship between the axial elastic modulus and the apparent density was found to be approximately linear, as reported by others for proximal tibial cancellous bone.     

The underestimation of Young's modulus in compressive testing of cancellous bone specimens

In order to determine the accuracy of measurements of Young's modulus of cancellous bone by conventional compression testing, two independent strain measurements were made simultaneously during non-destructive uniaxial compression to 0.8% strain of rectangular specimens (n = 18). Strain was measured by an extensometer attached to the compression anvils close to the specimen and by an optical system covering the central half of the specimens. Mean Young's modulus determined by the extensometer technique was 689 MPa, but was 871 MPa when determined by the optical technique (mean difference = 182 MPa, SED = 50 MPa, p less than 0.002). Uneven strain distribution due to lack of support of cut vertical trabeculae at the anvil-specimen interface is believed to be causing the underestimation of Young's modulus measured by the extensometer technique. The influence of friction at the specimen-anvil interface was studied by performing a finite element analysis. It is concluded that Young's modulus of specimens of the chosen geometry on average is underestimated by about 20% by conventional compressing testing. The underestimation seems not to be dependent upon specimen density.     

Effects of damage on the orthotropic material symmetry of bovine tibial trabecular bone

The macroscopic mechanical properties of trabecular bone can be predicted by its architecture using theoretical relationships between the elastic and architectural properties. Microdamage caused by overloading or fatigue decreases the apparent elastic moduli of trabecular bone requiring these relationships to be modified to predict the damaged elastic properties. In the case of isotropic damage, the apparent level elastic properties could be determined by multiplying all of the elastic constants by a single scalar factor. If the damage is anisotropic, the elastic constants may change by differing factors and the material coordinate system could become misaligned with the fabric coordinate system. High-resolution finite element models were used to simulate damage overloading on seven trabecular bone specimens subjected to pure shear strain in two planes. Comparison of the apparent elastic moduli of the specimens before and after damage showed that the reduction of the elastic moduli was anisotropic. This suggests that the microdamage within the specimens was inhomogeneous. However, after damage the specimens exhibited nearly orthotropic material symmetry as they did before damage. Changes in the orientation of the orthotropic material coordinate system were also small and occurred primarily in the transverse plane. Thus, while damage in trabecular bone is anisotropic, the material coordinate system remains aligned with the fabric tensor.     

Heterogeneous linear elastic trabecular bone modelling using micro-CT attenuation data and experimentally measured heterogeneous tissue properties

High-resolution voxel-based finite element software, such as FEEBE developed at the NCBES, is widely used for studying trabecular bone at the micro-scale. A new approach to determine heterogeneous bone tissue material properties for computational models was proposed in this study. The specimen-specific range of tissue moduli across strut width was determined from nanoindentation testing. This range was mapped directly using linear interpolation to that specimen's micro-computed tomography (microCT) grey value range as input material properties for finite element analysis. The method was applied to cuboid trabecular bone samples taken from eight, 4-year-old (skeletally mature) ovine L5 vertebrae. Before undergoing experimental uniaxial compression tests, the samples were microCT scanned and 30 microm resolution finite element models were generated. The linear elastic finite element models were compressed to 1% strain. This material property assignment method for computational models accurately reproduced the experimentally determined apparent modulus and concentrations of stress at locations of failure.     

Apparent Young's modulus of vertebral cortico-cancellous bone specimens

Up to now, due to cortical thickness and imaging resolution, it is not possible to derive subject-specific mechanical properties on the 'vertebral shell' from imaging modalities applicable in vivo. As a first step, the goal of this study was to assess the apparent Young's modulus of vertebral cortico-cancellous bone specimens using an inverse method. A total of 22 cortico-cancellous specimens were harvested from 22 vertebral bodies. All specimens were tested in compression until failure. To compute the apparent Young's modulus of the specimen from the inverse method, the boundary conditions of the biomechanical experiments were faithfully reproduced in a finite element model (FEM), and an optimisation routine was used. The results showed a mean of the apparent Young's modulus of 374?±?208?MPa, ranging from 87 to 791?MPa. By computing an apparent Young's modulus of a cortico-cancellous medium, this study gives mechanical data for an FEM of an entire vertebra including an external shell combining both bone tissues.     

Elastic modulus of trabecular bone material

An ultrasonic technique was used to measure both the elastic modulus (Young's modulus) of trabecular bone material and the elastic modulus of the cancellous structure. The average trabecular modulus, measured on specimens obtained from three human and one bovine distal femora, was 13.0 GPa (S.D. 1.47) and 10.9 GPa (S.D. 1.57), respectively. On human specimens the structural elastic modulus was found to be related to the structural (apparent) density raised to the 1.88 power. The elastic modulus from the bovine specimens showed a more linear relationship with the density of the cancellous structure (density raised to the 1.57 power).     

Glenoid cancellous bone strength and modulus.

The objectives of this study were to determine the strength and modulus of glenoid cancellous bone, including regional variations. The motivations were: to select a suitable bone substitute for standardized testing of glenoid prosthesis loosening, to assist in shoulder prosthesis design and to provide input data for finite element analyses. Ten glenoids from eight cadavers (mean age, 81) were tested by in situ indentation. Mean strength ranged from 6.7 to 17 MPa for the ten glenoids, the overall mean being 10.3 MPa. Mean E moduli ranged from 67 to 171 MPa for the individual glenoids, the overall mean being 99 MPa. These values are likely at the lower end of what would be expected for normal bone since strength and modulus decrease with age and the available specimens were older. These values may be appropriate for prosthesis design, however, since mechanical properties are reduced in rheumatoid arthritic bone. Regional trends were very similar for modulus and strength. The strongest region was postero-superior. The central column, correlating with the keel position in many glenoid components, was weaker than both the anterior and posterior regions but deeper. A large drop in strength and modulus below the subchondral layer emphasizes the importance of maintaining this layer during prosthetic replacement.     

NACOB presentation to ASB Young Scientist Award: Postdoctoral. The impact of boundary conditions and mesh size on the accuracy of cancellous bone tissue modulus determination using large-scale finite-element modeling. North American Congress on Biomechanics.

The apparent properties of cancellous bone are determined by a combination of both hard tissue properties and microstructural organization. A method is desired to extract the underlying hard tissue properties from simple mechanical tests, free from the complications of microstructure. It has been suggested that microCT voxel-based large-scale finite element models could be employed to accomplish this goal (van Rietbergen et al., 1995, Journal of Biomechanics, 28, 69-81). This approach has recently been implemented and it is becoming increasingly popular as finite element models increase in size and sophistication (Fyhrie et al., 1997, Proceedings of the 43rd Annual Meeting of the Orthopaedic Research Society, San Francisco, CA, p. 815; van Rietbergen et al., 1997, Proceedings of the 43rd Annual Meeting of the Orthopaedic Research Society, San Francisco, CA, p. 62). However, no direct quantitative measurements of the accuracy of this method applied to porous structures such as cancellous bone have been made. This project demonstrates the feasibility of this approach by quantifying its best-case accuracy in determining the trabecular hard tissue modulus of analogues fabricated of a material with known material properties determined independently by direct testing. In addition we were able to assess the impact of mesh size and boundary conditions on accuracy. We found that the assumption of a frictionless boundary condition in the parallel plate compression loading configuration was a significant source of error that could be overcome with the use of rigid end-caps similar to those used by Keaveny et al. (1997 Journal of Orthopaedic Research, 15(1), 101-110). In conclusion, we found that this approach is an effective method for determining the average trabecular hard tissue properties of human cancellous bone with an expected practical accuracy level better than 5%.     

Mineral heterogeneity has a minor influence on the apparent elastic properties of human cancellous bone: a SRμCT-based finite element study

At the tissue level, the local material properties of human cancellous bone are heterogeneous due to constant remodelling. Since standard high-resolution computed tomography scanning methods are unable to capture this heterogeneity in detail, local differences in mineralisation are normally not incorporated in computational models. To investigate the effects of heterogeneous mineral distribution on the apparent elastic properties, 40 cancellous bone samples from the human femoral neck were scanned by means of synchrotron radiation microcomputed tomography (SRμCT). SRμCT-based micromechanical finite element models that accounted for mineral heterogeneity were compared with homogeneous models. Evaluation of the apparent stiffness tensor of both model types revealed that homogeneous models led to a minor but significant (p < 0.05) overestimation of the elastic properties of heterogeneous models by 2.18 ± 1.89%. Variation of modelling parameters did not affect the overestimation to a great extent. It was concluded that the heterogeneous mineralisation has only a minor influence on the apparent elastic properties of human cancellous bone.     

Convergence behavior of high-resolution finite element models of trabecular bone

The convergence behavior of finite element models depends on the size of elements used, the element polynomial order, and on the complexity of the applied loads. For high-resolution models of trabecular bone, changes in architecture and density may also be important. The goal of this study was to investigate the influence of these factors on the convergence behavior of high-resolution models of trabecular bone. Two human vertebral and two bovine tibial trabecular bone specimens were modeled at four resolutions ranging from 20 to 80 microns and subjected to both compressive and shear loading. Results indicated that convergence behavior depended on both loading mode (axial versus shear) and volume fraction of the specimen. Compared to the 20 microns resolution, the differences in apparent Young's modulus at 40 microns resolution were less than 5 percent for all specimens, and for apparent shear modulus were less than 7 percent. By contrast, differences at 80 microns resolution in apparent modulus were up to 41 percent, depending on the specimen tested and loading mode. Overall, differences in apparent properties were always less than 10 percent when the ratio of mean trabecular thickness to element size was greater than four. Use of higher order elements did not improve the results. Tissue level parameters such as maximum principal strain did not converge. Tissue level strains converged when considered relative to a threshold value, but only if the strains were evaluated at Gauss points rather than element centroids. These findings indicate that good convergence can be obtained with this modeling technique, although element size should be chosen based on factors such as loading mode, mean trabecular thickness, and the particular output parameter of interest.     

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.     

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.     

Mechanisms of uniformity of yield strains for trabecular bone

Variations in yield strains for trabecular bone within a specific anatomic site are only a small fraction of the substantial variations that exist for elastic modulus and strength, and yet the source of this uniformity is not known. Our goal was to investigate the underlying mechanisms by using high-resolution, materially nonlinear finite element models of 12 human femoral neck trabecular bone specimens. The finite element models, used to obtain apparent yield strains in both tension and compression, assumed that the tissue-level yield strains were the same across all specimens. Comparison of the model predictions with the experimental data therefore enabled us to isolate the combined roles of volume fraction and architecture from the role of tissue material properties. Results indicated that, for both tensile and compressive loading, natural variations in volume fraction and architecture produced a negligible coefficient of variation (less than 3%) in apparent yield strains. Analysis of tissue-level strains showed that while bending of individual trabeculae played only a minor role in the apparent elastic behavior, the combined effects of this bending and tissue-level strength asymmetry produced apparent-level failure strains in compression that were 14% lower than those at the tissue level. By contrast, tissue and apparent-level yield strains were equivalent for tensile loading. We conclude that the uniformity of apparent yield strains is primarily the result of the highly oriented architecture that minimizes bending. Most of the variation that does occur is the result of the non-uniformity of the tissue-level yield strains.     

Trabecular bone modulus-density relationships depend on anatomic site

One outstanding issue regarding the relationship between elastic modulus and density for trabecular bone is whether the relationship depends on anatomic site. To address this, on-axis elastic moduli and apparent densities were measured for 142 specimens of human trabecular bone from the vertebra (n=61), proximal tibia (n=31), femoral greater trochanter (n=23), and femoral neck (n=27). Specimens were obtained from 61 cadavers (mean+/-SD age=67+/-15 years). Experimental protocols were used that minimized end-artifact errors and controlled for specimen orientation. Tissue moduli were computed for a subset of 18 specimens using high-resolution linear finite element analyses and also using two previously developed theoretical relationships (Bone 25 (1999) 481; J. Elasticity 53 (1999) 125). Resultant power law regressions between modulus and density did depend on anatomic site, as determined via an analysis of covariance. The inter-site differences were among the leading coefficients (p<0.02), but not the exponents (p>0.08), which ranged 1.49-2.18. At a given density, specimens from the tibia had higher moduli than those from the vertebra (p=0.01) and femoral neck (p=0.002); those from the trochanter had higher moduli than the vertebra (p=0.02). These differences could be as large as almost 50%, and errors in predicted values of modulus increased by up to 65% when site-dependence was ignored. These results indicate that there is no universal modulus-density relationship for on-axis loading. Tissue moduli computed using methods that account for inter-site architectural variations did not differ across site (p>0.15), suggesting that the site-specificity in apparent modulus-density relationships may be attributed to differences in architecture.     

Long bone torsion: II. A combined experimental and computational method for determining an effective shear modulus

A technique is established which allows an effective torsional shear modulus to be determined for long bones, while remaining nondestructive to whole bone specimens. Strain gages are bonded to the diaphysis of the bone. Strains are then recorded under pure torsional loads. Theoretical stress predictions are combined with experimental strain recordings to arrive at a modulus value. Shear modulus calculations for four canine radii are reported using theoretical stress predictions from circular, elliptical and finite element models of the transverse bone geometry. The effective shear modulus, obtained from an average of the shear moduli determined at strain gage locations, serves to average the heterogeneous shear modulus distribution over the cross section. The shear modulus obtained is that associated with the "circumferential" direction in transverse planes.