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.     

Static and dynamic finite element analyses of an idealized structural model of vertebral trabecular bone.

An idealized three-dimensional finite element model of a rodlike trabecular bone structure was developed to study its static and dynamic responses under compressive loading, considering the effects of bone marrow and apparent density. Static analysis of the model predicted hydraulic stiffening of trabecular bone due to the presence of bone marrow. The predicted power equation relating trabecular bone apparent elastic modulus to its apparent density was in good agreement with those of the reported experimental investigations. The ratio of the maximum stress in the trabecular bone tissue to its apparent stress had a high value, decreasing with increasing bone apparent density. Frequency analyses of the model predicted higher natural frequencies for the bone without marrow than those for the bone with marrow. Adding a mass relatively large compared to that of bone rendered a single-degree-of-freedom response. In this case, the resonant frequency was higher for the bone with marrow than that for the bone without marrow. The predicted vibrational measurement of apparent modulus was in good agreement with that of the static measurement, suggesting vibrational testing as a method for nondestructive measurement of trabecular bone elastic moduli.     

Biomechanical properties across trabeculae from the proximal femur of normal and ovariectomised sheep

The elastic behaviour of trabecular bone is a function not only of bone volume and architecture, but also of tissue material properties. Variation in tissue modulus can have a substantial effect on the biomechanical properties of trabecular bone. However, the nature of tissue property variation within a single trabecula is poorly understood. This study uses nanoindentation to determine the mechanical properties of bone tissue in individual trabeculae. Using an ovariectomised ovine model, the modulus and hardness distribution across trabeculae were measured. In both normal and ovariectomised bone, the modulus and hardness were found to increase towards the core of the trabeculae. Across the width of the trabeculae, the modulus was significantly less in the ovariectomised bone than in the control bone. However, in contrast to this hardness was found not to differ significantly between the two groups. This study provides valuable information on the variation of mechanical material properties in healthy and diseased trabecular bone tissue. The results of the current study will be useful in finite element modelling where more accurate values of trabecular bone modulus will enable the prediction of the macroscale behaviour of trabecular 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.     

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.     

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.     

Incorporating tissue anisotropy and heterogeneity in finite element models of trabecular bone altered predicted local stress distributions

Trabecular bone is composed of organized mineralized collagen fibrils, which results in heterogeneous and anisotropic mechanical properties at the tissue level. Recently, biomechanical models computing stresses and strains in trabecular bone have indicated a significant effect of tissue heterogeneity on predicted stresses and strains. However, the effect of the tissue-level mechanical anisotropy on the trabecular bone biomechanical response is unknown. Here, a computational method was established to automatically impose physiologically relevant orientation inherent in trabecular bone tissue on a trabecular bone microscale finite element model. Spatially varying tissue-level anisotropic elastic properties were then applied according to the bone mineral density and the local tissue orientation. The model was used to test the hypothesis that anisotropy in both homogeneous and heterogeneous models alters the predicted distribution of stress invariants. Linear elastic finite element computations were performed on a 3 mm cube model isolated from a microcomputed tomography scan of human trabecular bone from the distal femur. Hydrostatic stress and von Mises equivalent stress were recorded at every element, and the distributions of these values were analyzed. Anisotropy reduced the range of hydrostatic stress in both tension and compression more strongly than the associated increase in von Mises equivalent stress. The effect of anisotropy was independent of the spatial redistribution high compressive stresses due to tissue elastic heterogeneity. Tissue anisotropy and heterogeneity are likely important mechanisms to protect bone from failure and should be included for stress analyses in trabecular bone.     

Side-artifact errors in yield strength and elastic modulus for human trabecular bone and their dependence on bone volume fraction and anatomic site

In the context of reconciling the mechanical properties of trabecular bone measured from in vitro mechanical testing with the true in situ behavior, recent attention has focused on the "side-artifact" which results from interruption of the trabecular network along the sides of machined specimens. The objective of this study was to compare the magnitude of the side-artifact error for measurements of elastic modulus vs. yield stress and to determine the dependence of these errors on anatomic site and trabecular micro-architecture. Using a series of parametric variations on micro-CT-based finite element models of trabecular bone from the human vertebral body (n=24) and femoral neck (n=10), side-artifact correction factors were quantified as the ratio of the side-artifact-free apparent mechanical property to the corresponding property measured in a typical experiment. The mean (+/-SD) correction factors for yield stress were 1.32+/-0.17 vs. 1.20+/-0.11 for the vertebral body and femoral neck (p<0.05), respectively, and the corresponding factors for modulus were 1.24+/-0.09 vs. 1.10+/-0.04 (p<0.0001). Correction factors were greater for yield stress than modulus (p<0.003), but no anatomic site effect was detected (p>0.29) after accounting for variations in bone volume fraction (BV/TV). Approximately 30-55% of the variation in the correction factors for modulus and yield stress could be accounted for by BV/TV or micro-architecture, representing an appreciable systematic component of the error. Although some scatter in the correction factor-BV/TV relationships may confound accurate correction of modulus and yield stress for individual specimens, side-artifact correction is nonetheless essential for obtaining accurate mean estimates of modulus and yield stress for a cohort of specimens. We conclude that appreciation and correction for the differential effects of the side-artifact in modulus vs. yield stress and their dependence on BV/TV may improve the interpretation of measured elastic and failure properties for trabecular bone.     

Effects of freeze–thaw and micro-computed tomography irradiation on structure–property relations of porcine trabecular bone

We study the effects of freeze–thaw and irradiation on structure–property relations of trabecular bone. We measure the porosity, apparent density, mineral content, trabecular orientation, trabecular thickness, fractal dimension, surface area, and connectivity of trabecular bone using micro-computed tomography (micro-CT) and relate them to Young?s modulus and ultimate strength measured by uniaxial compression testing. The analysis is done on six-month porcine trabecular bone from femoral heads. The effects of freeze–thaw are studied by using bones from three different groups: fresh bone and bones frozen for one and five years. We find that the porosity and apparent density have most dominant influence on the elastic modulus and strength of fresh bone. Also, five years of freezing lowers both Young?s modulus and ultimate strength of trabecular bone. Additionally, the effects of radiation are investigated by comparing Young?s modulus before and after micro-CT exposure. We find that the micro-CT irradiation has a negligible effect on the Young?s modulus of trabecular bone. These findings provide insights on the effects of tissue preservation and imaging on properties of trabecular bone.     

Mechanical properties of femoral trabecular bone in dogs

Background  

Studying mechanical properties of canine trabecular bone is important for a better understanding of fracture mechanics or bone disorders and is also needed for numerical simulation of canine femora. No detailed data about elastic moduli and degrees of anisotropy of canine femoral trabecular bone has been published so far, hence the purpose of this study was to measure the elastic modulus of trabecular bone in canine femoral heads by ultrasound testing and to assess whether assuming isotropy of the cancellous bone in femoral heads in dogs is a valid simplification.     

Nonlinear behavior of trabecular bone at small strains

Study of the behavior of trabecular bone at strains below 0.40 percent is of clinical and biomechanical importance. The goal of this work was to characterize, with respect to anatomic site, loading mode, and apparent density, the subtle concave downward stress-strain nonlinearity, that has been observed recently for trabecular bone at these strains. Using protocols designed to minimize end-artifacts, 155 cylindrical cores from human vertebrae, proximal tibiae, proximal femora, and bovine proximal tibiae were mechanically tested to yield at 0.50 percent strain per second in tension or compression. The nonlinearity was quantified by the reduction in tangent modulus at 0.20 percent and 0.40 percent strain as compared to the initial modulus. For the pooled data, the mean +/- SD percentage reduction in tangent modulus at 0.20 percent strain was 9.07+/- 3.24 percent in compression and 13.8 +/- 4.79 percent in tension. At 0.40 percent strain, these values were 23.5 +/- 5.71 and 35.7+/- 7.10 percent, respectively. The magnitude of the nonlineari't depended on both anatomic site (p < 0.001) and loading mode (p < 0.001), and in tension was positively correlated with density. Calculated values of elastic modulus and yield properties depended on the strain range chosen to define modulus via a linear curve fit (p < 0.005). Mean percent differences in 0.20 percent offset yield strains were as large as 10.65 percent for some human sites. These results establish that trabecular bone exhibits nonlinearity at low strains, and that this behavior can confound intersite comparisons of mechanical properties. A nonlinear characterization of the small strain behavior of trabecular bone was introduced to characterize the initial stress-strain behavior more thoroughly.     

Fast and accurate specimen-specific simulation of trabecular bone elastic modulus using novel beam-shell finite element models

Elastic modulus and strength of trabecular bone are negatively affected by osteoporosis and other metabolic bone diseases. Micro-computed tomography-based beam models have been presented as a fast and accurate way to determine bone competence. However, these models are not accurate for trabecular bone specimens with a high number of plate-like trabeculae. Therefore, the aim of this study was to improve this promising methodology by representing plate-like trabeculae in a way that better reflects their mechanical behavior. Using an optimized skeletonization and meshing algorithm, voxel-based models of trabecular bone samples were simplified into a complex structure of rods and plates. Rod-like and plate-like trabeculae were modeled as beam and shell elements, respectively, using local histomorphometric characteristics. To validate our model, apparent elastic modulus was determined from simulated uniaxial elastic compression of 257 cubic samples of trabecular bone (4mm×4mm×4mm; 30μm voxel size; BIOMED I project) in three orthogonal directions using the beam-shell models and using large-scale voxel models that served as the gold standard. Excellent agreement (R(2)=0.97) was found between the two, with an average CPU-time reduction factor of 49 for the beam-shell models. In contrast to earlier skeleton-based beam models, the novel beam-shell models predicted elastic modulus values equally well for structures from different skeletal sites. It allows performing detailed parametric analyses that cover the entire spectrum of trabecular bone microstructures.     

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.     

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.     

Effects of thresholding techniques on microCT-based finite element models of trabecular bone

Microimaging based finite element analysis is widely used to predict the mechanical properties of trabecular bone. The choice of thresholding technique, a necessary step in converting grayscale images to finite element models, can significantly influence the predicted bone volume fraction and mechanical properties. Therefore, we investigated the effects of thresholding techniques on microcomputed tomography (micro-CT) based finite element models of trabecular bone. Three types of thresholding techniques were applied to 16-bit micro-CT images of trabecular bone to create three different models per specimen. Bone volume fractions and apparent moduli were predicted and compared to experimental results. In addition, trabecular tissue mechanical parameters and morphological parameters were compared among different models. Our findings suggest that predictions of apparent mechanical properties and structural properties agree well with experimental measurements regardless of the choice of thresholding methods or the format of micro-CT images.     

A new multiscale micromechanical model of vertebral trabecular bones

A new three-dimensional (3D) multiscale micromechanical model has been suggested as adept at predicting the overall linear anisotropic mechanical properties of a vertebral trabecular bone (VTB) highly porous microstructure. A nested 3D modeling analysis framework spanning the multiscale nature of the VTB is presented herein. This hierarchical analysis framework employs the following micromechanical methods: the 3D parametric high-fidelity generalized method of cells (HFGMC) as well as the 3D sublaminate model. At the nanoscale level, the 3D HFGMC method is applied to obtain the effective elastic properties of a representative unit cell (RUC) representing the mineral collagen fibrils composite. Next at the submicron scale level, the 3D sublaminate model is used to generate the effective elastic properties of a repeated stack of multilayered lamellae demonstrating the nature of the trabeculae (bone-wall). Thirdly, at the micron scale level, the 3D HFGMC method is used again on a RUC of the highly porous VTB microstructure. The VTB-RUC geometries are taken from microcomputed tomography scans of VTB samples harvested from different vertebrae of human cadavers \((n=10)\). The predicted anisotropic overall elastic properties for native VTBs are, then, examined as a function of age and sex. The predicted results of the VTBs longitudinal Young’s modulus are compared to reported values found in the literature. The proposed 3D nested modeling analysis framework provides a good agreement with reported values of Young’s modulus of single trabeculae as well as for VTB-RUC in the literature.     

A morphological model of vertebral trabecular bone

In their micro-structures, typical natural cellular materials such as vertebral trabecular bone have a network of doubly tapered struts, thickening near the strut joints. However, past analytical models for vertebral trabecular bone do not take account of the effect of strut taper on the mechanical properties.This paper presents an analytical cell model comprised of doubly tapered struts to predict the global mechanical properties of vertebral trabecular bone. The predicted results for male, female, and both sexes fit the experimental data well. By considering several strut taper geometries, it is shown that the horizontal Young's modulus and the horizontal uniaxial collapse stress are, in some cases, approximately 1.8- and 2.2-fold higher, respectively, than those of the uniform strut model. This finding illustrates the importance of increased trabecular thickening near the strut joints (i) for improving the accuracy of calculating the mechanical properties and (ii) for the effective treatment of aged bone using drug therapy. It also highlights the need to combine trabecular architecture measurements with information about the morphology near the strut joints.     

Modeling the onset and propagation of trabecular bone microdamage during low-cycle fatigue

Relatively small amounts of microdamage have been suggested to have a major effect on the mechanical properties of bone. A significant reduction in mechanical properties (e.g. modulus) can occur even before the appearance of microcracks. This study uses a novel non-linear microdamaging finite-element (FE) algorithm to simulate the low-cycle fatigue behavior of high-density trabecular bone. We aimed to investigate if diffuse microdamage accumulation and concomitant modulus reduction, without the need for complete trabecular strut fracture, may be an underlining mechanism for low-cycle fatigue failure (defined as a 30% reduction in apparent modulus). A microCT constructed FE model was subjected to a single cycle monotonic compression test, and constant and variable amplitude loading scenarios to study the initiation and accumulation of low-cycle fatigue microdamage. Microcrack initiation was simulated using four damage criteria: 30%, 40%, 50% and 60% reduction in bone element modulus (el-MR). Evaluation of structural (apparent) damage using the four different tissue level damage criteria resulted in specimen fatigue failure at 72, 316, 969 and 1518 cycles for the 30%, 40%, 50% and 60% el-MR models, respectively. Simulations based on the 50% el-MR model were consistent with previously published experimental findings. A strong, significant non-linear, power law relationship was found between cycles to failure (N) and effective strain (Deltasigma/E(0)): N=1.394x10(-25)(Deltasigma/E(0))(-12.17), r(2)=0.97, p<0.0001. The results suggest that microdamage and microcrack propagation, without the need for complete trabecular strut fracture, are mechanisms for high-density trabecular bone failure. Furthermore, the model is consistent with previous numerical fatigue simulations indicating that microdamage to a small number of trabeculae results in relatively large specimen modulus reductions and rapid failure.     

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.     

How morphology predicts mechanical properties of trabecular structures depends on intra-specimen trabecular thickness variations

Two observations underlie this work. First, that the architecture of trabecular bone can accurately predict the mechanical stiffness characteristics of bone specimens when considering the combination of volume fraction and fabric, which is a measure of architectural anisotropy. Second, that the same morphological measures could not accurately predict the mechanical properties of porous structures in general. We hypothesize that this discrepancy can be explained by the special nature of trabecular bone as a structure in remodeling equilibrium relative to the external loads. We tested this hypothesis using a generic model of trabecular bone. Five series of 153 different architectures were created with this model. Each architecture was subjected to morphological analysis, and four different fabric measures were calculated to evaluate their effectiveness in characterizing the architecture. Relationships were determined relating morphology to the elastic constants. The quality of these relationships was tested by correlating the predicted elastic constants with those determined from finite element analysis. We found that the four fabric measures used could estimate the mechanical properties almost equally well. So the suggestion that fabric measures based on trabecular bone volume better represent the architecture than mean intercept length could not be affirmed. We conclude that for structures with equally sized elliptical voids the mechanical properties can be predicted well only if trabecular thickness variations within each structure are limited. These structures closely resemble previously developed models of trabecular bone. Furthermore, they are stiff in the principal fabric direction, hence, according to Cowin (J. Biomech. Eng. (108) (1986) 83), they are in remodeling equilibrium. These structures are also stiff over a large range of loading orientations, hence, are relatively insensitive to deviations in direction of loading.