The effects of side-artifacts on the elastic modulus of trabecular bone

Determining accurate density-mechanical property relationships for trabecular bone is critical for correct characterization of this important structure-function relation. When testing any excised specimen of trabecular bone, an unavoidable experimental artifact originates from the sides of the specimen where peripheral trabeculae lose their vertical load-bearing capacity due to interruption of connectivity, a phenomenon denoted here as the 'side-artifact'. We sought in this study to quantify the magnitude of such side-artifact errors in modulus measurement and to do so as a function of the trabecular architecture and specimen size. Using parametric computational analysis of high-resolution micro-CT-based finite-element models of cores of elderly human vertebral trabecular bone, a specimen-specific correction factor for the side-artifact was quantified as the ratio of the side-artifact-free apparent modulus (Etrue) to the apparent modulus that would be measured in a typical experiment (Emeasured). We found that the width over which the peripheral trabeculae were mostly unloaded was between 0.19 and 0.58 mm. The side-artifact led to an underestimation error in Etrue of over 50% in some specimens, having a mean (+/-SD) of 27+/-11%. There was a trend for the correction factor to linearly increase as volume fraction decreased (p=0.001) and as mean trabecular separation increased (p<0.001). Further analysis indicated that the error increased substantially as specimen size decreased. Two methods used for correcting for the side-artifact were both successful in bringing Emeasured into statistical agreement with Etrue. These findings have important implications for the interpretation of almost all literature data on trabecular bone mechanical properties since they indicate that such properties need to be adjusted to eliminate the substantial effects of side-artifacts in order to provide more accurate estimates of in situ behavior.     

Dependence of yield strain of human trabecular bone on anatomic site

Understanding the dependence of human trabecular bone strength behavior on anatomic site provides insight into structure-function relationships and is essential to the increased success of site-specific finite element models of whole bones. To investigate the hypothesis that the yield strains of human trabecular bone depend on anatomic site, the uniaxial tensile and compressive yield properties were compared for cylindrical specimens from the vertebra (n=61), proximal tibia (n=31), femoral greater trochanter (n=23), and femoral neck (n=27) taken from 61 donors (67+/-15years). Test protocols were used that minimized end artifacts and loaded specimens along the main trabecular orientation. Yield strains by site (mean+/-S.D.) ranged from 0.70+/-0.05% for the trochanter to 0.85+/-0.10% for the femoral neck in compression, from 0.61+/-0.05% for the trochanter to 0.70+/-0.05% for the vertebra in tension, and were always higher in compression than tension (p<0.001). The compressive yield strain was higher for the femoral neck than for all other sites (p<0.001), as was the tensile yield strain for the vertebra (p<0.007). Analysis of covariance, with apparent density as the covariate, indicated that inter-site differences existed in yield stress even after adjusting statistically for density (p<0.035). Coefficients of variation in yield strain within each site ranged from only 5-12%, consistent with the strong linear correlations (r(2)=0.94-0.98) found between yield stress and modulus. These results establish that the yield strains of human trabecular bone can differ across sites, but that yield strain may be considered uniform within a given site despite substantial variation in elastic modulus and yield stress.     

Heterogeneity of yield strain in low-density versus high-density human trabecular bone

Understanding the off-axis behavior of trabecular yield strains may lend unique insight into the etiology of fractures since yield strains provide measures of failure independent of elastic behavior. We sought to address anisotropy of trabecular yield strains while accounting for variations in both density and anatomic site and to determine the mechanisms governing this behavior. Cylindrical specimens were cored from vertebral bodies (n=22, BV/TV=0.11±0.02) and femoral necks (n=28, BV/TV=0.22±0.06) with the principal trabecular orientation either aligned along the cylinder axis (on-axis, n=22) or at an oblique angle of 15° or 45° (off-axis, n=28). Each specimen was scanned with micro-CT, mechanically compressed to failure, and analysed with nonlinear micro-CT-based finite element analysis. Yield strains depended on anatomic site (p=0.03, ANOVA), and the effect of off-axis loading was different for the two sites (p=0.04)—yield strains increased for off-axis loading of the vertebral bone (p=0.04), but were isotropic for the femoral bone (p=0.66). With sites pooled together, yield strains were positively correlated with BV/TV for on-axis loading (R²=58%, p<0.0001), but no such correlation existed for off-axis loading (p=0.79). Analysis of the modulus-BV/TV and strength-BV/TV relationships indicated that, for the femoral bone, the reduction in strength associated with off-axis loading was greater than that for modulus, while the opposite trend occurred for the vertebral bone. The micro-FE analyses indicated that these trends were due to different failure mechanisms for the two types of bone and the different loading modes. Taken together, these results provide unique insight into the failure behavior of human trabecular bone and highlight the need for a multiaxial failure criterion that accounts for anatomic site and bone volume fraction.     

Comparison of trabecular bone behavior in core and whole bone samples using high-resolution modeling of a vertebral body

Computational analysis of trabecular bone normally involves the modeling of (experimental tests of) cored samples. However, the lack of constraint on the sides of the extracted trabecular bone samples limits the information that can be inferred regarding true in situ behavior. Here, the element-by-element voxel-based finite element method was applied via, a custom-written software suite (FEEBE), to a 72 μm resolution model of an ovine vertebra. The difference between the apparent modulus of eight concentric core cylinders when modeled as part of the whole bone (containing 84 × 10⁶ degrees of freedom) and independent of the whole bone was investigated. The results showed that cored trabecular bone apparent modulus depended significantly on the core diameter when modeled as an extracted core (r ² = 0.975) and as part of a whole bone (r ² = 0.986). The cause of this result was separated into the side-artifact effect and bone volume fraction (BV/TV) effect. For the independently modeled cores, the apparent modulus of an inner core region of interest varied with increasing thickness of the outer annulus. This was attributed to the side-artifact effect, given that the BV/TV of the core region was constant. Within the whole trabecular structure, the side artifact was eliminated as the entire bone structure was modeled. However, a BV/TV effect influenced the apparent modulus depending on the size of the core selected for determining apparent modulus. Changing the size of the core varied the overall BV/TV of the core, and this significantly (r ² = 0.999) influences the apparent modulus. Therefore, determining a ‘true’ apparent modulus for trabecular bone was not achievable. The independently modeled cores consistently under-predict the in vivo apparent modulus. It is recommended that if a ‘true’ apparent modulus is required, the BV/TV at which it is required needs to be first determined. Apparent modeling of entire bones at microscale resolution allowed regions of low and high tissue strains to be identified, consistent with patterns of trabecular bone remodeling and resorption reported in literature. The basivertebral vein cavity underwent the highest strains within the entire vertebral body, suggesting that failure might initiate here, despite containing visibly thicker struts and plate trabeculae. Although computationally expensive, analysis of the entire vertebral body provided a full picture of in situ trabecular bone deformation.     

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.     

Intrinsic mechanical properties of trabecular calcaneus determined by finite-element models using 3D synchrotron microtomography

To determine intrinsic mechanical properties (elastic and failure) of trabecular calcaneus bone, chosen as a good predictor of hip fracture, we looked for the influence of image's size on a numerical simulation. One cubic sample of cancellous bone (9 x 9 x 9 mm(3)) was removed from the body of the calcaneus (6 females, 6 males, 79+/-9 yr). These samples were tested under compressive loading. Before compressive testing, these samples were imaged at 10.13 microm resolution using a 3D microcomputed tomography (muCT) (ESRF, France). The muCT images were converted to finite-element models. Depending on the bone density values (BV/TV), we compared two different finite element models: a linear hexahedral and a linear beam finite element models. Apparent experimental Young's modulus (E(app)(exp)) and maximum apparent experimental compressive stress (sigma(max)(exp)) were significantly correlated with bone density obtained by Archimedes's test (E(app)(exp)=236+/-231 MPa [19-742 MPa], sigma(max)(exp)=2.61+/-1.97 MPa [0.28-5.81 MPa], r>0.80, p<0.001). Under threshold at 40 microm, the size of the numerical samples (5.18(3) and 6.68(3)mm(3)) seems to be an important parameter on the accuracy of the results. The numerical trabecular Young's modulus was widely higher (E(trabecular)(num)=34,182+/-22,830 MPa [9700-87,211 MPa]) for the larger numerical samples and high BV/TV than those found classically by other techniques (4700-15,000 MPa). For rod-like bone samples (BV/TV<12%, n=7), Young's modulus, using linear beam element (E(trabecular)(num-skeleton): 10,305+/-5500 MPa), were closer to the Young's modulus found by other techniques. Those results show the limitation of hexahedral finite elements at 40 microm, mostly used, for thin trabecular structures.     

Gender specific LRP5 influences on trabecular bone structure and strength

A mutation in LRP5 (low-density lipoprotein receptor-related protein 5) has been shown to increase bone mass and density in humans and animals. Transgenic mice expressing the LRP5 mutation (G171V) demonstrate an increase in bone mass as compared to non-transgenic (NTG) littermates. This study evaluated LRP5 gene and gender-related influences on the structural and biomechanical strength properties of trabecular and cortical bone in femurs and vertebrae (L5) of 17-week-old mice. Micro-computed tomography was used to evaluate the trabecular bone structure of distal femurs and vertebrae ex vivo. Mechanical testing of the trabecular bone in the distal femur was done to determine biomechanical strength. Differences due to genotype and gender were tested using two-way ANOVA at a significance level of p<0.05. Trabecular bone structural parameters (BV/TV, trabecular thickness, number, etc.) at the distal femur, femoral neck, and vertebral body sites were greater in the transgenic as compared to the NTG mice. In addition, vertebral cortical thickness and trabecular strength parameters (ultimate and yield loads, stiffness, ultimate and yield stresses) in the distal femur were greater in the transgenic mice as compared to NTG. The increasing trends of cortical thickness were also noted in the transgenic mice as compared to NTG. Within LRP5 (G171V) mutant mice, there were significant gender-related differences in some of the trabecular bone structural parameters at all the sites (distal femur, femoral neck, and vertebral body). However, unlike trabecular structural parameters, the gender-specific differences were not found in the trabecular strength of LRP5 transgenic mice. In summary, these findings suggest that the LRP5 (G171V) mutation results in greater trabecular bone structure and strength at both the distal femurs and vertebral bodies as compared to NTG. In addition, only the trabecular structure parameters were affected by gender within the LRP5 (G171V) mutation.     

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.     

Volume to density relation in adult human bone tissue

Uniformity of tissue mineralisation is a strongly debated issue, due to its relation with bone mechanical behaviour. Bone mineral density (BMD) is measured in the clinical practice and is applied in computational application to derive material proprieties of bone tissue. However, BMD cannot identify if the variation in bone density is related to a modification of tissue mineral density (TMD), a change in bone volume or a combination of the two. This study was aimed to investigate whether TMD can be assumed as a constant in adult human bone (trabecular and cortical).A total number of 115 cylindrical bone specimens were collected. An inter-site analysis (96 specimens, 2 donors) was performed on cortical and trabecular specimens extracted from different anatomical sites. An intra-site study (19 specimens, 19 donors) was performed on specimens extracted from femoral heads. Bone volume fraction (BV/TV) was computed by means of a micro-computed tomography. Furthermore, ash density (ρ_ash) was measured. TMD was computed as the ratio between ρ_ash and BV/TV.It was found that the TMD of trabecular (1.24±0.16 g/cm³) and cortical (1.19±0.06 g/cm³) bone were not statistically different (p=0.31). Furthermore, the linear regression between ρ_ash and BV/TV was statistically significant (r²=0.99, p<0.001). Intra- and inter-site analyses demonstrated that the mineral distribution was independent of the extraction site.The present study suggests that TMD can be assumed reasonably constant in non-pathological adult bone tissue. Consequently, it is suggested that TMD can be managed as a constant in computational models, varying only BV in relation to clinical densitometric analysis.     

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.     

Dependence of mechanical compressive strength on local variations in microarchitecture in cancellous bone of proximal human femur

Human cancellous bone is a heterogeneous material. Despite this, most of the published studies report correlations between mechanical properties and morphometric parameters averaged on the whole specimen. This work investigated whether local variations in morphometric parameters were linked to the localized failure regions of cancellous bone. Additionally, it was examined whether local values of morphometric parameters can predict the ultimate stress better than the average bone volume fraction (BV/TV). Cylindrical cancellous bone specimens extracted along the primary compressive group of human femoral heads were studied. These were microCT-imaged to assess the morphometric parameters, compressed to determine the ultimate stress, and rescanned by microCT to visualize the failure region. Failure involved slightly less than half of the free height of the specimens. Significant differences were found in the morphometric parameters calculated in the failure and in the non-failure regions. The cross-sections containing minimum BV/TV values were those most often located inside the failure region (83%, p<0.001). Regression analysis confirmed that variations in BV/TV best describe variations in ultimate stress (R2=0.84) out of the averaged morphometric parameters. The prediction of ultimate stress increased when minimum or maximum values of the morphometric parameters were taken, with the highest prediction found by considering the minimum BV/TV (R2=0.95). In conclusion, due to the heterogeneity of cancellous bone, there may exist regions characterized by a different microarchitecture, where the bone is weaker and consequently is more likely to fail. These regions mostly contain minimum values in BV/TV, which were found to predict ultimate stress better than average BV/TV.     

The role of fabric in the large strain compressive behavior of human trabecular bone

Osteoporosis-related vertebral body fractures involve large compressive strains of trabecular bone. The small strain mechanical properties of the trabecular bone such as the elastic modulus or ultimate strength can be estimated using the volume fraction and a second order fabric tensor, but it remains unclear if similar estimations may be extended to large strain properties. Accordingly, the aim of this work is to identify the role of volume fraction and especially fabric in the large strain compressive behavior of human trabecular bone from various anatomical locations. Trabecular bone biopsies were extracted from human T12 vertebrae (n=31), distal radii (n=43), femoral head (n=44), and calcanei (n=30), scanned using microcomputed tomography to quantify bone volume fraction (BV/TV) and the fabric tensor (M), and tested either in unconfined or confined compression up to very large strains (～70%). The mechanical parameters of the resulting stress-strain curves were analyzed using regression models to examine the respective influence of BV/TV and fabric eigenvalues. The compressive stress-strain curves demonstrated linear elasticity, yielding with hardening up to an ultimate stress, softening toward a minimum stress, and a steady rehardening followed by a rapid densification. For the pooled experiments, the average minimum stress was 1.89 ± 1.77 MPa, while the corresponding mean strain was 7.15 ± 1.84%. The minimum stress showed a weaker dependence with fabric as the elastic modulus or ultimate strength. For the confined experiments, the stress at a logarithmic strain of 1.2 was 8.08 ± 7.91 MPa, and the dissipated energy density was 5.67 ± 4.42 MPa. The latter variable was strongly related to the volume fraction (R(2)=0.83) but the correlation improved only marginally with the inclusion of fabric (R(2)=0.84). The influence of fabric on the mechanical properties of human trabecular bone decreases with increasing strain, while the role of volume fraction remains important. In particular, the ratio of the minimum versus the maximum stress, i.e., the relative amount of softening, decreases strongly with fabric, while the dissipated energy density is dominated by the volume fraction. The collected results will prove to be useful for modeling the softening and densification of the trabecular bone using the finite element method.     

Damage in trabecular bone at small strains

Evidence that damage decreases bone quality, increases fracture susceptibility, and serves as a remodeling stimulus motivates further study of what loading magnitudes induce damage in trabecular bone. In particular, whether damage occurs at the smaller strains characteristic of habitual, as opposed to traumatic, loading is not known. The overall goal of this study was to characterize damage accumulation in trabecular bone at small strains (0.20 - 0.45% strain). A continuum damage mechanics approach was taken whereby damage was quantified by changes in modulus and residual strain. Human vertebral specimens (n = 7) were tested in compression using a multi-cycle load - unload protocol in which the maximum applied strain for each cycle, epsilonmax, was increased incrementally from epsilonmax = 0.20% on the first loading cycle to epsilonmax = 0.45% on the last cycle. Modulus and residual strain were measured for each cycle. Both changes in modulus and residual strains commenced at small strains, beginning as early as 0.24 and 0.20% strain, respectively. Strong correlations between changes in modulus and residual strains were observed (r = 0.51 - 0.98). Fully nonlinear, high-resolution finite element analyses indicated that even at small apparent strains, tissue-level strains were sufficiently high to cause local yielding. These results demonstrate that damage in trabecular bone occurs at apparent strains less than half the apparent compressive yield strain reported previously for human vertebral trabecular bone. Further, these findings imply that, as a consequence of the highly porous trabecular structure, tissue yielding can initiate at very low apparent strains and that this local failure has detectable and negative consequences on the apparent mechanical properties of trabecular bone.     

Normal masticatory function partially protects the rat mandibular bone from estrogen-deficiency induced osteoporosis

Background/aim

In a previous study we showed that mandibular alveolar (trabecular) bone appears to be less sensitive to estrogen deficiency than the proximal tibia spongiosa. We hypothesized that the mechanical loading of the alveolar process during mastication may protect the alveolar bone from the detrimental effects observed in other skeletal sites. To test this hypothesis we compared the effect of ovariectomy on the mandibular alveolar bone and the proximal tibia spongiosa of rats fed either a normal (hard) or a soft diet.

Methods

Forty six-month-old female Sprague–Dawley rats underwent trans-abdominal ovariectomy (OVX) or sham operation (SHAM). Half of the animals received their food in the usual form of pellets (hard consistency), while the other half received a soft, porridge-like, isocaloric diet of identical composition (soft consistency). Micro-computed tomographic histomorphometry was used to evaluate the trabecular micro-architecture. A two-factor analysis of variance was used to test for effects and interaction of ovariectomy and/or soft diet.

Results

OVX had a significantly negative effect on the proximal tibia spongiosa (all parameters under study except trabecular thickness; p<0.001) and on the mandibular alveolar bone (trabecular number and spacing; p<0.05). Soft diet led to a further decrease of mandibular BV/TV (p<0.01), trabecular thickness (p<0.05) and number (p<0.05), as well as increase of separation (p<0.001). A significant interaction was observed between OVX and soft diet concerning the mandibular BV/TV, as well as trabecular thickness and spacing (p<0.05).

Conclusion

Normal (hard) diet limited significantly the negative effects of estrogen deficiency on mandibular alveolar bone micro-architecture four months after ovariectomy.     

Trabecular bone response to mechanical loading in ovariectomized Sprague-Dawley rats depends on baseline bone quantity

Mechanical loading is one of the determining factors for bone modulation, and is therefore frequently used to treat or prevent bone loss; however, there appears to be no data on the effects of baseline bone quantity on this response. This study aimed to verify whether baseline bone quantity affects osteoporotic trabecular bone adaptive response to mechanical stimulation. Twenty-four female Sprague-Dawley (SD) rats were ovariectomized (OVX). After 3 weeks of OVX, rats were divided into a high bone quantity and a low bone quantity group, and rats in each group were then subdivided into 4 groups that were exposed to different loading strategies. In the loading groups, tibiae were stimulated through axial loading at 2000με of strain, for 1500 cycles each of 75s, 150s, or 250s. The sham treatment groups received no loading. Changes in BV/TV for trabecular bone in the tibia were measured at the baseline (before loading), and at 3 weeks and 6 weeks after loading. BV/TVs in loading groups of the low baseline bone quantity group were significantly increased at 6 weeks, compared with those in the no-loading groups (p<0.05), while those in the high quantity groups were not increased (p>0.05). A significant negative correlation was observed between baseline BV/TV and its relative variations at 3 weeks or 6 weeks (p<0.05). These results indicate that adaptive responses of osteoporotic trabecular bone to mechanical loading depend on baseline bone quantity.     

Combination of topological parameters and bone volume fraction better predicts the mechanical properties of trabecular bone

Trabecular bone structure may complement bone volume/total volume fraction (BV/TV) in the prediction of the mechanical properties. Nonetheless, the direct in vivo use of information pertaining to trabecular bone structure necessitates some predictive analytical model linking structure measures to mechanical properties. In this context, the purpose of this study was to combine BV/TV and topological parameters so as to better estimate the mechanical properties of trabecular bone. Thirteen trabecular bone mid-sagittal sections were imaged by magnetic resonance (MR) imaging at the resolution of 117 x 117x 300 microm(3). Topological parameters were evaluated in applying the 3D-line skeleton graph analysis (LSGA) technique to the binary MR images. The same images were used to estimate the elastic moduli by finite element analysis (FEA). In addition to the mid-sagittal section, two cylindrical samples were cored from each vertebra along vertical and horizontal directions. Monotonic compression tests were applied to these samples to measure both vertical and horizontal ultimate stresses. BV/TV was found as a strong predictor of the mechanical properties, accounting for 89-94% of the variability of the elastic moduli and for 69-86% of the variability of the ultimate stresses. Topological parameters and BV/TV were combined following two analytical formulations, based on: (1) the normalization of the topological parameters; and on (2) an exponential fit-model. The normalized parameters accounted for 96-98% of the variability of the elastic moduli, and the exponential model accounted for 80-95% of the variability of the ultimate stresses. Such formulations could potentially be used to increase the prediction of the mechanical properties of trabecular bone.     

Trabecular shear stress in human vertebral cancellous bone: intra- and inter-individual variations

Correlation of the mean and standard deviation of trabecular stresses has been proposed as a mechanism by which a strong relationship between the apparent strength and stiffness of cancellous bone can be achieved. The current study examined whether the relationship between the mean and standard deviation of trabecular von Mises stresses can be generalized for any group of cancellous bone. Cylindrical human vertebral cancellous bone specimens were cut in the infero-superior direction from T12 of 23 individuals (inter-individual group). Thirty nine additional specimens were prepared similarly from the T4-T12 and L2-L5 vertebrae of a 63 year old male (intra-individual group). The specimens were scanned by micro-computed tomography (microCT) and trabecular von Mises stresses were calculated using finite element modeling. The expected value, standard deviation and coefficient of variation of the von Mises stress were calculated form a three-parameter Weibull function fitted to von Mises stress data from each specimen. It was found that the average and standard deviation of trabecular von Mises shear stress were: (i) correlated with each other, supporting the idea that high correlation between the apparent strength and stiffness of cancellous bone can be achieved through controlling the trabecular level shear stress variations, (ii) dependent on anatomical site and sample group, suggesting that the variation of stresses are correlated to the mean stress to different degrees between vertebrae and individuals, and (iii) dependent on bone volume fraction, consistent with the idea that shear stress is less well controlled in bones with low BV/TV. The conversion of infero-superior loading into trabecular von Mises stresses was maximum for the tissue at the junction of the thoracic and lumbar spine (T12-L1) consistent with this junction being a common site of vertebral fracture.     

Application of homogenization theory to the study of trabecular bone mechanics.

It is generally accepted that the strength and stiffness of trabecular bone is strongly affected by trabecular microstructure. It has also been hypothesized that stress induced adaptation of trabecular bone is affected by trabecular tissue level stress and/or strain. At this time, however, there is no generally accepted (or easily accomplished) technique for predicting the effect of microstructure on trabecular bone apparent stiffness and strength or estimating tissue level stress or strain. In this paper, a recently developed mechanics theory specifically designed to analyze microstructured materials, called the homogenization theory, is presented and applied to analyze trabecular bone mechanics. Using the homogenization theory it is possible to perform microstructural and continuum analyses separately and then combine them in a systematic manner. Stiffness predictions from two different microstructural models of trabecular bone show reasonable agreement with experimental results, depending on metaphyseal region, (R2 greater than 0.5 for proximal humerus specimens, R2 less than 0.5 for distal femur and proximal tibia specimens). Estimates of both microstructural strain energy density (SED) and apparent SED show that there are large differences (up to 30 times) between apparent SED (as calculated by standard continuum finite element analyses) and the maximum microstructural or tissue SED. Furthermore, a strut and spherical void microstructure gave very different estimates of maximum tissue SED for the same bone volume fraction (BV/TV). The estimates from the spherical void microstructure are between 2 and 20 times greater than the strut microstructure at 10-20% BV/TV.     

The modified super-ellipsoid yield criterion for human trabecular bone

Despite the importance of multiaxial failure of trabecular bone in many biomechanical applications, to date no complete multiaxial failure criterion for human trabecular bone has been developed. By using experimentally validated nonlinear high-resolution, micromechanical finite-element models as a surrogate for multiaxial loading experiments, we determined the three-dimensional normal strain yield surface and all combinations of the two-dimensional normal-shear strain yield envelope. High-resolution finite-element models of three human femoral neck trabecular bone specimens obtained through microcomputed tomography were used. In total, 889 multiaxial-loading cases were analyzed, requiring over 41,000 CPU hours on parallel supercomputers. Our results indicated that the multiaxial yield behavior of trabecular bone in strain space was homogeneous across the specimens and nearly isotropic. Analysis of stress-strain curves along each axis in the 3-D normal strain space indicated uncoupled yield behavior whereas substantial coupling was seen for normal-shear loading. A modified super-ellipsoid surface with only four parameters fit the normal strain yield data very well with an arithmetic error +/-SD less than -0.04 +/- 5.1%. Furthermore, the principal strains associated with normal-shear loading showed excellent agreement with the yield surface obtained for normal strain loading (arithmetic error +/- SD < 2.5 +/- 6.5%). We conclude that the four-parameter "Modified Super-Ellipsoid" yield surface presented here describes the multiaxial failure behavior of human femoral neck trabecular bone very well.     

Similarity in the fatigue behavior of trabecular bone across site and species

Within the context of improving knowledge of the structure-function relations for trabecular bone for cyclic loading, we hypothesized that the S-N curve for cyclic compressive loading of trabecular bone, after accounting for differences in monotonic strength behavior, does not depend on either site or species. Thirty-five cores of fresh-frozen elderly human vertebral trabecular bone, harvested from nine donors (mean+/-S.D., age=74+/-17 years), were biomechanically tested in compression at sigma/E(0) values (ratio of applied stress to pre-fatigue elastic modulus) ranging from 0.0026 to 0.0070, and compared against literature data (J. Biomech. Eng. 120 (1998) 647-654) for young bovine tibial trabecular bone (n=37). As reported for the bovine bone, the number of cycles to failure for the human vertebral bone was related to sigma/E(0) by a power-law relation (r(2)=0.54, n=35). Quantitative comparison of these data against those reported for the bovine bone supported our hypothesis. Namely, when the differences in mean monotonic yield strain between the two types of bone were accounted for, a single S-N curve worked well for the pooled data (r(2)=0.75, n=72). Since elderly human vertebral and young bovine tibial trabecular bone represent two very different types of trabecular bone in terms of volume fraction and architecture, these findings suggest that the dominant failure mechanisms in trabecular bone for cyclic loading occur at the ultrastructural level.