Trabecular bone strains around a dental implant and associated micromotions—A micro-CT-based three-dimensional finite element study

The first objective of this computational study was to assess the strain magnitude and distribution within the three-dimensional (3D) trabecular bone structure around an osseointegrated dental implant loaded axially. The second objective was to investigate the relative micromotions between the implant and the surrounding bone. The work hypothesis adopted was that these virtual measurements would be a useful indicator of bone adaptation (resorption, homeostasis, formation).In order to reach these objectives, a μCT-based finite element model of an oral implant implanted into a Berkshire pig mandible was developed along with a robust software methodology. The finite element mesh of the 3D trabecular bone architecture was generated from the segmentation of μCT scans. The implant was meshed independently from its CAD file obtained from the manufacturer. The meshes of the implant and the bone sample were registered together in an integrated software environment. A series of non-linear contact finite element (FE) analyses considering an axial load applied to the top of the implant in combination with three sets of mechanical properties for the trabecular bone tissue was devised. Complex strain distribution patterns are reported and discussed. It was found that considering the Young’s modulus of the trabecular bone tissue to be 5, 10 and 15 GPa resulted in maximum peri-implant bone microstrains of about 3000, 2100 and 1400. These results indicate that, for the three sets of mechanical properties considered, the magnitude of maximum strain lies within an homeostatic range known to be sufficient to maintain/form bone. The corresponding micro-motions of the implant with respect to the bone microstructure were shown to be sufficiently low to prevent fibrous tissue formation and to favour long-term osseointegration.     

Elastic constants of composites formed from PMMA bone cement and anisotropic bovine tibial cancellous bone

An ultrasonic pulse-transit time technique is used to determine the nine orthotropic engineering constants of 32 cement-cancellous bone composites as a function of volume fractions of bone ranging from 0.0 to 0.4. The composites are manufactured using well-aligned bovine cancellous bone from the proximal end of the tibia and low viscosity bone cement. Selected composites are also subjected to mechanical compression tests to compare with the ultrasonic results. There is excellent correlation between the dynamic or ultrasonically determined moduli and the static or mechanically determined moduli; the dynamic moduli are approximately twice the static moduli and this difference is thought to be due to the effect of strain rate. An orthotropic model is assumed requiring nine independent elastic constants to be determined. The dynamic Young's modulus in the direction of major trabecular alignment, E1, increases linearly from 4.9 to 10.4 GPa as bone volume fraction increases from 0 to 0.4; dynamic E2 and E3 values increase from 4.9 to 7 GPa as bone volume fractions increase from 0 to 0.4, with E2 being slightly higher than E3. The dynamic shear modulus, G12, increases from 1.8 to 3.0 GPa, and G31 and G23 increase slightly from 1.8 to 2.2 GPa as bone volume fractions increase from 0 to 0.4. The Poisson's ratios are more sensitive than the Young's moduli and shear moduli to experimental error in the velocity measurements. The mechanically tested modulus (static modulus) in the direction of major trabecular alignment, E1, increases with volume fraction of bone from 2.4 to 4.4 GPa as the bone volume fraction increases from 0 to 0.25; static E2 and E3 values are either equal to or lower than that of pure PMMA.     

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.     

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.     

The effect of specimen geometry on the mechanical behaviour of trabecular bone specimens.

The effect of specimen geometry on the mechanical behaviour of trabecular bone specimens was studied by non-destructive uniaxial compression to 0.4% strain using cylindrical specimens with different sizes and length-to-diameter ratios, and by comparing cubic and cylindrical specimens with the same cross-sectional area. Both the length and the cross-sectional area of the specimen had a highly significant influence on the mechanical behaviour (p less than 0.0001). Within the actual range of length (2.75-11.0 mm) the normalized stiffness (Young's modulus) was related nearly linearly to the specimen length. This dependency on specimen length is suggested to be caused mainly by structural disintegrity of the trabecular specimens near the surface. The normalized stiffness (Young's modulus) was also positively correlated to the cross-sectional area. This dependency on cross-sectional area is probably due to friction-induced stress inhomogeneity at the platen-specimen interface. A cube with side length 6.5 mm or a cylindrical specimen with 7.5 mm diameter and 6.5 mm length are suggested as standard specimens for comparative studies on trabecular bone mechanics.     

Stiffness of compact bone: effects of porosity and density

Stiffness of compact bone is found to be highly and nonlinearly dependent on its porosity, its complement, bone volume fraction and apparent density. Elastic modulus decreases as a power (0.55) of increasing porosity and increases both as a power of increasing bone tissue volume (10.92) and increasing apparent density (7.4). These data indicate that small changes in the amount or density of compact bone tissue exert a more pronounced influence on its stiffness than would similar changes in trabecular bone.     

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.     

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.     

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.     

Fast and accurate specimen-specific simulation of trabecular bone elastic modulus using novel beam-shell finite element models

Elastic modulus and strength of trabecular bone are negatively affected by osteoporosis and other metabolic bone diseases. Micro-computed tomography-based beam models have been presented as a fast and accurate way to determine bone competence. However, these models are not accurate for trabecular bone specimens with a high number of plate-like trabeculae. Therefore, the aim of this study was to improve this promising methodology by representing plate-like trabeculae in a way that better reflects their mechanical behavior. Using an optimized skeletonization and meshing algorithm, voxel-based models of trabecular bone samples were simplified into a complex structure of rods and plates. Rod-like and plate-like trabeculae were modeled as beam and shell elements, respectively, using local histomorphometric characteristics. To validate our model, apparent elastic modulus was determined from simulated uniaxial elastic compression of 257 cubic samples of trabecular bone (4mm×4mm×4mm; 30μm voxel size; BIOMED I project) in three orthogonal directions using the beam-shell models and using large-scale voxel models that served as the gold standard. Excellent agreement (R(2)=0.97) was found between the two, with an average CPU-time reduction factor of 49 for the beam-shell models. In contrast to earlier skeleton-based beam models, the novel beam-shell models predicted elastic modulus values equally well for structures from different skeletal sites. It allows performing detailed parametric analyses that cover the entire spectrum of trabecular bone microstructures.     

European Society of Biomechanics S.M. Perren Award 2008: using temporal trends of 3D bone micro-architecture to predict bone quality

In longitudinal studies, three-dimensional (3D) bone images are acquired at sequential time points essentially resulting in four-dimensional (4D) data for an individual. Based on the 4D data, we propose to calculate temporal trends and project these trends to estimate future bone architecture. Multiple consecutive deformation fields, calculated with Demons deformable image registration algorithm, were extrapolated on a voxel-by-voxel basis. Test data were from in vivo micro-computed tomography (microCT) scans of the proximal tibia of Wistar rats that were either ovariectomized (OVX; N=5) or sham operated (SHAM; N=6). Measurements performed at baseline, 4 and 8 weeks were the basis to predict the 12 week data. Predicted and actual 12 week data were compared using qualitative (3D rendering) and quantitative (geometry, morphology and micro-finite element, microFE) methods. The results indicated a voxel-based linear extrapolation scheme yielded mean geometric errors that were smaller than the voxel size of 15 microm. Key morphological parameters that were estimated included bone volume ratio (BV/TV; mean error 0.4%, maximum error 9%), trabecular thickness (Tb.Th; -1.1%, 11%), connectivity density (Conn.D; 9.0%, 18.5%) and the apparent Young's modulus (E(1); 6.0%, 32%). These data demonstrated a promising and novel approach for quantitatively capturing in vivo bone dynamics at the local trabecular level. The method does not require an a priori understanding of the diseases state, and can provide information about the trends of the bone remodeling process that may be used for better monitoring and treatment of diseases such as osteoporosis.     

High-resolution finite element models with tissue strength asymmetry accurately predict failure of trabecular bone

The ability to predict trabecular failure using microstructure-based computational models would greatly facilitate study of trabecular structure–function relations, multiaxial strength, and tissue remodeling. We hypothesized that high-resolution finite element models of trabecular bone that include cortical-like strength asymmetry at the tissue level, could predict apparent level failure of trabecular bone for multiple loading modes. A bilinear constitutive model with asymmetric tissue yield strains in tension and compression was applied to simulate failure in high-resolution finite element models of seven bovine tibial specimens. Tissue modulus was reduced by 95% when tissue principal strains exceeded the tissue yield strains. Linear models were first calibrated for effective tissue modulus against specimen-specific experimental measures of apparent modulus, producing effective tissue moduli of (mean±S.D.) 18.7±3.4 GPa. Next, a parameter study was performed on a single specimen to estimate the tissue level tensile and compressive yield strains. These values, 0.60% strain in tension and 1.01% strain in compression, were then used in non-linear analyses of all seven specimens to predict failure for apparent tensile, compressive, and shear loading. When compared to apparent yield properties previously measured for the same type of bone, the model predictions of both the stresses and strains at failure were not statistically different for any loading case (p>0.15). Use of symmetric tissue strengths could not match the experimental data. These findings establish that, once effective tissue modulus is calibrated and uniform but asymmetric tissue failure strains are used, the resulting models can capture the apparent strength behavior to an outstanding level of accuracy. As such, these computational models have reached a level of fidelity that qualifies them as surrogates for destructive mechanical testing of real specimens.     

Convergence behavior of high-resolution finite element models of trabecular bone

The convergence behavior of finite element models depends on the size of elements used, the element polynomial order, and on the complexity of the applied loads. For high-resolution models of trabecular bone, changes in architecture and density may also be important. The goal of this study was to investigate the influence of these factors on the convergence behavior of high-resolution models of trabecular bone. Two human vertebral and two bovine tibial trabecular bone specimens were modeled at four resolutions ranging from 20 to 80 microns and subjected to both compressive and shear loading. Results indicated that convergence behavior depended on both loading mode (axial versus shear) and volume fraction of the specimen. Compared to the 20 microns resolution, the differences in apparent Young's modulus at 40 microns resolution were less than 5 percent for all specimens, and for apparent shear modulus were less than 7 percent. By contrast, differences at 80 microns resolution in apparent modulus were up to 41 percent, depending on the specimen tested and loading mode. Overall, differences in apparent properties were always less than 10 percent when the ratio of mean trabecular thickness to element size was greater than four. Use of higher order elements did not improve the results. Tissue level parameters such as maximum principal strain did not converge. Tissue level strains converged when considered relative to a threshold value, but only if the strains were evaluated at Gauss points rather than element centroids. These findings indicate that good convergence can be obtained with this modeling technique, although element size should be chosen based on factors such as loading mode, mean trabecular thickness, and the particular output parameter of interest.     

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.     

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.     

Synchrotron radiation micro-CT at the micrometer scale for the analysis of the three-dimensional morphology of microcracks in human trabecular bone

Bone quality is an important concept to explain bone fragility in addition to bone mass. Among bone quality factors, microdamage which appears in daily life is thought to have a marked impact on bone strength and plays a major role in the repair process. The starting point for all studies designed to further our understanding of how bone microdamage initiate or dissipate energy, or to investigate the impact of age, gender or disease, remains reliable observation and measurement of microdamage. In this study, 3D Synchrotron Radiation (SR) micro-CT at the micrometric scale was coupled to image analysis for the three-dimensional characterization of bone microdamage in human trabecular bone specimens taken from femoral heads. Specimens were imaged by 3D SR micro-CT with a voxel size of 1.4 μm. A new tailored 3D image analysis technique was developed to segment and quantify microcracks. Microcracks from human trabecular bone were observed in different tomographic sections as well as from 3D renderings. New 3D quantitative measurements on the microcrack density and morphology are reported on five specimens. The 3D microcrack density was found between 3.1 and 9.4/mm3 corresponding to a 2D density between 0.55 and 0.76 /mm2. The microcrack length and width measured in 3D on five selected microcrack ranged respectively from 164 μm to 209 μm and 100 μm to 120 μm. This is the first time that various microcracks in unloaded human trabecular bone--from the simplest linear crack to more complex cross-hatch cracks--have been examined and quantified by 3D imaging at this scale. The suspected complex morphology of microcracks is here considerably more evident than in the 2D observations. In conclusion, this technique opens new perspective for the 3D investigation of microcracks and the impact of age, disease or treatment.     

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.     

Static and dynamic finite element analyses of an idealized structural model of vertebral trabecular bone.

An idealized three-dimensional finite element model of a rodlike trabecular bone structure was developed to study its static and dynamic responses under compressive loading, considering the effects of bone marrow and apparent density. Static analysis of the model predicted hydraulic stiffening of trabecular bone due to the presence of bone marrow. The predicted power equation relating trabecular bone apparent elastic modulus to its apparent density was in good agreement with those of the reported experimental investigations. The ratio of the maximum stress in the trabecular bone tissue to its apparent stress had a high value, decreasing with increasing bone apparent density. Frequency analyses of the model predicted higher natural frequencies for the bone without marrow than those for the bone with marrow. Adding a mass relatively large compared to that of bone rendered a single-degree-of-freedom response. In this case, the resonant frequency was higher for the bone with marrow than that for the bone without marrow. The predicted vibrational measurement of apparent modulus was in good agreement with that of the static measurement, suggesting vibrational testing as a method for nondestructive measurement of trabecular bone elastic moduli.     

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.