Biomechanical properties across trabeculae from the proximal femur of normal and ovariectomised sheep

The elastic behaviour of trabecular bone is a function not only of bone volume and architecture, but also of tissue material properties. Variation in tissue modulus can have a substantial effect on the biomechanical properties of trabecular bone. However, the nature of tissue property variation within a single trabecula is poorly understood. This study uses nanoindentation to determine the mechanical properties of bone tissue in individual trabeculae. Using an ovariectomised ovine model, the modulus and hardness distribution across trabeculae were measured. In both normal and ovariectomised bone, the modulus and hardness were found to increase towards the core of the trabeculae. Across the width of the trabeculae, the modulus was significantly less in the ovariectomised bone than in the control bone. However, in contrast to this hardness was found not to differ significantly between the two groups. This study provides valuable information on the variation of mechanical material properties in healthy and diseased trabecular bone tissue. The results of the current study will be useful in finite element modelling where more accurate values of trabecular bone modulus will enable the prediction of the macroscale behaviour 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.     

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

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

Trabecular surface remodeling simulation for cancellous bone using microstructural voxel finite element models

A computational simulation method for three-dimensional trabecular surface remodeling was proposed, using voxel finite element models of cancellous bone, and was applied to the experimental data. In the simulation, the trabecular microstructure was modeled based on digital images, and its morphological changes due to surface movement at the trabecular level were directly expressed by removing/adding the voxel elements from/to the trabecular surface. A remodeling simulation at the single trabecular level under uniaxial compressive loading demonstrated smooth morphological changes even though the trabeculae were modeled with discrete voxel elements. Moreover, the trabecular axis rotated toward the loading direction with increasing stiffness, simulating functional adaptation to the applied load. In the remodeling simulation at the trabecular structural level, a cancellous bone cube was modeled using a digital image obtained by microcomputed tomography (microCT), and was uniaxially compressed. As a result, the apparent stiffness against the applied load increased by remodeling, in which the trabeculae reoriented to the loading direction. In addition, changes in the structural indices of the trabecular architecture coincided qualitatively with previously published experimental observations. Through these studies, it was demonstrated that the newly proposed voxel simulation technique enables us to simulate the trabecular surface remodeling and to compare the results obtained using this technique with the in vivo experimental data in the investigation of the adaptive bone remodeling phenomenon.     

Quantification of trabecular bone microdamage using the virtual internal bond model and the individual trabeculae segmentation technique

Trabecular bone microdamage significantly influences the skeletal integrity and bone remodelling process. In this paper a novel constitutive model, called the virtual internal bond model (VIB), was adopted for simulating the damage behaviour of bone tissue. A unique 3D image analysis technique, named individual trabeculae segmentation, was used to analyse the effects of microarchitectures on the damage behaviours of trabecular bone. We demonstrated that the process of initiation and accumulation of microdamage in trabecular bone samples can be captured by the VIB-embedded finite-element method simulation without a separate fracture criterion. Our simulation results showed that the microdamage can occur at as early as about 0.2–0.4% apparent strain, and a large volume of microdamage was accumulated around the apparent yield strain. In addition we found that the plate-like trabeculae, especially the longitudinal ones, take crucial roles in the microdamage behaviours of trabecular bone.     

Effect of trabecular curvature on the stiffness of trabecular bone

Simplified structural models of trabecular bone have been used to model various forms of trabecular variability. The structural effects of variability of direction, length and thickness of the trabeculae have been studied using 'lattice-type' finite element models. However, many of the trabeculae are not perfectly straight, and have a small degree of curvature. The objective of this study is to quantify the influence of small curvatures of the trabeculae on the effective modulus of trabecular bone, in the principal material direction. An analytical analysis of the effect of curvature on a single trabecula is performed, utilizing the concept of cellular-solid models. Closed-form expressions are derived for the effect of curvature on the flexibility in the principal material direction. For comparison, expressions are derived for the flexibility of a straight oblique element, representing angular variability. A quantitative comparison is presented, which is dependent on the thickness of the trabeculae. It was found that small curvatures have a large effect on the stiffness of the trabecular structure. This effect is largest for thin trabeculae, and decreases for thick trabeculae. The stiffness of the trabecular structure can be reduced by a factor of up to four for thin trabeculae and up to two for thick trabeculae, even for small curvatures. The flexibility of curved elements is found to be larger than the flexibility of oblique elements with similar eccentricities. Thus it seems that curvature might play a role in determining the effective modulus of trabecular bone.     

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.     

Heterogeneous linear elastic trabecular bone modelling using micro-CT attenuation data and experimentally measured heterogeneous tissue properties

High-resolution voxel-based finite element software, such as FEEBE developed at the NCBES, is widely used for studying trabecular bone at the micro-scale. A new approach to determine heterogeneous bone tissue material properties for computational models was proposed in this study. The specimen-specific range of tissue moduli across strut width was determined from nanoindentation testing. This range was mapped directly using linear interpolation to that specimen's micro-computed tomography (microCT) grey value range as input material properties for finite element analysis. The method was applied to cuboid trabecular bone samples taken from eight, 4-year-old (skeletally mature) ovine L5 vertebrae. Before undergoing experimental uniaxial compression tests, the samples were microCT scanned and 30 microm resolution finite element models were generated. The linear elastic finite element models were compressed to 1% strain. This material property assignment method for computational models accurately reproduced the experimentally determined apparent modulus and concentrations of stress at locations of failure.     

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.     

The role of an effective isotropic tissue modulus in the elastic properties of cancellous bone.

Conceptually, the elastic characteristics of cancellous bone could be predicted directly from the trabecular morphology--or architecture--and by the elastic properties of the tissue itself. Although hardly any experimental evidence exists, it is often implicitly assumed that tissue anisotropy has a negligible effect on the apparent elastic properties of cancellous bone. The question addressed in this paper is whether this is actually true. If it is, then micromechanical finite element analysis (micro-FEA) models, representing trabecular architecture, using an 'effective isotropic tissue modulus' should be able to predict apparent elastic properties of cancellous bone. To test this, accurate multi-axial compressive mechanical tests of 29 whale bone specimens were simulated with specimen-specific micro-FEA computer models built from true three-dimensional reconstructions. By scaling the micro-FEA predictions by a constant tissue modulus, 92% of the variation of Young's moduli determined experimentally could be explained. The correlation even increased to 95% when the micro-FEA moduli were scaled to the isotropic tissue moduli of individual specimens. Excellent agreement was also found in the elastic symmetry axes and anisotropy ratios. The prediction of Poisson's ratios was somewhat less precise at 85% correlation. The results support the hypothesis; for practical purposes, the concept of an 'effective isotropic tissue modulus' concept is a viable one. They also suggest that the value of such a modulus for individual cases might be inferred from the average tissue density, hence the degree of mineralization. Future studies must clarify how specific the tissue modulus should be for different types of bone if adequate predictions of elastic behavior are to be made in this way.     

A New Computational Efficient Approach for Trabecular Bone Analysis using Beam Models Generated with Skeletonized Graph Technique

Micro-finite element (FE) analysis is a well established technique for the evaluation of the elastic properties of trabecular bone, but is limited in its application due to the large number of elements that it requires to represent the complex internal structure of the bone. In this paper, we present an alternative FE approach that makes use of a recently developed 3D-Line Skeleton Graph Analysis (LSGA) technique to represent the complex internal structure of trabecular bone as a network of simple straight beam elements in which the beams are assigned geometrical properties of the trabeculae that they represent. Since an enormous reduction of cputime can be obtained with this beam modeling approach, ranging from approximately 1,200 to 3,600 for the problems investigated here, we think that the FE modeling technique that we introduced could potentially constitute an interesting alternative for the evaluation of the elastic mechanical properties of trabecular bone.     

A new computational efficient approach for trabecular bone analysis using beam models generated with skeletonized graph technique

Micro-finite element (FE) analysis is a well established technique for the evaluation of the elastic properties of trabecular bone, but is limited in its application due to the large number of elements that it requires to represent the complex internal structure of the bone. In this paper, we present an alternative FE approach that makes use of a recently developed 3D-Line Skeleton Graph Analysis (LSGA) technique to represent the complex internal structure of trabecular bone as a network of simple straight beam elements in which the beams are assigned geometrical properties of the trabeculae that they represent. Since an enormous reduction of cputime can be obtained with this beam modeling approach, ranging from approximately 1,200 to 3,600 for the problems investigated here, we think that the FE modeling technique that we introduced could potentially constitute an interesting alternative for the evaluation of the elastic mechanical properties of trabecular bone.     

Fast tool for evaluation of iliac crest tissue elastic properties using the reduced-basis methods

Computationally expensive finite element (FE) methods are generally used for indirect evaluation of tissue mechanical properties of trabecular specimens, which is vital for fracture risk prediction in the elderly. This work presents the application of reduced-basis (RB) methods for rapid evaluation of simulation results. Three cylindrical transiliac crest specimens (diameter: 7.5 mm, length: 10-12 mm) were obtained from healthy subjects (20 year-old, 22 year-old, and 24 year-old females) and scanned using microcomputed tomography imaging. Cubic samples of dimensions 5×5×5 mm(3) were extracted from the core of the cylindrical specimens for FE analysis. Subsequently, a FE solution library (test space) was constructed for each of the specimens by varying the material property parameters: tissue elastic modulus and Poisson's ratio, to develop RB algorithms. The computational speed gain obtained by the RB methods and their accuracy relative to the FE analysis were evaluated. Speed gains greater than 4000 times, were obtained for all three specimens for a loss in accuracy of less than 1% in the maxima of von-Mises stress with respect to the FE-based value. The computational time decreased from more than 6 h to less than 18 s. RB algorithms can be successfully utilized for real-time reliable evaluation of trabecular bone elastic properties.     

Geometrical and material parameters to assess the macroscopic mechanical behaviour of fresh cranial bone samples

The present study aims at providing quantitative data for the personalisation of geometrical and mechanical characteristics of the adult cranial bone to be applied to head FE models. A set of 351 cranial bone samples, harvested from 21 human skulls, were submitted to three-point bending tests at 10 mm/min. For each of them, an apparent elastic modulus was calculated using the beam's theory and a density-dependant beam inertia. Thicknesses, apparent densities and percentage of ash weight were also measured. Distributions of characteristics among the different skull bones show their symmetry and their significant differences between skull areas. A data analysis was performed to analyse potential relationship between thicknesses, densities and the apparent elastic modulus. A specific regression was pointed out to estimate apparent elastic modulus from the product of thickness by apparent density. These results offer quantitative tools in view of personalising head FE models and thus improve definition of local injury criteria for this body part.     

Nanostructure and elastic modulus of single trabecula in bovine cancellous bone

We aimed to investigate the elastic modulus of trabeculae using tensile tests and assess the effects of nanostructure at the hydroxyapatite (HAp) crystal scale on the elastic modulus. In the experiments, 18 trabeculae that were at least 3 mm in length in the proximal epiphysis of three adult bovine femurs were used. Tensile tests were conducted using a small tensile testing device coupled with microscopy under air-dried condition. The c-axis orientation of HAp crystals and the degree of orientation were measured by X-ray diffraction. To observe the deformation behavior of HAp crystals under tensile loading, the same tensile tests were conducted in X-ray diffraction measurements. The mineral content of specimens was evaluated using energy dispersive X-ray spectrometry. The elastic modulus of a single trabecula varied from 4.5 to 23.6 GPa, and the average was 11.5±5.0 GPa. The c-axis of HAp crystals was aligned with the trabecular axis and the crystals were lineally deformed under tensile loading. The ratio of the HAp crystal strain to the tissue strain (strain ratio) had a significant correlation with the elastic modulus (r=0.79; P<0.001). However, the mineral content and the degree of orientation did not vary widely and did not correlate with the elastic modulus in this study. It suggests that the strain ratio may represent the nanostructure of a single trabecula and would determine the elastic modulus as well as mineral content and orientation.     

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

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

Simulation of vertebral trabecular bone loss using voxel finite element analysis

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

Micro-CT检测中等强度跑台运动对去卵巢大鼠腰椎微结构的影响

目的用micro-CT方法,评估中等强度跑台运动对去卵巢大鼠腰椎微结构的影响。方法将30只3月龄雌性SD大鼠按体重分层后随机分为假手术、去卵巢静止和去卵巢运动三个组。运动组每周进行4次45min、速度18 m/min、坡度5°的跑台训练。正式运动处理14周时,取第2腰椎检测骨密度,取第4腰椎行micro-CT分析及三维结构重建;取第3腰椎椎体进行椎体压缩实验。结果去卵巢运动组第2腰椎骨密度、第3腰椎最大载荷、最大应力和弹性模量以及第4腰椎骨小梁体积和骨小梁数目显著高于去卵巢静止组,骨小梁分离度显著低于去卵巢静止组,而骨小梁厚度无显著变化。结论中等强度跑台运动能改善去卵巢大鼠腰椎的微结构。     

A global relationship between trabecular bone morphology and homogenized elastic properties.

An alternative concept of the relationship between morphological and elastic properties of trabecular bone is presented and applied to human tissue from several anatomical locations using a digital approach. The three-dimensional morphology of trabecular bone was assessed with a microcomputed tomography system and the method of directed secants as well as the star volume procedure were used to compute mean intercept length (MIL) and average bone length (ABL) of 4 mm cubic specimens. Assuming isotropic elastic properties for the trabecular tissue, the general elastic tensors of the bone specimens were determined using the homogenization method and the closest orthotropic tensors were calculated with an optimization algorithm. The assumption of orthotropy for trabecular bone was found to improve with specimen size and hold within 6.1 percent for a 4 mm cube size. A strong global relationship (r2 = 0.95) was obtained between fabric and the orthotropic elastic tensor with a minimal set of five constants. Mean intercept length and average bone length provided an equivalent power of prediction. These results support the hypothesis that the elastic properties of human trabecular bone from an arbitrary anatomical location can be estimated from an approximation of the anisotropic morphology and a prior knowledge of tissue properties.     

Pullout strength of suture anchors: Effect of mechanical properties of trabecular bone

This study investigated the relationships between trabecular microstructure and elastic modulus, compressive strength, and suture anchor pullout strength. Twelve fresh-frozen humeri underwent mechanical testing followed by micro-computed tomography (μCT). Either compression testing of cylindrical bone samples or pullout testing using an Arthrex 5 mm Corkscrew was performed in synthetic sawbone or at specific locations in the humerus such as the greater tuberosity, lesser tuberosity, and humeral head. Synthetic sawbone underwent identical mechanical testing and μCT analysis. Bone volume fraction (BVF), structural model index (SMI), trabecular thickness (TbTh), trabecular spacing (TbSp), trabecular number (TbN), and connectivity density were compared against modulus, compressive strength, and pullout strength in both materials. In cadaveric bone, modulus showed correlations to all of the microstructural properties, while compressive and pullout strength were only correlated to BVF, SMI, and TbSp. The microstructure of synthetic bone differed from cadaveric bone as SMI and TbTh showed little variation across the densities tested. Therefore, SMI and TbTh were the only microstructural properties that did not show correlations to the mechanical properties tested in synthetic bone. This study helps identify key microstructure–property relationships in cadaveric and synthetic bone as well as illustrate the similarities and differences between cadaveric and synthetic bone as biomechanical test materials.