The sensitivity of nonlinear computational models of trabecular bone to tissue level constitutive model

Microarchitectural finite element models have become a key tool in the analysis of trabecular bone. Robust, accurate, and validated constitutive models would enhance confidence in predictive applications of these models and in their usefulness as accurate assays of tissue properties. Human trabecular bone specimens from the femoral neck (n = 3), greater trochanter (n = 6), and lumbar vertebra (n = 1) of eight different donors were scanned by μ-CT and converted to voxel-based finite element models. Unconfined uniaxial compression and shear loading were simulated for each of three different constitutive models: a principal strain-based model, Drucker–Lode, and Drucker–Prager. The latter was applied with both infinitesimal and finite kinematics. Apparent yield strains exhibited minimal dependence on the constitutive model, differing by at most 16.1%, with the kinematic formulation being influential in compression loading. At the tissue level, the quantities and locations of yielded tissue were insensitive to the constitutive model, with the exception of the Drucker–Lode model, suggesting that correlation of microdamage with computational models does not improve the ability to discriminate between constitutive laws. Taken together, it is unlikely that a tissue constitutive model can be fully validated from apparent-level experiments alone, as the calculations are too insensitive to identify differences in the outcomes. Rather, any asymmetric criterion with a valid yield surface will likely be suitable for most trabecular bone models.     

On the interrelationship of permeability and structural parameters of vertebral trabecular bone: a parametric computational study

To characterise the flow of a fluid through a uniform porous medium, the medium may be completely described by its permeability, a measure of flow resistance. Fluid flow in the intertrabecular spaces of cancellous bone has been recognised as an important factor in a number of physical phenomena. In order to investigate the interdependence of permeability, porosity and the structural parameters, we adapted a morphological model and systematically varied its structural parameters. By simulating a viscous Stokes flow regime, we were able to estimate the anisotropic permeability tensor and performed an extensive, stepwise multivariate regression analysis to establish empirical relationships between the morphological parameters and the permeability for the anatomical directions individually. The regression analysis indicated high values of determination coefficients [0.88 < R² < 0.89 (transversal directions) and R² = 0.60 (longitudinal direction), porosity-based prediction and R² = 0.98 for all directions and information presented to the regression model]. We conclude that a pooled set of structural parameters may explain up to 98% of the permeability variability, the regression model predicts permeability values that match experimental data, and a good prediction performance could be achieved by only incorporating the porosity and either the degree of anisotropy (0.89 < R² < 0.91) or the trabecular spacing predictor (0.96 < R² < 0.97). These conclusions imply that trabecular thickness and shape parameters only play a minor role in the determination of vertebral trabecular bone permeability. However, a major limitation of the model is that it reflects an idealisation of the real, regionally varying structure of trabecular bone. Therefore, the goodness-of-fit estimates we are presenting should be considered as an upper bound limitation regarding the prediction performance.     

Validation of composite finite elements efficiently simulating elasticity of trabecular bone

Patient-specific analyses of the mechanical properties of bones become increasingly important for the management of patients with osteoporosis. The potential of composite finite elements (CFEs), a novel FE technique, to assess the apparent stiffness of vertebral trabecular bone is investigated in this study. Segmented volumes of cylindrical specimens of trabecular bone are compared to measured volumes. Elasticity under uniaxial loading conditions is simulated; apparent stiffnesses are compared to experimentally determined values. Computational efficiency is assessed and recommendations for simulation parameters are given. Validating apparent uniaxial stiffnesses results in concordance correlation coefficients 0.69 ≤ r_𝒸 ≤ 0.92 for resolutions finer than 168 μm, and an average error of 5.8% between experimental and numerical results at 24 μm resolution. As an application, the code was used to compute local, macroscopic stiffness tensors for the trabecular structure of a lumbar vertebra. The presented technique allows for computing stiffness using smooth FE meshes at resolutions that are well achievable in peripheral high resolution quantitative CT. Therefore, CFEs could be a valuable tool for the patient-specific assessment of bone stiffness.     

Fluid–structure interactions in micro-interlocked regions of the cement–bone interface

Experimental tests and computational modelling were used to explore the fluid dynamics at the trabeculae–cement interlock regions found in the tibial component of total knee replacements. A cement–bone construct of the proximal tibia was created to simulate the immediate post-operative condition. Gap distributions along nine trabeculae–cement regions ranged from 0 to 50.4 μm (mean = 12 μm). Micro-motions ranged from 0.56 to 4.7 μm with a 1 MPa compressive load to the cement. Fluid–structure analysis between the trabeculae and the cement used idealised models with parametric evaluation of loading direction, gap closing fraction (GCF), gap thickness, loading frequency and fluid viscosity. The highest fluid shear stresses (926 Pa) along the trabecular surface were found for conditions with very thin and large GCFs, much larger than reported physiological levels (～1–5 Pa). A second fluid–structure model was created with a provision for bone resorption using a constitutive model with resorption velocity proportional to fluid shear rate. A lower cut-off was used, below which bone resorption would not occur (50 s^? 1). Results showed that there was initially high shear rates (>1000 s^? 1) that diminished after initial trabecular resorption. Resorption continued in high shear rate regions, resulting in a final shape with bone left deep in the cement layer, and is consistent with morphology found in post-mortem retrievals. Small gaps between the trabecular surface and the cement in the immediate post-operative state produce fluid flow conditions that appear to be supra-physiologic; these may cause fluid-induced lysis of trabeculae in the micro-interlock regions.     

Influence of anisotropic bone properties on the biomechanical behavior of the acetabular cup implant: a multiscale finite element study

Although the biomechanical behavior of the acetabular cup (AC) implant is determinant for the surgical success, it remains difficult to be assessed due to the multiscale and anisotropic nature of bone tissue. The aim of the present study was to investigate the influence of the anisotropic properties of peri-implant trabecular bone tissue on the biomechanical behavior of the AC implant at the macroscopic scale. Thirteen bovine trabecular bone samples were imaged using micro-computed tomography (μCT) with a resolution of 18 μm. The anisotropic biomechanical properties of each sample were determined at the scale of the centimeter based on a dedicated method using asymptotic homogenization. The material properties obtained with this multiscale approach were used as input data in a 3D finite element model to simulate the macroscopic mechanical behavior of the AC implant under different loading conditions. The largest stress and strain magnitudes were found around the equatorial rim and in the polar area of the AC implant. All macroscopic stiffness quantities were significantly correlated (R² > 0.85, p < 6.5 e-6) with BV/TV (bone volume/total volume). Moreover, the maximum value of the von Mises stress field was significantly correlated with BV/TV (R² > 0.61, p < 1.6 e-3) and was always found at the bone-implant interface. However, the mean value of the microscopic stress (at the scale of the trabeculae) decrease as a function of BV/TV for vertical and torsional loading and do not depend on BV/TV for horizontal loading. These results highlight the importance of the anisotropic properties of bone tissue.     

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.     

Bone and Gastric Bypass Surgery: Effects of Dietary Calcium and Vitamin D

Objective:To examine bone mass and metabolism in women who had previously undergone Roux‐en‐Y gastric bypass (RYGB) and determine the effect of supplementation with calcium (Ca) and vitamin D. Research Methods and Procedures: Bone mineral density and bone mineral content (BMC) were examined in 44 RYGB women (≥3 years post‐surgery; 31% weight loss; BMI, 34 kg/m²) and compared with age‐ and weight‐matched control (CNT) women (n = 65). In a separate analysis, RYGB women who presented with low bone mass (n = 13) were supplemented to a total 1.2 g Ca/d and 8 μg vitamin D/d over 6 months and compared with an unsupplemented CNT group (n = 13). Bone mass and turnover and serum parathyroid hormone (PTH) and 25‐hydroxyvitamin D were measured. Results:Bone mass did not differ between premenopausal RYGB and CNT women (42 ± 5 years), whereas postmenopausal RYGB women (55 ± 7 years) had higher bone mineral density and BMC at the lumbar spine and lower BMC at the femoral neck. Before and after dietary supplementation, bone mass was similar, and serum PTH and markers of bone resorption were higher (p < 0.001) in RYGB compared with CNT women and did not change significantly after supplementation. Discussion: Postmenopausal RYGB women show evidence of secondary hyperparathyroidism, elevated bone resorption, and patterns of bone loss (reduced femoral neck and higher lumbar spine) similar to other subjects with hyperparathyroidism. Although a modest increase in Ca or vitamin D does not suppress PTH or bone resorption, it is possible that greater dietary supplementation may be beneficial.     

The influence of dietary sodium on bone development in growing rats

The present study investigated the effects of dietary sodium on bone growth in young rats. Five-week-old rats were fed one of three different diets for 60 days: low sodium (NaCl, 0.32 g/kg diet), normal sodium (NaCl, 2.6 g/kg) and high sodium (NaCl, 20 g/kg). The proximal tibial metaphysis (PTM), the fifth lumbar vertebra (LV₅) and the middle part of the tibia shaft (TX) were analysed by bone histomorphometry. The expression of three osteogenesis genes, Runx2, osteopontin and osteocalcin, was determined by RT-PCR in bone samples from the skull. In both the PTM and LV₅, trabecular area and thickness were increased by the low-sodium diet, while the high-sodium diet decreased trabecular area in LV₅. Dynamic data revealed that sodium restriction increased bone formation parameters in the PTM and LV₅, but decreased bone resorption in LV₅. In TX, endosteal bone formation was enhanced by the low-sodium diet and depressed by the high-sodium diet compared to the normal sodium group. But there were no statistically changes in the cortical bone area of TX. Low-sodium intake significantly enhanced the expression of all three osteogenesis genes compared to the normal sodium group, while high-sodium intake suppressed osteogenic gene expression. Our results suggest that sodium restriction in growing rats promotes bone development by influencing both bone formation and resorption.     

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.     

Accounting for spatial variation of trabecular anisotropy with subject-specific finite element modeling moderately improves predictions of local subchondral bone stiffness at the proximal tibia

IntroductionPreviously, a finite element (FE) model of the proximal tibia was developed and validated against experimentally measured local subchondral stiffness. This model indicated modest predictions of stiffness (R² = 0.77, normalized root mean squared error (RMSE%) = 16.6%). Trabecular bone though was modeled with isotropic material properties despite its orthotropic anisotropy. The objective of this study was to identify the anisotropic FE modeling approach which best predicted (with largest explained variance and least amount of error) local subchondral bone stiffness at the proximal tibia.MethodsLocal stiffness was measured at the subchondral surface of 13 medial/lateral tibial compartments using in situ macro indentation testing. An FE model of each specimen was generated assuming uniform anisotropy with 14 different combinations of cortical- and tibial-specific density-modulus relationships taken from the literature. Two FE models of each specimen were also generated which accounted for the spatial variation of trabecular bone anisotropy directly from clinical CT images using grey-level structure tensor and Cowin’s fabric-elasticity equations. Stiffness was calculated using FE and compared to measured stiffness in terms of R² and RMSE%.ResultsThe uniform anisotropic FE model explained 53–74% of the measured stiffness variance, with RMSE% ranging from 12.4 to 245.3%. The models which accounted for spatial variation of trabecular bone anisotropy predicted 76–79% of the variance in stiffness with RMSE% being 11.2–11.5%.ConclusionsOf the 16 evaluated finite element models in this study, the combination of Synder and Schneider (for cortical bone) and Cowin’s fabric-elasticity equations (for trabecular bone) best predicted local subchondral bone stiffness.     

Hispidulin Exerts Anti-osteoporotic Activity in Ovariectomized Mice via Activating AMPK Signaling Pathway

To investigate the effect of Hispidulin on ovariectomy (OVX)-induced bone loss in mice. Female mice subjected to OVX were treated with Hispidulin for 8 weeks. The total body bone mineral density was measured at the beginning and after the OVX at a time interval of 4 weeks. Micro-computed tomography of the tibia, bone histomorphometric analysis of the femur, and biomechanical analysis of tibia, vertebra, and femoral head were performed to fully evaluate the anti-osteoporotic effect of Hispidulin. Western blot analysis was performed to determine the level of activated AMPK. Hispidulin treatment effectively prevented OVX-induced body weight loss and attenuated OVX-induced bone loss. Hispidulin treatment also decreased trabecular spacing in OVX mice. The suppressing effect of Hispidulin on osteoclast surface and number was also found via histomorphometric analysis. Western blot analysis revealed that Hispidulin significantly elevated the activated AMPK levels. Our findings suggest that Hispidulin exerts anti-osteoporotic and bone resorption attenuating effects via activating the AMPK signaling pathway.     

RANKL from bone marrow adipose lineage cells promotes osteoclast formation and bone loss

Receptor activator of NF‐κB ligand (RANKL) is essential for osteoclast formation and bone remodeling. Nevertheless, the cellular source of RANKL for osteoclastogenesis has not been fully uncovered. Different from peripheral adipose tissue, bone marrow (BM) adipose lineage cells originate from bone marrow mesenchymal stromal cells (BMSCs). Here, we demonstrate that adiponectin promoter‐driven Cre expression (Adipoq^Cre ) can target bone marrow adipose lineage cells. We cross the Adipoq^Cre mice with rankl^fl/fl mice to conditionally delete RANKL from BM adipose lineage cells. Conditional deletion of RANKL increases cancellous bone mass of long bones in mice by reducing the formation of trabecular osteoclasts and inhibiting bone resorption but does not affect cortical bone thickness or resorption of calcified cartilage. Adipoq^Cre; rankl^fl/fl mice exhibit resistance to estrogen deficiency and rosiglitazone (ROS)‐induced trabecular bone loss but show bone loss induced by unloading. BM adipose lineage cells therefore represent an essential source of RANKL for the formation of trabecula osteoclasts and resorption of cancellous bone during remodeling under physiological and pathological conditions. Targeting bone marrow adiposity is a promising way of preventing pathological bone loss.     

Interactive graph-cut segmentation for fast creation of finite element models from clinical ct data for hip fracture prediction

In this study, we propose interactive graph cut image segmentation for fast creation of femur finite element (FE) models from clinical computed tomography scans for hip fracture prediction. Using a sample of N = 48 bone scans representing normal, osteopenic and osteoporotic subjects, the proximal femur was segmented using manual (gold standard) and graph cut segmentation. Segmentations were subsequently used to generate FE models to calculate overall stiffness and peak force in a sideways fall simulations. Results show that, comparable FE results can be obtained with the graph cut method, with a reduction from 20 to 2–5 min interaction time. Average differences between segmentation methods of 0.22 mm were not significantly correlated with differences in FE derived stiffness (R² = 0.08, p = 0.05) and weakly correlated to differences in FE derived peak force (R² = 0.16, p = 0.01). We further found that changes in automatically assigned boundary conditions as a consequence of small segmentation differences were significantly correlated with FE derived results. The proposed interactive graph cut segmentation software MITK-GEM is freely available online at https://simtk.org/home/mitk-gem.     

Impaired bone formation and osteopenia in heterozygous β<Superscript>IVSII-654</Superscript> knockin thalassemic mice

β-thalassemia caused by the C→T mutation at nucleotide 654 of the intron 2 (β^IVSII-654) results in aberrant splicing of β-globin RNA, leading to an almost absence of β-globin synthesis. Although trabecular and cortical bone loss was previously reported in β-thalassemic mice with deletion of β-globin gene, the microscopic changes in trabecular structure in β^IVSII-654 thalassemic mice remained elusive. Here, we investigated the macroscopic and microscopic bone changes in 12-week-old β^IVSII-654 knockin thalassemic mice by dual-energy X-ray absorptiometry (DXA) and histomorphometric analysis, respectively. DXA revealed a decrease in bone mineral density in the lumbar vertebrae and tibial metaphysis, but not in the femoral diaphysis, suggesting that β^IVSII-654 thalassemia predominantly led to osteopenia at the trabecular site, but not the cortical site. Further histomorphometric analysis of the tibial secondary spongiosa showed that trabecular bone volume was significantly decreased with the expansion of marrow cavity. Decreases in osteoblast surface, osteoid surface, mineral apposition rate, mineralizing surface, and mineralized volume were also observed. Moreover, trabecular bone resorption was markedly enhanced as indicated by increases in the osteoclast surface and eroded surface. It could be concluded that β^IVSII-654 thalassemia impaired bone formation and enhanced bone resorption, thereby leading to osteopenia especially at the trabecular sites, such as the tibial metaphysis.     

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.     

Principal trabecular structural orientation predicted by quantitative ultrasound is strongly correlated with \upmu FEA determined anisotropic apparent stiffness

The microarchitecture and alignment of trabecular bone adapts to the particular mechanical milieu applied to it. Due to this anisotropic mechanical property, measurement orientation has to be taken into consideration when assessing trabecular bone quality and fracture risk prediction. Quantitative ultrasound (QUS) has demonstrated the ability in predicting the principal structural orientation (PSO) of trabecular bone. Although the QUS prediction for PSO is very close to that of \(\upmu \) CT, certain angle differences still exist. It remains unknown whether this angle difference can induce significant differences in mechanical properties or not. The objective of this study was to evaluate the mechanical properties in different PSOs predicted using different methods, QUS and \(\upmu \) CT, thus to investigate the ability of QUS as a means to predict the PSO of trabecular bone noninvasively. By validating the ability of QUS to predict the PSO of trabecular bone, it is beneficial for future QUS applications because QUS measurements in the PSO can provide information more correlated with the mechanical properties than with other orientations. In this study, seven trabecular bone balls from distal bovine femurs were used to generate finite element models based on the 3-dimensional \(\upmu \) CT images. Uniaxial compressive loading was performed on the bone ball models in the finite element analysis (FEA) in six different orientations (three anatomical orientations, two PSOs predicted by QUS and the longest vector of mean intercept length (MIL) tensor calculated by \(\upmu \) CT). The stiffness was calculated based on the reaction force of the bone balls under loading, and the von Mises stress results showed that both the mechanical properties in the PSOs predicted by QUS are significantly higher than the anatomical orientations and comparatively close to the longest vector of MIL tensor. The stiffness in the PSOs predicted by QUS is also highly correlated with the stiffness in the MIL tensor orientation (ATTmax vs. MIL, \(R^{2}\)  = 0.98, \(p<001\) ; UVmax vs. MIL, \(R^{2}\)  = 0.92, \(p<001\) ). These results were validated by in vitro mechanical testing on the bone ball samples. This study demonstrates that the PSO of trabecular bone predicted by QUS has an equally strong apparent stiffness with the orientation predicted by \(\upmu \) CT.     

Contribution to FE modeling for intraoperative pedicle screw strength prediction

Although the use of pedicle screws is considered safe, mechanical issues still often occur. Commonly reported issues are screw loosening, screw bending and screw fracture. The aim of this study was to develop a Finite Element (FE) model for the study of pedicle screw biomechanics and for the prediction of the intraoperative pullout strength. The model includes both a parameterized screw model and a patient-specific vertebra model. Pullout experiments were performed on 30 human cadaveric vertebrae from ten donors. The experimental force-displacement data served to evaluate the FE model performance. μCT images were taken before and after screw insertion, allowing the creation of an accurate 3D-model and a precise representation of the mechanical properties of the bone. The experimental results revealed a significant positive correlation between bone mineral density (BMD) and pullout strength (Spearman ρ = 0.59, p < 0.001) as well as between BMD and pullout stiffness (Spearman ρ = 0.59, p < 0.001). A high positive correlation was also found between the pullout strength and stiffness (Spearman ρ = 0.84, p < 0.0001). The FE model was able to reproduce the linear part of the experimental force-displacement curve. Moreover, a high positive correlation was found between numerical and experimental pullout stiffness (Pearson ρ = 0.96, p < 0.005) and strength (Pearson ρ = 0.90, p < 0.05). Once fully validated, this model opens the way for a detailed study of pedicle screw biomechanics and for future adjustments of the screw design.     

Finite element analysis of trabecular bone structure: a comparison of image-based meshing techniques

In this study, we investigate if finite element (FE) analyses of human trabecular bone architecture based on 168 μm images can provide relevant information about the bone mechanical characteristics. Three human trabecular bone samples, one taken from the femoral head, one from the iliac crest, and one from the lumbar spine, were imaged with micro-computed tomography (micro-CT) using a 28 μm resolution. After reconstruction the resolution was coarsened to 168 μm. First, all reconstructions were thresholded and directly converted to FE-models built of hexahedral elements. For the coarser resolutions of two samples, this resulted in a loss of trabecular connections and a subsequent loss of stiffness. To reduce this effect, a tetrahedral element meshing based on the marching cubes algorithm, as well as a modified hexahedron meshing, which thresholds the image such that load carrying bone mass is preserved, were employed. For each sample elastic moduli and tissue Von Mises stresses of the three different 168 μm models were compared to those from the hexahedron 28 μm model. For one sample the hexahedron meshing at 168 μm produced excellent results. For the other two samples the results obtained from the hexahedral models at 168 μm resolution were poor. Considerably better results were attained for these samples when using the mass-compensated or tetrahedron meshing techniques. We conclude that the accuracy of the FE-models at 168 μm strongly depends on the bone morphology, in particular its trabecular thickness. A substantial loss of trabecular connections during the hexahedron meshing process indicates that poor FE results will be obtained. In this case the tetrahedron or mass-compensated hexahedron meshing techniques can reduce the loss of connections and produce better results than the plain hexahedron meshing techniques.     

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.