Quantification of trabecular spatial orientation from low-resolution images

No accepted methodology exists to assess trabecular bone orientation from clinical CT scans. The aim of this study was to test the hypothesis that the distribution of grey values in clinical CT images is related to the underlying trabecular architecture and that this distribution can be used to identify the principal directions and local anisotropy of trabecular bone. Fourteen trabecular bone samples were extracted from high-resolution (30 μm) micro-CT scans of seven human femoral heads. Trabecular orientations and local anisotropy were calculated using grey-level deviation (GLD), a novel method providing a measure of the three-dimensional distribution of image grey values. This was repeated for different image resolutions down to 300 μm and for volumes of interest (VOIs) ranging from 1 to 7 mm. Outcomes were compared with the principal mechanical directions and with mean intercept length (MIL) as calculated for the segmented 30-μm images. For the 30-μm images, GLD predicted the mechanical principal directions equally well as MIL. For the 300-μm images, which are resolutions that can be obtained in vivo using clinical CT, only a small increase (3°–6°) in the deviation from the mechanical orientations was found. VOIs of 5 mm resulted in a robust quantification of the orientation. We conclude that GLD can quantify structural bone parameters from low-resolution CT images.     

The effect of resorption cavities on bone stiffness is site dependent

Resorption cavities formed during the bone remodelling cycle change the structure and thus the mechanical properties of trabecular bone. We tested the hypotheses that bone stiffness loss due to resorption cavities depends on anatomical location, and that for identical eroded bone volumes, cavities would cause more stiffness loss than homogeneous erosion. For this purpose, we used beam–shell finite element models. This new approach was validated against voxel-based FE models. We found an excellent agreement for the elastic stiffness behaviour of individual trabeculae in axial compression (R² = 1.00) and in bending (R²>0.98), as well as for entire trabecular bone samples to which resorption cavities were digitally added (R² = 0.96, RMSE = 5.2%). After validation, this new method was used to model discrete cavities, with dimensions taken from a statistical distribution, on a dataset of 120 trabecular bone samples from three anatomical sites (4th lumbar vertebra, femoral head, iliac crest). Resorption cavities led to significant reductions in bone stiffness. The largest stiffness loss was found for samples from the 4th lumbar vertebra, the lowest for femoral head samples. For all anatomical sites, resorption cavities caused significantly more stiffness loss than homogeneous erosion did. This novel technique can be used further to evaluate the impact of resorption cavities, which are known to change in several metabolic bone diseases and due to treatment, on bone competence.     

Estimation of the effective yield properties of human trabecular bone using nonlinear micro-finite element analyses

Micro-finite element (\(\upmu \)FE) analyses are often used to determine the apparent mechanical properties of trabecular bone volumes. Yet, these apparent properties depend strongly on the applied boundary conditions (BCs) for the limited size of volumes that can be obtained from human bones. To attenuate the influence of the BCs, we computed the yield properties of samples loaded via a surrounding layer of trabecular bone (“embedded configuration”). Thirteen cubic volumes (10.6 mm side length) were collected from \(\upmu \)CT reconstructions of human vertebrae and femora and converted into \(\upmu \)FE models. An isotropic elasto-plastic material model was chosen for bone tissue, and nonlinear \(\upmu \)FE analyses of six uniaxial, shear, and multi-axial load cases were simulated to determine the yield properties of a subregion (5.3 mm side length) of each volume. Three BCs were tested. Kinematic uniform BCs (KUBCs: each boundary node is constrained with uniform displacements) and periodicity-compatible mixed uniform BCs (PMUBCs: each boundary node is constrained with a uniform combination of displacements and tractions mimicking the periodic BCs for an orthotropic material) were directly applied to the subregions, while the embedded configuration was achieved by applying PMUBCs on the larger volumes instead. Yield stresses and strains, and element damage at yield were finally compared across BCs. Our findings indicate that yield strains do not depend on the BCs. However, KUBCs significantly overestimate yield stresses obtained in the embedded configuration (+43.1 ± 27.9%). PMUBCs underestimate (?10.0 ± 11.2%), but not significantly, yield stresses in the embedded situation. Similarly, KUBCs lead to higher damage levels than PMUBCs (+51.0 ± 16.9%) and embedded configurations (+48.4 ± 15.0%). PMUBCs are better suited for reproducing the loading conditions in subregions of the trabecular bone and deliver a fair estimation of their effective (asymptotic) yield properties.     

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.     

Validation of a voxel-based FE method for prediction of the uniaxial apparent modulus of human trabecular bone using macroscopic mechanical tests and nanoindentation

Assessment of the mechanical properties of trabecular bone is of major biological and clinical importance for the investigation of bone diseases, fractures and their treatments. Finite element (FE) methods are getting increasingly popular for quantifying the elastic and failure properties of trabecular bone. In particular, voxel-based FE methods have been previously used to calculate the effective elastic properties of trabecular microstructures. However, in most studies, bone tissue moduli were assumed or back-calculated to match the apparent elastic moduli from experiments, which often lead to surprisingly low values when compared to nanoindentation results. In this study, voxel-based FE analysis of trabecular bone is combined with physical measures of volume fraction, micro-CT (microCT) reconstructions, uniaxial mechanical tests and specimen-specific nanoindentation tests for proper validation of the method. Cylindrical specimens of cancellous bone were extracted from human femurs and their volume fraction determined with Archimede's method. Uniaxial apparent modulus of the specimens was measured with an improved tension-compression testing protocol that minimizes boundary artefacts. Their microCT reconstructions were segmented to match the measured bone volume fraction and used to create full-size voxel models with 30-45 microm element size. For each specimen, linear isotropic elastic material properties were defined based on specific nanoindentation measurements of its embedded bone tissue. Linear FE analyses were finally performed to simulate the uniaxial mechanical tests. Additional parametric analyses were performed to evaluate the potential errors on the predicted apparent modulus arising from variations in segmentation threshold, tissue modulus, and the use of 125-mm(3) cubic sub-regions. The results demonstrate an excellent correspondence between experimental measures and FE predictions of uniaxial apparent modulus. In conclusion, the adopted voxel-based FE approach is found to be a robust method to predict the linear elastic properties of human cancellous bone, provided segmentation of the microCT reconstructions is carefully calibrated, tissue modulus is known a priori and the entire region of interest is included in the analysis.     

A Three-Dimensional Finite Element Scheme to Investigate the Apparent Mechanical Properties of Trabecular Bone

A computational technique is described for investigating the apparent mechanical properties of trabecular bone based on tissue geometry obtained from the marching cubes volume rendering scheme. Using this scheme, a 3D representation of the trabecular bone was extracted from two-dimensional cross-sections of the tissue originating from a quantitative serial sectioning procedure. Surface information consists of node coordinates and polygon connectivity in a 3D space. A custom, adaptive mesh generation technique using a normal offset was used to prepare 3D finite element volume meshes (4-node tetrahedral elements) of variable mesh density from the extracted surface geometry. Nine target mesh resolutions (32 μm to 107 μm) were examined for a (1.5 mmx 1.5 mmx 2 mm) volume of trabecular bone. A mesh density of 50,000 elements/mm(3) of bone tissue was found to be adequate for convergence of apparent (bulk) modulus for 1% uniaxial compression. For this convergent case, the maximum local normal compressive tissue stress was 400 MPa which was six hundred-fold greater than the computed apparent stress. Variation in the apparent modulus was less than 5% when Poisson's ratio values were varied between 0.1 and 0.4. Poisson's ratio values greater than 0.4 had a more marked effect on the apparent modulus. Based upon these results, approximately 1 million, 4-node tetrahedral elements are required to analyze a continuum scale model of trabecular bone (5 mm cube).     

Towards validation of computational analyses of peri-implant displacements by means of experimentally obtained displacement maps

Micro-finite element (μFE) analysis has recently been introduced for the detailed quantification of the mechanical interaction between bone and implant. The technique has been validated at an apparent level. The aim of this study was to address the accuracy of μFE analysis at the trabecular level. Experimental displacement fields were obtained by deformable image registration, also known as strain mapping (SM), of dynamic hip screws implanted in three human femoral heads. In addition, displacement fields were calculated using μFE analysis. On a voxel-by-voxel basis, the coefficients of determination (R(2)) between experimental and μFE-calculated displacements ranged from 0.67 to 0.92. Linear regression of the mean displacements over nine volumes of interest yielded R(2) between 0.81 and 0.84. The lowest R(2) values were found in regions of very small displacements. In conclusion, we found that peri-implant bone displacements calculated with μFE analysis correlated well with displacements obtained from experimental SM.     

Micromechanical analyses of vertebral trabecular bone based on individual trabeculae segmentation of plates and rods

Trabecular plates play an important role in determining elastic moduli of trabecular bone. However, the relative contribution of trabecular plates and rods to strength behavior is still not clear. In this study, individual trabeculae segmentation (ITS) and nonlinear finite element (FE) analyses were used to evaluate the roles of trabecular types and orientations in the failure initiation and progression in human vertebral trabecular bone. Fifteen human vertebral trabecular bone samples were imaged using micro computed tomography (μCT), and segmented using ITS into individual plates and rods by orientation (longitudinal, oblique, and transverse). Nonlinear FE analysis was conducted to perform a compression simulation for each sample up to 1% apparent strain. The apparent and relative trabecular number and tissue fraction of failed trabecular plates and rods were recorded during loading and data were stratified by trabecular orientation. More trabecular rods (both in number and tissue fraction) failed at the initiation of compression (0.1–0.2% apparent strain) while more plates failed around the apparent yield point (>0.7% apparent strain). A significant correlation between plate bone volume fraction (pBV/TV) and apparent yield strength was found (r²=0.85). From 0.3% to 1% apparent strain, significantly more longitudinal trabecular plate and transverse rod failed than other types of trabeculae. While failure initiates at rods and rods fail disproportionally to their number, plates contribute significantly to the apparent yield strength because of their larger number and tissue volume. The relative failed number and tissue fraction at apparent yield point indicate homogeneous local failure in plates and rods of different orientations.     

Prediction of apparent trabecular bone stiffness through fourth-order fabric tensors

The apparent stiffness tensor is an important mechanical parameter for characterizing trabecular bone. Previous studies have modeled this parameter as a function of mechanical properties of the tissue, bone density, and a second-order fabric tensor, which encodes both anisotropy and orientation of trabecular bone. Although these models yield strong correlations between observed and predicted stiffness tensors, there is still space for reducing accuracy errors. In this paper, we propose a model that uses fourth-order instead of second-order fabric tensors. First, the totally symmetric part of the stiffness tensor is assumed proportional to the fourth-order fabric tensor in the logarithmic scale. Second, the asymmetric part of the stiffness tensor is derived from relationships among components of the harmonic tensor decomposition of the stiffness tensor. The mean intercept length (MIL), generalized MIL (GMIL), and fourth-order global structure tensor were computed from images acquired through microcomputed tomography of 264 specimens of the femur. The predicted tensors were compared to the stiffness tensors computed by using the micro-finite element method (\(\upmu \)FE), which was considered as the gold standard, yielding strong correlations (\(R^2\) above 0.962). The GMIL tensor yielded the best results among the tested fabric tensors. The Frobenius error, geodesic error, and the error of the norm were reduced by applying the proposed model by 3.75, 0.07, and 3.16 %, respectively, compared to the model by Zysset and Curnier (Mech Mater 21(4):243–250, 1995) with the second-order MIL tensor. From the results, fourth-order fabric tensors are a good alternative to the more expensive \(\upmu \)FE stiffness predictions.     

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.     

Stochastic multi-scale prediction on the apparent elastic moduli of trabecular bone considering uncertainties of biological apatite (BAp) crystallite orientation and image-based modelling

An assessment of the mechanical properties of trabecular bone is important in determining the fracture risk of human bones. Many uncertainty factors contribute to the dispersion of the estimated mechanical properties of trabecular bone. This study was undertaken in order to propose a computational scheme that will be able to predict the effective apparent elastic moduli of trabecular bone considering the uncertainties that are primarily caused by image-based modelling and trabecular stiffness orientation. The effect of image-based modelling which focused on the connectivity was also investigated. A stochastic multi-scale method using a first-order perturbation-based and asymptotic homogenisation theory was applied to formulate the stochastically apparent elastic properties of trabecular bone. The effective apparent elastic modulus was predicted with the introduction of a coefficient factor to represent the variation of bone characteristics due to inter-individual differences. The mean value of the predicted effective apparent Young's modulus in principal axis was found at approximately 460 MPa for respective 15.24% of bone volume fraction, and this is in good agreement with other experimental results. The proposed method may provide a reference for the reliable evaluation of the prediction of the apparent elastic properties of trabecular bone.     

The use of preserved tissue in finite element modelling of fresh ovine vertebral behaviour

The aim of this study was to investigate whether the predicted finite element (FE) stiffness of vertebral bone is altered when using images of preserved rather than fresh tissue to generate specimen-specific FE models. Fresh ovine vertebrae were used to represent embalmed (n = 3) and macerated dry-bone (n = 3) specimens and treated accordingly. Specimens were scanned pre- and post-treatment using micro-computed tomography. A constant threshold level derived from these images was used to calculate the respective bone volume fraction (BV/TV) from which the conversion factor validated for fresh tissue was used to determine material properties that were assigned to corresponding FE models. Results showed a definite change in the BV/TV between the fresh and the preserved bone. However, the changes in the predicted FE stiffness were not generally greater than the variations expected from assignment of loading and boundary conditions. In conclusion, images of preserved tissue can be used to generate FE models that are representative of fresh tissue with a tolerable level of error.     

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.     

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.     

Finite element calculated uniaxial apparent stiffness is a consistent predictor of uniaxial apparent strength in human vertebral cancellous bone tested with different boundary conditions

Strong correspondence between the uniaxial apparent strength and stiffness of cancellous bone allows the use of stiffness as a predictor of bone strength. Measured values of mechanical properties in cancellous bone can be different between experiments due to different experimental conditions. In the current study, bone volume fraction, experimentally determined and finite element (FE) predicted stiffness were examined as predictors of cancellous bone ultimate strength in two different groups each of which was tested using a different end constraint. It is demonstrated that, although always significant, the relationships of strength with bone volume fraction and experimentally determined stiffness are different between test groups. Apparent stiffness, estimated by FE modeling, predicts the ultimate strength of human cancellous bone consistently for all examined experimental protocols.     

Effects of suppression of bone turnover on cortical and trabecular load sharing in the canine vertebral body

The relative biomechanical effects of antiresorptive treatment on cortical thickness vs. trabecular bone microarchitecture in the spine are not well understood. To address this, T-10 vertebral bodies were analyzed from skeletally mature female beagle dogs that had been treated with oral saline (n=8 control) or a high dose of oral risedronate (0.5 mg/kg/day, n=9 RIS-suppressed) for 1 year. Two linearly elastic finite element models (36-μm voxel size) were generated for each vertebral body—a whole-vertebra model and a trabecular-compartment model—and subjected to uniform compressive loading. Tissue-level material properties were kept constant to isolate the effects of changes in microstructure alone. Suppression of bone turnover resulted in increased stiffness of the whole vertebra (20.9%, p=0.02) and the trabecular compartment (26.0%, p=0.01), while the computed stiffness of the cortical shell (difference between whole-vertebra and trabecular-compartment stiffnesses, 11.7%, p=0.15) was statistically unaltered. Regression analyses indicated subtle but significant changes in the relative structural roles of the cortical shell and the trabecular compartment. Despite higher average cortical shell thickness in RIS-suppressed vertebrae (23.1%, p=0.002), the maximum load taken by the shell for a given value of shell mass fraction was lower (p=0.005) for the RIS-suppressed group. Taken together, our results suggest that—in this canine model—the overall changes in the compressive stiffness of the vertebral body due to suppression of bone turnover were attributable more to the changes in the trabecular compartment than in the cortical shell. Such biomechanical studies provide an unique insight into higher-scale effects such as the biomechanical responses of the whole vertebra.     

Determining the Cellular Structure of Two Cereal Food Foams by X-ray Micro-Tomography

The cellular structure of two products, an extruded breakfast cereal and a short dough biscuit, was characterized by two different X-ray micro computed tomographic systems. Acquisitions were made by a compact desktop system Skyscan 1174 (Bruker μCT, Belgium) and at the European Synchrotron Radiation Facility (ESRF, beamline ID19, France) at different resolutions (voxel size of 6.5 μm, 7.5 μm, 16.2 μm and 25.8 μm). 3D images were processed for the density, the connectivity index and the granulometry of cells and cell walls. These experiments underlined the importance of the resolution for determination of quantitative measurements such as densities and thicknesses. The median width calculated for the cell walls distribution in the biscuit dropped from 141 to 50 μm when the voxel size changed from 25.8 to 7.5 μm. Images well showed that even though the food products had close values of porosity 0.6 and 0.7 for biscuit and extruded breakfast cereal respectively, their cellular structures were very different. The biscuit had small cells (median value of the distribution varied from 125 to 152 μm, according to resolution) and larger cell walls (50–141 μm) than the extrudate (32–109 μm) which, on the contrary, exhibited very large cells (307–400 μm). Beyond methodological issues, these differences could be clearly attributed to the differences of compositions and processes.     

Micro-finite element simulation of trabecular-bone post-yield behaviour – effects of material model,element size and type

Micro-finite element (micro-FE) analysis became a standard tool for the evaluation of trabecular bone mechanical properties. The accuracy of micro-FE models for linear analyses is well established. However, the accuracy of recently developed nonlinear micro-FE models for simulations of trabecular bone failure is not known. In this study, a trabecular bone specimen was compressed beyond the apparent yield point. The experiment was simulated using different micro-FE meshes with different element sizes and types, and material models based on cortical bone. The results from the simulations were compared with experimental results to study the effects of the different element and material models. It was found that a decrease in element size from 80 to 40 μm had little effect on predicted post-yield behaviour. Element type and material model had significant effects. Nevertheless, none of the established material models for cortical bone were able to predict the typical descent in the load-displacement curve seen during compression of trabecular bone.     

Mineral heterogeneity has a minor influence on the apparent elastic properties of human cancellous bone: a SRμCT-based finite element study

At the tissue level, the local material properties of human cancellous bone are heterogeneous due to constant remodelling. Since standard high-resolution computed tomography scanning methods are unable to capture this heterogeneity in detail, local differences in mineralisation are normally not incorporated in computational models. To investigate the effects of heterogeneous mineral distribution on the apparent elastic properties, 40 cancellous bone samples from the human femoral neck were scanned by means of synchrotron radiation microcomputed tomography (SRμCT). SRμCT-based micromechanical finite element models that accounted for mineral heterogeneity were compared with homogeneous models. Evaluation of the apparent stiffness tensor of both model types revealed that homogeneous models led to a minor but significant (p < 0.05) overestimation of the elastic properties of heterogeneous models by 2.18 ± 1.89%. Variation of modelling parameters did not affect the overestimation to a great extent. It was concluded that the heterogeneous mineralisation has only a minor influence on the apparent elastic properties of human cancellous bone.