Multiaxial strength characteristics of trabecular bone

Bovine trabecular bone specimens were tested in multiaxial stress, including pure shear, in a special test device. Shear strength was proportional to apparent density to the exponent 1.65, in approximate agreement with theoretical and experimental studies on the shear strengths of porous foams. The mean shear strength was 6.60 +/- 1.66 MPa, after normalizing for apparent density. This compares well with normalized shear strengths from Saha and Gorman's (1981) study on human femora. A scanning electron microscope study indicated random trabecular architecture and a complex fracture mechanism at the level of the individual trabecular struts. Hoffmann's (1967) 3-D isotropic failure criterion was applied to the multiaxial test data, along with data from uniaxial compression tests, indicating a compressive strength approximately three times the tensile strength.     

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.     

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.     

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.     

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.     

Mechanical properties of metaphyseal bone in the proximal femur

We used a three-point bending test to investigate the structural behavior of 123 rectangular flat plate specimens harvested from the metaphyseal shell of the cervical and intertrochanteric regions of five fresh/frozen human proximal femora. For comparison purposes, 36 specimens of similar geometry were also fabricated from bone of the femoral diaphysis. All specimens were oriented in either the local longitudinal or transverse directions. The mean longitudinal elastic modulus was 9650 +/- 2410 (SD) MPa and demonstrated a 24% decrease from that measured for the diaphysis (12500 +/- 2140 MPa) using the same testing technique. However, the transverse elastic moduli did not differ significantly between the proximal (5470 +/- 1720 MPa) and diaphyseal (5990 +/- 1520 MPa) specimens. The globally averaged values for the ultimate tensile strengths of the metaphyseal shell were 101 +/- 26 MPa in the longitudinal and 50 +/- 12 MPa in the transverse directions. These compared with diaphyseal values of 128 +/- 16 MPa and 47 +/- 12 MPa, respectively. While these differences were largely due to the reduced density of the proximal specimens, a slight decrease in transverse anisotropy for the proximal specimens was also noted by comparing the ratio of longitudinal to transverse moduli (1.76) and tensile strength (2.02) to the diaphyseal values (2.09 and 2.71, respectively). Use of these data should lead to improved performance of analytical models for the proximal femur, and thus help focus increased attention on the structural contribution of trabecular bone to the strength and rigidity of the proximal femur.     

The effects of damage and microcracking on the impact strength of bone

Microcracking has been shown to occur when bone is 'damaged' as shown by a loss of stiffness. The effect on bone's toughness of the types of damage produced at low losses of stiffness are not known. We loaded bovine bone specimens in bending and tension to stiffness losses of up to 27%, and examined the microcracking produced. The tensile specimens had diffuse arrays of microcracks of 2-20 microm in length, characteristic of tensile loading, on all surfaces. The bending specimens showed tensile microcracking on the tensile surface and characteristic long, straight, cross-hatched compression cracks on the compressive surface. Specimens were then broken in impact. Those that had been damaged in bending were divided into two groups, in one group the part of the specimen which had undergone compression damage was placed in tension, and in the other group the tensile damage was placed in tension. Tensile damage loaded in tension did not reduce the bone's energy-absorbing ability in impact until a modulus reduction of over 20%. However compression damage loaded in tension did severely reduce the bone's energy absorption capabilities (by an average of about 40%).     

Application of the Tsai-Wu quadratic multiaxial failure criterion to bovine trabecular bone

As a first step toward development of a multiaxial failure criterion for human trabecular bone, the Tsai-Wu quadratic failure criterion was modified as a function of apparent density and applied to bovine tibial trabecular bone. Previous data from uniaxial compressive, tensile, and torsion tests (n = 139 total) were combined with those from new triaxial tests (n = 17) to calibrate and then verify the criterion. Combinations of axial compression and radial pressure were used to produce the triaxial compressive stress states. All tests were performed with minimal end artifacts in the principal material coordinate system of the trabecular network. Results indicated that the stress interaction term F12 exhibited a strong nonlinear dependence on apparent density (r2 > 0.99), ranging from -0.126 MPa-2 at low densities (0.29 g/cm3) to 0.005 MPa-2 at high densities (0.63 g/cm3). After calibration and when used to predict behavior of new-specimens without any curve-fitting, the Tsai-Wu criterion had a mean (+/- SD) error of -32.6 +/- 10.6 percent. Except for the highest density triaxial specimens, most (15/17 specimens) failed at axial stresses close to their predicted uniaxial values, and some reinforcement for transverse loading was observed. We conclude that the Tsai-Wu quadratic criterion, as formulated here, is at best only a reasonable predictor of the multiaxial failure behavior of trabecular bone, and further work is required before it can be confidently applied to human bone.     

Heterogeneity of the mechanical properties of demineralized bone

Knowledge of the mechanical properties of the collagenous component of bone is required for composite modeling of bone tissue and for understanding the age- and disease-related reductions in the ductility and strength of bone. The overall goal of this study was to investigate the heterogeneity of the mechanical properties of demineralized bone which remains unexplained and may be due to differences in the collagen structure or organization or in experimental protocols. Uniaxial tension tests were conducted to measure the elastic and failure properties of demineralized human femoral (n = 10) and tibial (n = 13) and bovine humeral (n = 8) and tibial (n = 8) cortical bone. Elastic modulus differed between groups (p = 0.02), varying from 275 +/- 94 MPa (mean +/- SD) to 450 + 50 MPa. Similarly, ultimate stress varied across groups from 15 + 4.2 to 26 + 4.7 MPa (p = 0.03). No significant differences in strain-to-failure were observed between any groups in this study (pooled mean of 8.4 +/- 1.6%; p = 0.42). However, Bowman et al. (1996) reported an average ultimate strain of 12.3 +/- 0.5% for demineralized bovine humeral bone, nearly 40% higher than our value. Taken together, it follows that all the monotonic mechanical properties of demineralized bone can display substantial heterogeneity. Future studies directed at explaining such differences may therefore provide insight into aging and disease of bone tissue.     

Tensile and compressive properties of cancellous bone

The relationship between the mechanical properties of trabecular bone in tension and compression was investigated by non-destructive testing of the same specimens in tension and compression, followed by random allocation to a destructive test in either tension or compression. There was no difference between Young's modulus in tension and compression, and there was a strong positive correlation between the values (R = 0.97). Strength, ultimate strain and work to failure was significantly higher in tensile testing than in compressive testing.     

Trabecular bone microdamage and microstructural stresses under uniaxial compression

The balance between local remodeling and accumulation of trabecular bone microdamage is believed to play an important role in the maintenance of skeletal integrity. However, the local mechanical parameters associated with microdamage initiation are not well understood. Using histological damage labeling, micro-CT imaging, and image-based finite element analysis, regions of trabecular bone microdamage were detected and registered to estimated microstructural von Mises effective stresses and strains, maximum principal stresses and strains, and strain energy density (SED). Bovine tibial trabecular bone cores underwent a stepwise uniaxial compression routine in which specimens were micro-CT imaged following each compression step. The results indicate that the mode of trabecular failure observed by micro-CT imaging agreed well with the polarity and distribution of stresses within an individual trabecula. Analysis of on-axis subsections within specimens provided significant positive relationships between microdamage and each estimated tissue stress, strain and SED parameter. In a more localized analysis, individual microdamaged and undamaged trabeculae were extracted from specimens loaded within the elastic region and to the apparent yield point. As expected, damaged trabeculae in both groups possessed significantly higher local stresses and strains than undamaged trabeculae. The results also indicated that microdamage initiation occurred prior to apparent yield at local principal stresses in the range of 88-121 MPa for compression and 35-43 MPa for tension and local principal strains of 0.46-0.63% in compression and 0.18-0.24% in tension. These data provide an important step towards understanding factors contributing to microdamage initiation and establishing local failure criteria for normal and diseased trabecular bone.     

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.     

Bone creep-fatigue damage accumulation

Creep and fatigue tests were performed on human femoral cortical bone and the results were compared to a cumulative damage model for bone fracture. Fatigue tests in tension, compression, and reversed loading with a tensile mean stress were conducted at 2 Hz and 0.02 Hz. Load frequency had a strong influence on the number of cycles to failure but did not influence the total time to failure. Bone displayed poor creep-fracture properties in both tension and compression. The fracture surfaces of the tensile creep specimens are distinctly different than those of the compressive specimens. The results suggest that tensile cyclic loading creates primarily time-dependent damage and compressive cyclic loading creates primarily cycle-dependent damage. However, data for load histories involving both tensile and compressive loading indicate lower time to failure than predicted by a simple summation of time-dependent and cycle-dependent damage.     

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.     

Axial-shear interaction effects on microdamage in bovine tibial trabecular bone

Because many osteoporotic fractures occur during a fall, understanding the effect of off-axis loads on initiation and propagation of microdamage in trabecular bone should provide further insight into the biomechanics of age-related fractures. Fourteen on-axis cylindrical specimens were prepared from 12 bovine tibiae. Fluorescent stains were used to label the microdamage due to a sequence of compressive and torsional damaging loads. The mean decrease in Young's modulus was over four times greater than that in the shear modulus after the compressive overload, while there was no difference between the decrease in the axial and torsional stiffnesses after the torsional overload. The total microcrack density due to compression was uniform across the radius of the cylindrical specimens, while the mean density of microcracks due to torsional overloading increased from the axis of the cylindrical specimen to the circumference. The high density of microcracks near the axis of the specimen following torsional overloading was unexpected because of the low strains. Nearly 40% of the microcracks due to torsion propagated from pre-existing microcracks caused by axial compression, indicating that existing microcracks may extend at relatively low strain if the loading mode changes. The propagating microcracks were, on average, longer than the initiating microcracks due to either compressive or torsional loading. Damage due to axial compression appears to increase the susceptibility of trabecular bone to damage propagation during subsequent torsional loads, but it has little effect on the elastic properties in shear.     

A cellular solid criterion for predicting the axial-shear failure properties of bovine trabecular bone.

In a long-term effort to develop a complete multi-axial failure criterion for human trabecular bone, the overall goal of this study was to compare the ability of a simple cellular solid mechanistic criterion versus the Tsai-Wu, Principal Strain, and von Mises phenomenological criteria--all normalized to minimize effects of interspecimen heterogeneity of strength--to predict the on-axis axial-shear failure properties of bovine trabecular bone. The Cellular Solid criterion that was developed here assumed that vertical trabeculae failed due to a linear superposition of axial compression/tension and bending stresses, induced by the apparent level axial and shear loading, respectively. Twenty-seven bovine tibial trabecular bone specimens were destructively tested on-axis without end artifacts, loaded either in combined tension-torsion (n = 10), compression-torsion (n = 11), or uniaxially (n = 6). For compression-shear, the mean (+/- S.D.) percentage errors between measured values and criterion predictions were 7.7 +/- 12.6 percent, 19.7 +/- 23.2 percent, 22.8 +/- 18.9 percent, and 82.4 +/- 64.5 percent for the Cellular Solid, Tsai-Wu, Principal Strain, and von Mises criteria, respectively; corresponding mean errors for tension-shear were -5.2 +/- 11.8 percent, 14.3 +/- 12.5 percent, 6.9 +/- 7.6 percent, and 57.7 +/- 46.3 percent. Statistical analysis indicated that the Cellular Solid criterion was the best performer for compression-shear, and performed as well as the Principal Strain criterion for tension-shear. These data should substantially improve the ability to predict axial-shear failure of dense trabecular bone. More importantly, the results firmly establish the importance of cellular solid analysis for understanding and predicting the multiaxial failure behavior of trabecular 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.     

Mechanical testing of cancellous bone from the femoral head: experimental errors due to off-axis measurements

The aim of this study was to verify whether a misalignment between the testing direction and the trabecular main direction has a significant effect on the compressive behaviour of cancellous bone. Ten cylindrical specimens were extracted from femoral heads with a misalignment to the trabecular main direction of approximately 20 degrees. Each specimen was paired with a specimen extracted aligned with the main direction of the trabeculae on the basis of the closest bone volume fraction, obtaining two groups, one 'aligned' and one 'misaligned'. The average off-axis angle was 6.1 degrees and 21.6 degrees for the 'aligned' and 'misaligned' group, respectively. The specimens underwent micro-tomographic analysis, compressive testing, micro-indentation testing and ashing. No significant differences were found in histomorphometric parameters, hardness and ash density between the two groups, whereas significant differences were found in Young's modulus and ultimate stress: both parameters, measured for the 'misaligned' group, were about 40% lower than those measured for the 'aligned' group. These results demonstrate a great effect of the angle between the testing direction and the main direction of the trabecular structure on the compressive behaviour of cancellous bone. This angle should be reduced as much as possible (in the present work the average value was 6.6+/-3.3 degrees), in any case measured, and always reported together with the mechanical parameters of cancellous bone.     

Nonlinear tensile properties of bovine articular cartilage and their variation with age and depth

Tensile stiffness of articular cartilage is much greater than its compressive stiffness and plays an essential role even in compressive properties by increasing transient fluid pressures during physiological loading. Recent studies of nonlinear properties of articular cartilage in compression revealed several physiologically pertinent nonlinear behaviors, all of which required that cartilage tensile stiffness increase significantly with stretch. We therefore performed sequences of uniaxial tension tests on fresh bovine articular cartilage slices using a protocol that allowed several hours to attain equilibrium and measured longitudinal and transverse tissue strain. By testing bovine cartilage from different ages (6 months to 6 years) we found that equilibrium and transient tensile modulus increased significantly with maturation and age, from 0 to 15 MPa at equilibrium and from 10 to 28 MPa transiently. Our results indicate that cartilage stiffens with age in a manner similar to other highly hydrated connective tissues, possibly due to age-dependent content of enzymatic and nonenzymatic collagen cross links. The long relaxation period used in our tests (5-10 hours) was necessary in order to attain equilibrium and avoid a very significant overestimation of equilibrium modulus that occurs when much shorter times are used (15-30 minutes). We also found that equilibrium and transient tensile modulus increased nonlinearly when cartilage is stretched from 0 to 10% strain without any previous tare load. Although our results estimate a nonlinear increase in tensile stiffness with stretch that is an order of magnitude lower than that required to predict nonlinear properties in compression, they are in agreement with previous results from other uniaxial tension tests of collagenous materials. We therefore speculate that biaxial tensile moduli may be much higher and thereby more compatible with observed nonlinear compressive properties.     

Indirect determination of trabecular bone effective tissue failure properties using micro-finite element simulations

Trabecular bone strength is marked not only by the onset of local yielding, but also by post-yield behavior. To study and predict trabecular bone elastic and yield properties, micro-finite element (micro-FE) models were successfully applied. However, trabecular bone strength predictions require micro-FE models incorporating post-yield behavior of trabecular bone tissue. Due to experimental difficulties, such data is currently not available. Here we used micro-FE modeling to determine failure behavior of trabecular bone tissue indirectly, by iteratively fitting FE simulation to experimental results. Failure parameters were fitted to an isotropic plasticity model based on Hill's yield function, using materially and geometrically nonlinear micro-FE models of seven bovine trabecular bone specimens. The predictive value of the averaged effective tissue properties was subsequently tested. The results showed that compression softening had to be included on the tissue level in order to accurately describe the apparent-level behavior of the bone specimens. A sensitivity study revealed that the simulated response was less sensitive to variations in the post-yield properties of the bone tissue than variations in the elastic and yield properties. Due to fitting of the tissue properties, apparent-level behavior could be accurately reproduced for each specimen separately. Predictions based on the averaged and fixed tissue properties were less accurate, due to inter-specimen variations in the tissue properties.