Simulation of vertebral trabecular bone loss using voxel finite element analysis

Trabecular bone loss in human vertebral bone is characterised by thinning and eventual perforation of the horizontal trabeculae. Concurrently, vertical trabeculae are completely lost with no histological evidence of significant thinning. Such bone loss results in deterioration in apparent modulus and strength of the trabecular core. In this study, a voxel-based finite element program was used to model bone loss in three specimens of human vertebral trabecular bone. Three sets of analyses were completed. In Set 1, strain adaptive resorption was modelled, whereby elements which were subject to the lowest mechanical stimulus (principal strain) were removed. In Set 2, both strain adaptive and microdamage mechanisms of bone resorption were included. Perforation of vertical trabeculae occurred due to microdamage resorption of elements with strains that exceeded a damage threshold. This resulted in collapse of the trabecular network under compression loading for two of the specimens tested. In Set 3, the damage threshold strain was gradually increased as bone loss progressed, resulting in reduced levels of microdamage resorption. This mechanism resulted in trabecular architectures in which vertical trabeculae had been perforated and which exhibited similar apparent modulus properties compared to experimental values reported in the literature. Our results indicate that strain adaptive remodelling alone does not explain the deterioration in mechanical properties that have been observed experimentally. Our results also support the hypothesis that horizontal trabeculae are lost principally by strain adaptive resorption, while vertical trabeculae may be lost due to perforation from microdamage resorption followed by rapid strain adaptive resorption of the remaining unloaded trabeculae.     

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.     

A novel stress distribution organ in the arthropod exoskeleton

The vertebrate endoskeleton possesses a massive internal network of load-distributing trabeculae that in most locations accounts for the vast majority of bone cross sectional area. In contrast, arthropods rely on the external cuticle and its intermittent outpocketings to distribute the daily stresses of physiological loading. One of the constraints of the arthropod exoskeleton is the necessity to house the musculature involved in locomotion, feeding and etc. Because of this lack of an extensive internal load-distributing trabecular network, any load-distributing mechanism in arthropods would necessarily have to incorporate the exoskeleton. Several authors have identified structural apophysi whose functions presumably have mechanical significance, but few have been identified using quantitative analyses. This study investigates a novel stress-reducing structure arising from the articulation sites in the exoskeleton of the blue crab, Callinectes sapidus. During dissection of the merus-carpus joint and leg cuticle of the blue crab, an unique system of internal strut-like members was found radiating, both longitudinally and laterally, from the articular surface of the proximal merus segment, tapering into the diaphyseal region. This strut system, an internal outpocketing of the exoskeleton and semi-circular in cross section, mirrors the trabecular pattern seen radiating from vertebrate joint surfaces. Earlier reports of this structural system described it as a muscle attachment site and made little or no reference to potential load distribution properties. Finite element analysis (FEA) models confirm the efficacy of stress distributing properties of this articular strut system in the blue crab leg. In the models, the struts significantly reduce stress concentrations, reduce localized strains and minimize the risk of failure via buckling. Models lacking this strut system generate 94.7% larger peak von Mises stress at the articulation site, 37% higher peak displacement and 4% greater equivalent strain. The model with the struts is capable of withstanding an applied physiological load of up to 16.6 N prior to buckling, more than twice that of the model without struts (7.8 N). We suggest that this novel arthropod strut system is likely utilized at many joint surfaces at locations of high skeletal stress concentrations, is an adaptation for minimizing skeletal failure via localized buckling, and may be present in other arthropod taxa.     

A 3D elastic micropolar model of vertebral trabecular bone from lattice homogenization of the bone microstructure

A 3D anisotropic micropolar continuum model of vertebral trabecular bone is presently developed accounting for the influence of microstructure-related scale effects on the macroscopic effective properties. Vertebral trabecular bone is modeled as a cellular material with an idealized periodic structure made of open 3D cells. The micromechanical approach relies on the discrete homogenization technique considering lattice microrotations as additional degrees of freedom at the microscale. The effective elastic properties of 3D lattices made of articulated beams taking into account axial, transverse shearing, flexural, and torsional deformations of the cell struts are derived as closed form expressions of the geometrical and mechanical microparameters. The scaling laws of the effective moduli versus density are determined in situations of low and high effective densities to assess the impact of the transverse shear deformation. The classical and micropolar effective moduli and the internal flexural and torsional lengths are identified versus the same microparameters. A finite element model of the local architecture of the trabeculae gives values of the effective moduli that are in satisfactory agreement with the homogenized moduli.     

Comparison of the linear finite element prediction of deformation and strain of human cancellous bone to 3D digital volume correlation measurements

The mechanical properties of cancellous bone and the biological response of the tissue to mechanical loading are related to deformation and strain in the trabeculae during function. Due to the small size of trabeculae, their motion is difficult to measure. To avoid the need to measure trabecular motions during loading the finite element method has been used to estimate trabecular level mechanical deformation. This analytical approach has been empirically successful in that the analytical models are solvable and their results correlate with the macroscopically measured stiffness and strength of bones. The present work is a direct comparison of finite element predictions to measurements of the deformation and strain at near trabecular level. Using the method of digital volume correlation, we measured the deformation and calculated the strain at a resolution approaching the trabecular level for cancellous bone specimens loaded in uniaxial compression. Smoothed results from linearly elastic finite element models of the same mechanical tests were correlated to the empirical three-dimensional (3D) deformation in the direction of loading with a coefficient of determination as high as 97% and a slope of the prediction near one. However, real deformations in the directions perpendicular to the loading direction were not as well predicted by the analytical models. Our results show, that the finite element modeling of the internal deformation and strain in cancellous bone can be accurate in one direction but that this does not ensure accuracy for all deformations and strains.     

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.     

Physical and mechanical properties of calf lumbosacral trabecular bone.

The physical and mechanical properties of calf lumbar and sacral trabecular bone were determined and compared with those of human trabecular bone. The mean tissue density (1.66 +/- 0.12 g cm-3), equivalent mineral density (169 +/- 36 mg cm-3), apparent density (453 +/- 89 mg cm-3), ash density (194 +/- 59 mg cm-3), ash content (0.6 +/- 0.05%), compressive strength (7.1 +/- 3.0 MPa) and compressive modulus (173 +/- 97 MPa) of calf trabecular bone are similar to those of young human. There were moderate, positive linear correlations between apparent density and equivalent mineral density, ash density, and compressive strength; and between compressive strength and equivalent mineral density (R2 ranging from 0.35 to 0.48, p less than 0.001). Apparent density, ash density, and equivalent mineral density did not differ significantly in different regions. In contrast to humans, the compressive strength increased from posterior, near the facet, to the anterior vertebral body. These comparisons of physical and mechanical properties, as well as anatomical comparisons by others, indicate that the calf spine is a good model of the young non-osteoporotic human spine and thus useful for the testing of spinal instrumentation.     

The mechanical properties of trabecular bone: Dependence on anatomic location and function

In 1961, Evans and King documented the mechanical properties of trabecular bone from multiple locations in the proximal human femur. Since this time, many investigators have cataloged the distribution of trabecular bone material properties from multiple locations within the human skeleton to include femur, tibia, humerus, radius, vertebral bodies, and iliac crest. The results of these studies have revealed tremendous variations in material properties and anisotropy. These variations have been attributed to functional remodeling as dictated by Wolff's Law. Both linear and power functions have been found to explain the relationship between trabecular bone density and material properties. Recent studies have re-emphasized the need to accurately quantify trabecular bone architecture proposing several algorithms capable of determining the anisotropy, connectivity and morphology of the bone. These past studies, as well as continuing work, have significantly increased the accuracy of analytical and experimental models investigating bone, and bone/implant interfaces as well as enhanced our perspective towards understanding the factors which may influence bone formation or resorption.     

Finite element modeling of trabecular bone damage

This paper presents a finite element-based, computational model for analysis of structural damage to trabecular bone tissues. A modulus reduction method was formulated from elasto-plasticity theory, and was used to account for site-specific trabecular bone tissue damage. Trabecular bone tissue damage is illustrated using a large-scale, anatomically accurate, two-dimensional, microstructural finite element model of a human thoracic vertebral body. Four models with varying specifications for damage accumulation were subjected to compressive loading and unloading cycles. The numerical results and experimental validation demonstrated that the modulus reduction method reproduced the non-linear mechanical behaviour of vertebal trabecular bone. The iterative computational approach presented provides a methodology to study trabecular bone damage, and should provide researchers with a computational approach to study bone fracture and repair and to predict vertebral fragility.     

Examination of continuum and micro-structural properties of human vertebral cancellous bone using combined cellular solid models

Two- and three-dimensional structural models of the vertebral body have been used to estimate the mechanical importance of parameters that are difficult to quantify experimentally such as lattice disorder, trabecular thickness, trabecular spacing, connectivity, and fabric. Many of the models that investigate structure–function relationships of the vertebral body focus only on the trabecular architecture and neglect solid–fluid interactions. We developed a cellular solid model composed of two idealized unit cell geometries to investigate the continuum and micro-structural properties of human vertebral cancellous bone in a mathematically tractable model. Using existing histomorphological data we developed structure–function relationships for the mechanical properties of the solid phase, estimated the micro-structural strains, and predicted the fluid flow characteristics. We found that the micro-structural strains may be 1.7 to 2.2 times higher than the continuum level strains between the ages of 40 and 80. In addition, the predicted permeability agrees well with the experimental data.     

Trabecular bone structure and strength - remodelling and repair

The strength of the spinal trabecular bone declines by a factor of 4-5 from the age of 20 to 80 years. At the same time, the volumetric (apparent) density declines by a factor of only 2. This discrepancy can be explained by the known power relationship between density and strength; this power relationship is based on the fact that trabecular bone is a porous material. To date, it has not been possible to determine or quantify the influence other factors may have in determining the strength of a loadbearing trabecular network. However, it is known that with age: 1) There is a loss of connectivity through osteoclastic perforations of horizontal struts. 2) There is an increase in anisotropy - again due to loss of horizontal struts, and perhaps also due to micro-modelling drift or to thickening of some vertical trabeculae. 3) The changes in the network can lead to the slenderness ratio between vertical and horizontal struts reaching a certain magnitude and thereby inducing buckling under compression. 4) Microdamage and microfractures will occur - mainly in these very loaded vertical struts. The microfractures will be repaired by microcallus formation, and these calluses will later be removed by the remodelling process. 5) Bone material quality will slightly change, leading to a decrease in collagen content and a relative increase in the degree of mineralisation. But, it is not known how these factors will influence the power relationship between density and strength. Nor is it known how different treatment regimens will affect the 'natural' power relationship: will the same curve be followed, but in the opposite direction? Or will the curve be less or more steep? Will the gain in bone strength be larger if treatment is started early - on the steep part of the curve? Furthermore, as trabecular bone can never be isolated in vivo, other factors need to be investigated: The interplay between the cortical shell and the trabecular network; transmission of load; the interplay between soft tissues (cartilage, connective tissue, muscle) and bone; the shock absorbing capacity of the discs; and the hydraulic effect of the bone marrow. In order to answer these questions, more in vitro and in vivo studies on human bone in relation to aging, to immobilisation, to exercise and in relation to different treatment regimens are needed.     

Analysis of crack growth in a 3D Voronoi structure: a model for fatigue in low density trabecular bone

Both creep and crack growth contribute to the reduction in modulus associated with fatigue loading in bone. Here we simulate crack growth and subsequent strut failure in fatigue in an open-cell, three-dimensional Voronoi structure which is similar to that of low density, osteoporotic bone. The model indicates that sequential failure of struts leads to a precipitous drop in modulus: the failure of 1% of the struts leads to about a 10% decrease in modulus. A parametric study is performed to assess the influence of normalized stress range, relative density, initial crack size, crack shape and cell geometry on the fatigue life. The fatigue life is most sensitive to the relative density and the initial crack length. The results lead to a quantitative expression for the fatigue life associated with crack growth. Data for the fatigue life of trabecular bone are compared with the crack growth model described in this paper as well as with a previous model for creep of a three-dimensional Voronoi structure. In our models, creep dominates the fatigue behavior in low cycle fatigue while crack growth dominates in high cycle fatigue, consistent with previous observations on cortical bone. The large scatter in the trabecular bone fatigue data make it impossible to identify a transition between creep dominated fatigue and crack growth dominated fatigue. The parametric study of the crack growth model indicates that variations in relative density among specimens, initial crack size within trabeculae and crack shape could easily produce such variability in the test results.     

Failure modelling of trabecular bone using a non-linear combined damage and fracture voxel finite element approach

Trabecular bone tissue failure can be considered as consisting of two stages: damage and fracture; however, most failure analyses of 3D high-resolution trabecular bone samples are confined to damage mechanisms only, that is, without fracture. This study aims to develop a computational model of trabecular bone consisting of an explicit representation of complete failure, incorporating damage criteria, fracture criteria, cohesive forces, asymmetry and large deformation capabilities. Following parameter studies on a test specimen, and experimental testing of bone sample to complete failure, the asymmetric critical tissue damage and fracture strains of ovine vertebral trabecular bone were calibrated and validated to be compression damage ?1.16 %, tension damage 0.69 %, compression fracture ?2.91 % and tension fracture 1.98 %. Ultimate strength and post–ultimate strength softening were captured by the computational model, and the failure of individual struts in bending and shear was also predicted. This modelling approach incorporated a cohesive parameter that provided a facility to calibrate ductile–brittle behaviour of bone tissue in this non-linear geometric and non-linear constitutive property analyses tool. Finally, the full accumulation of tissue damage and tissue fracture has been monitored from range of small magnitude (normal daily loading) through to specimen yielding, ultimate strength and post–ultimate strength softening.     

Mechanical loading of mouse caudal vertebrae increases trabecular and cortical bone mass-dependence on dose and genotype

Most in vivo studies addressing the skeletal responses of mice to mechanical loading have targeted cortical bone. To investigate trabecular bone responses also we have developed a caudal vertebral axial compression device (CVAD) that transmits mechanical loads to compress the fifth caudal vertebra via stainless steel pins inserted into the forth and sixth caudal vertebral bodies. Here, we used the CVAD in C57BL/6 (B6) and C3H/Hej (C3H) female mice (15 weeks of age) to investigate whether the effect of regular bouts of mechanical stimulation on bone modeling and bone mass was dependent on dose and genotype. A combined micro-computed tomographic and dynamic histomorphometric analysis was carried out at the end of a 4-week loading regimen (3,000 cycles, 10 Hz, 3 x week) for load amplitudes of 0N, 2N, 4N and 8N. Significant increases in trabecular bone mass of 9 and 21% for loads of 4N and 8N, respectively, were observed in B6 mice. A significant increase of 10% in trabecular bone mass occurred for a load of 8N in the C3H strain. For other loads, no significant increases were detected. Both mouse strains exhibited substantial increases in trabecular bone formation rates for all loads, B6: 111% (2N), 86% (4N), 164% (8N), C3H: 41% (2N), 38% (4N), 141% (8N). Significant decreases in osteoclast number of 146 and 93% for a load of 8N were detected in B6 and C3H mice, respectively. These findings demonstrate that the effect of loading on the structural and functional parameters of bone is dose and genotype dependent. The caudal vertebral loading model established here is proposed for further studies addressing the molecular processes involved in the skeletal responses to mechanical stimuli.     

Dependence of mechanical properties of trabecular bone on plate–rod microstructure determined by individual trabecula segmentation (ITS)

Individual trabecula segmentation (ITS) technique can decompose the trabecular bone network into individual trabecular plates and rods and is capable of quantifying the plate/rod-related microstructural characteristics of trabecular bone. This novel technique has been shown to be able to provide in-depth insights into micromechanics and failure mechanisms of human trabecular bone, as well as to distinguish the fracture status independent of area bone mineral density in clinical applications. However, the plate/rod microstructural parameters from ITS have never been correlated to experimentally determined mechanical properties of human trabecular bone. In this study, on-axis cylindrical trabecular bone samples from human proximal tibia (n=22), vertebral body (n=10), and proximal femur (n=21) were harvested, prepared, scanned using micro computed-tomography (µCT), analyzed with ITS and mechanically tested. Regression analyses showed that the plate bone volume fraction (pBV/TV) and axial bone volume fraction (aBV/TV) calculated by ITS analysis correlated the best with elastic modulus (R²=0.96–0.97) and yield strength (R²=0.95–0.96). Trabecular plate-related microstructural parameters correlated highly with elastic modulus and yield strength, while most rod-related parameters were found inversely and only moderately correlated with the mechanical properties. In addition, ITS analysis also identified that trabecular bone at human femoral neck had the highest trabecular plate-related parameters while the other sites were similar with each other in terms of plate–rod microstructure.     

Mechanics of the taper integrated screwed-in (TIS) abutments used in dental implants

The tapered implant-abutment interface is becoming more popular due to the mechanical reliability of retention it provides. Consequently, understanding the mechanical properties of the tapered interface with or without a screw at the bottom has been the subject of a considerable amount of studies involving experiments and finite element (FE) analysis. This paper focuses on the tapered implant-abutment interface with a screw integrated at the bottom of the abutment. The tightening and loosening torques are the main factors in determining the reliability and the stability of the attachment. Analytical formulas are developed to predict tightening and loosening torque values by combining the equations related to the tapered interface with screw mechanics equations. This enables the identification of the effects of the parameters such as friction, geometric properties of the screw, the taper angle, and the elastic properties of the materials on the mechanics of the system. In particular, a relation between the tightening torque and the screw pretension is identified. It was shown that the loosening torque is smaller than the tightening torque for typical values of the parameters. Most of the tightening load is carried by the tapered section of the abutment, and in certain combinations of the parameters the pretension in the screw may become zero. The calculations performed to determine the loosening torque as a percentage of tightening torque resulted in the range 85-137%, depending on the values of taper angle and the friction coefficient.     

The effects of tensile-compressive loading mode and microarchitecture on microdamage in human vertebral cancellous bone

The amount of microdamage in bone tissue impairs mechanical performance and may act as a stimulus for bone remodeling. Here we determine how loading mode (tension vs. compression) and microstructure (trabecular microarchitecture, local trabecular thickness, and presence of resorption cavities) influence the number and volume of microdamage sites generated in cancellous bone following a single overload. Twenty paired cylindrical specimens of human vertebral cancellous bone from 10 donors (47–78 years) were mechanically loaded to apparent yield in either compression or tension, and imaged in three dimensions for microarchitecture and microdamage (voxel size 0.7×0.7×5.0 μm³). We found that the overall proportion of damaged tissue was greater (p=0.01) for apparent tension loading (3.9±2.4%, mean±SD) than for apparent compression loading (1.9±1.3%). Individual microdamage sites generated in tension were larger in volume (p<0.001) but not more numerous (p=0.64) than sites in compression. For both loading modes, the proportion of damaged tissue varied more across donors than with bone volume fraction, traditional measures of microarchitecture (trabecular thickness, trabecular separation, etc.), apparent Young?s modulus, or strength. Microdamage tended to occur in regions of greater trabecular thickness but not near observable resorption cavities. Taken together, these findings indicate that, regardless of loading mode, accumulation of microdamage in cancellous bone after monotonic loading to yield is influenced by donor characteristics other than traditional measures of microarchitecture, suggesting a possible role for tissue material properties.     

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.     

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.