Constitutive relationships of fabric, density, and elastic properties in cancellous bone architecture

High-resolution local trabecular strain within trabecular structure under cyclic loading

Trabecular bone structure is a complex microstructure consisting of rods and plates, which poses challenges for its mechanical characterization. Digital image correlation (DIC) offers the possibility to characterize the strain response on the surface of trabecular bone. This study employed DIC equipped with a telecentric lens to investigate the strain state of individual trabeculae within their trabecular structure by assessing the longitudinal strain of the trabeculae at both the middle and near the edges of the trabeculae. Due to the high-resolution of the used DIC system, local surface strain of trabeculae was analyzed too. Lastly, the correlation between longitudinal trabecular strain and the orientation and slenderness of the trabeculae was investigated. The results showed that the strain magnification close to the edge of the trabeculae was higher and reached up to 8-folds the strain along the middle of the trabeculae. On the contrary, no strain magnification was found for most of the trabeculae between the longitudinal trabecular strain along the middle of the trabeculae and the globally applied strain. High-resolution full-field strain maps were obtained on the surface of trabeculae showing heterogeneous strain distribution with increasing load. No significant correlation was found between longitudinal trabecular strain and its orientation or slenderness. These findings and the applied methodology can be used to broaden our understanding of the deformation mechanisms of trabeculae within the trabecular network.

High fat diet causes inferior vertebral structure and function without disc degeneration in RAGE-KO mice

Back pain and spinal pathologies are associated with obesity in juveniles and adults, yet studies identifying causal relationships are lacking and none investigate sex differences. This study determined if high fat (HF) diet causes structural and functional changes to vertebrae and intervertebral discs (IVDs); if these changes are modulated in mice with systematic ablation for the receptor for advanced glycation endproducts (RAGE-KO); and if these changes are sex-dependent. Wild-type (WT) and RAGE-KO mice were fed a low fat (LF) or HF diet for 12 weeks starting at 6 weeks, representing the juvenile population. HF diet led to weight/fat gain, glucose intolerance, and increased cytokine levels (IL-5, MIG, and RANTES); with less fat gain in RAGE-KO females. Most importantly, HF diet reduced vertebral trabecular bone volume fraction and compressive and shear moduli, without a modifying effect of RAGE-KO, but with a more pronounced effect in females. HF diet caused reduced cortical area fraction only in WT males. Neither HF diet nor RAGE-KO affected IVD degeneration grade. Biomechanical properties of coccygeal motion segments were affected by RAGE-KO but not diet, with some interactions identified. In conclusion, HF diet resulted in inferior vertebral structure and function with some sex differences, no IVD degeneration, and few modifying effects of RAGE-KO. These structural and functional deficiencies with HF diet provide further evidence that diet can affect spinal structures and may increase the risk for spinal injury and degeneration with aging and additional stressors. Back pain and spinal pathologies are associated with obesity in juveniles and adults, yet studies identifying causal relationships are lacking and none investigate sex differences.

Determination of anisotropic elastic parameters from morphological parameters of cancellous bone for osteoporotic lumbar spine

In biomechanics, large finite element models with macroscopic representation of several bones or joints are necessary to analyze implant failure mechanisms. In order to handle large simulation models of human bone, it is crucial to homogenize the trabecular structure regarding the mechanical behavior without losing information about the realistic material properties. Accordingly, morphology and fabric measurements of 60 vertebral cancellous bone samples from three osteoporotic lumbar spines were performed on the basis of X-ray microtomography (μCT) images to determine anisotropic elastic parameters as a function of bone density in the area of pedicle screw anchorage. The fabric tensor was mapped in cubic bone volumes by a 3D mean-intercept-length method. Fabric measurements resulted in a high degree of anisotropy (DA = 0.554). For the Young’s and shear moduli as a function of bone volume fraction (BV/TV, bone volume/total volume), an individually fit function was determined and high correlations were found (97.3 ≤ R² ≤ 99.1,p < 0.005). The results suggest that the mathematical formulation for the relationship between anisotropic elastic constants and BV/TV is applicable to current μCT data of cancellous bone in the osteoporotic lumbar spine. In combination with the obtained results and findings, the developed routine allows determination of elastic constants of osteoporotic lumbar spine. Based on this, the elastic constants determined using homogenization theory can enable efficient investigation of human bone using finite element analysis (FEA).

Graphical Abstract

Can DXA image-based deep learning model predict the anisotropic elastic behavior of trabecular bone?

3D image-based finite element (FE) and bone volume fraction (BV/TV)/fabric tensor modeling techniques are currently used to determine the apparent stiffness tensor of trabecular bone for assessing its anisotropic elastic behavior. Inspired by the recent success of deep learning (DL) techniques, we hypothesized that DL modeling techniques could be used to predict the apparent stiffness tensor of trabecular bone directly using dual-energy X-ray absorptiometry (DXA) images. To test the hypothesis, a convolutional neural network (CNN) model was trained and validated to predict the apparent stiffness tensor of trabecular bone cubes using their DXA images. Trabecular bone cubes obtained from human cadaver proximal femurs were used to obtain simulated DXA images as input, and the apparent stiffness tensor of the trabecular cubes determined by using micro-CT based FE simulations was used as output (ground truth) to train the DL model. The prediction accuracy of the DL model was evaluated by comparing it with the micro-CT based FE models, histomorphometric parameter based multiple linear regression models, and BV/TV/fabric tensor based multiple linear regression models. The results showed that DXA image-based DL model achieved high fidelity in predicting the apparent stiffness tensor of trabecular bone cubes (R² = 0.905-0.973), comparable to or better than the histomorphometric parameter based multiple linear regression and BV/TV/fabric tensor based multiple linear regression models, thus supporting the hypothesis of this study. The outcome of this study could be used to help develop DXA image-based DL techniques for clinical assessment of bone fracture risk.

Loss of longitudinal superiority marks the microarchitecture deterioration of osteoporotic cancellous bones

Osteoporosis (OP), a skeletal disease making bone mechanically deteriorate and easily fracture, is a global public health issue due to its high prevalence. It has been well recognized that besides bone loss, microarchitecture degradation plays a crucial role in the mechanical deterioration of OP bones, but the specific role of microarchitecture in OP has not been well clarified and quantified from mechanics perspective. Here, we successfully decoupled and identified the specific roles of microarchitecture, bone mass and tissue property in the failure properties of cancellous bones, through μCT-based digital modeling and finite element method simulations on bone samples from healthy and ovariectomy-induced osteoporotic mice. The results show that the microarchitecture of healthy bones exhibits longitudinal superiority in mechanical properties such as the effective stiffness, strength and toughness, which fits them well to bearing loads along their longitudinal direction. OP does not only reduce bone mass but also impair the microarchitecture topology. The former is mainly responsible for the mechanical degradation of bones in magnitude, wherever the latter accounts for the breakdown of their function-favorable anisotropy, the longitudinal superiority. Hence, we identified the microarchitecture-deterioration-induced directional mismatch between material and loading as a hazardous feature of OP and defined a longitudinal superiority index as measurement of the health status of bone microarchitecture. These findings provide useful insights and guidelines for OP diagnosis and treat assessment.

Trabecular Architecture and Mechanical Heterogeneity Effects on Vertebral Body Strength

PURPOSE OF REVIEW

We aimed to synthesize the recent work on the intra-vertebral heterogeneity in density, trabecular architecture and mechanical properties, its implications for fracture risk, its association with degeneration of the intervertebral discs, and its implications for implant design.

RECENT FINDINGS

As compared to the peripheral regions of the centrum, the central region of the vertebral body exhibits lower density and more sparse microstructure. As compared to the anterior region, the posterior region shows higher density. These variations are more pronounced in vertebrae from older persons and in those adjacent to degenerated discs. Mixed results have been reported in regard to variation along the superior-inferior axis and to relationships between the heterogeneity in density and vertebral strength and fracture risk. These discrepancies highlight that, first, despite the large amount of study of the intra-vertebral heterogeneity in microstructure, direct study of that in mechanical properties has lagged, and second, more measurements of vertebral loading are needed to understand how the heterogeneity affects distributions of stress and strain in the vertebra. These future areas of study are relevant not only to the question of spine fractures but also to the design and selection of implants for spine fusion and disc replacement. The intra-vertebral heterogeneity in microstructure and mechanical properties may be a product of mechanical adaptation as well as a key determinant of the ability of the vertebral body to withstand a given type of loading.

A method for fracture toughness measurement in trabecular bone using computed tomography, image correlation and finite element methods

The fracture resistance of load-bearing trabecular bone is adversely affected by diseases such as osteoporosis. However, there are few published measurements of trabecular bone fracture toughness due to the difficulty of conducting reliable tests in small specimens of this highly porous material. A new approach is demonstrated that uses digital volume correlation of X-ray computed tomographs to measure 3D displacement fields in which the crack shape and size can be objectively identified using a phase congruency analysis. The criteria for crack propagation, i.e. fracture toughness, can then be derived by finite element simulation, with knowledge of the elastic properties.

Empirical relationships between bone density and ultimate strength: A literature review

INTRODUCTION

Ultimate strength-density relationships for bone have been reported with widely varying results. Reliable bone strength predictions are crucial for many applications that aim to assess bone failure. Bone density and bone morphology have been proposed to explain most of the variance in measured bone strength. If this holds true, it could lead to the derivation of a single ultimate strength-density-morphology relationship for all anatomical sites.

METHODS

All relevant literature was reviewed. Ultimate strength-density relationships derived from mechanical testing of human bone tissue were included. The reported relationships were translated to ultimate strength-apparent density relationships and normalized with respect to strain rate. Results were grouped based on bone tissue type (cancellous or cortical), anatomical site, and loading mode (tension vs. compression). When possible, the relationships were compared to existing ultimate strength-density-morphology relationships.

RESULTS

Relationships that considered bone density and morphology covered the full spectrum of eight-fold inter-study difference in reported compressive ultimate strength-density relationships for trabecular bone. This was true for studies that tested specimens in different loading direction and tissue from different anatomical sites. Sparse data was found for ultimate strength-density relationships in tension and for cortical bone properties transverse to the main loading axis of the bone.

CONCLUSIONS

Ultimate strength-density-morphology relationships could explain measured strength across anatomical sites and loading directions. We recommend testing of bone specimens in other directions than along the main trabecular alignment and to include bone morphology in studies that investigate bone material properties. The lack of tensile strength data did not allow for drawing conclusions on ultimate strength-density-morphology relationships. Further studies are needed. Ideally, these studies would investigate both tensile and compressive strength-density relationships, including morphology, to close this gap and lead to more accurate evaluation of bone failure.

Multiscale design of artificial bones with biomimetic elastic microstructures

Cancellous bone is a highly porous, heterogeneous, and anisotropic material which can be found at the epiphyses of long bones and in the vertebral bodies. The hierarchical architecture makes cancellous bone a prime example of a lightweight natural material that combines strength with toughness. Better understanding the mechanics of cancellous bone is of interest for the diagnosis of bone diseases, the evaluation of the risk of fracture, and for the design of artificial bones and bone scaffolds for tissue engineering. A multiscale optimization method to maximize the stiffness of artificial bones using biomimetic cellular microstructures described by a finite set of geometrical micro-parameters is presented here. The most outstanding characteristics of its implementation are the use of: an interior point optimization algorithm, a precalculated response surface methodology for the evaluation of the elastic tensor of the microstructure as an analytical function of the micro-parameters, and the adjoint method for the computation of the sensitivity of the macroscopic mechanical response to the variation of the micro-parameters. The performance and effectiveness of the tool are evaluated by solving a problem that consists in finding the optimal distribution of the microstructures for a proximal end of a femur subjected to physiological loads. Two strategies for the specification of the solid volume fraction constraints are assessed. The results are compared with data of a computed tomography study of an actual human bone. The model successfully predicts the main features of the spatial arrangement of the trabecular and cortical microstructures of the natural bone.

Effect of fabric on the accuracy of computed tomography-based finite element analyses of the vertebra

Quantitative computed tomography (QCT)-based finite element (FE) models of the vertebra are widely used in studying spine biomechanics and mechanobiology, but their accuracy has not been fully established. Although the models typically assign material properties based only on local bone mineral density (BMD), the mechanical behavior of trabecular bone also depends on fabric. The goal of this study was to determine the effect of incorporating measurements of fabric on the accuracy of FE predictions of vertebral deformation. Accuracy was assessed by using displacement fields measured via digital volume correlation-applied to time-lapse microcomputed tomography (μCT)-as the gold standard. Two QCT-based FE models were generated from human L1 vertebrae (n = 11): the entire vertebral body and a cuboid-shaped portion of the trabecular centrum [dimensions: (20-30) × (15-20) × (15-20) mm3]. For axial compression boundary conditions, there was no difference (p = 0.40) in the accuracy of the FE-computed displacements for models using material properties based on local values of BMD versus those using material properties based on local values of fabric and volume fraction. However, when using BMD-based material properties, errors were higher for the vertebral-body models (8.4-50.1%) than cuboid models (1.5-19.6%), suggesting that these properties are inaccurate in the peripheral regions of the centrum. Errors also increased when assuming that the cuboid region experienced uniaxial loading during axial compression of the vertebra. These findings indicate that a BMD-based constitutive model is not sufficient for the peripheral region of the vertebral body when seeking accurate QCT-based FE modeling of the vertebra.

Cancellous bone and theropod dinosaur locomotion. Part I-an examination of cancellous bone architecture in the hindlimb bones of theropods

This paper is the first of a three-part series that investigates the architecture of cancellous ('spongy') bone in the main hindlimb bones of theropod dinosaurs, and uses cancellous bone architectural patterns to infer locomotor biomechanics in extinct non-avian species. Cancellous bone is widely known to be highly sensitive to its mechanical environment, and has previously been used to infer locomotor biomechanics in extinct tetrapod vertebrates, especially primates. Despite great promise, cancellous bone architecture has remained little utilized for investigating locomotion in many other extinct vertebrate groups, such as dinosaurs. Documentation and quantification of architectural patterns across a whole bone, and across multiple bones, can provide much information on cancellous bone architectural patterns and variation across species. Additionally, this also lends itself to analysis of the musculoskeletal biomechanical factors involved in a direct, mechanistic fashion. On this premise, computed tomographic and image analysis techniques were used to describe and analyse the three-dimensional architecture of cancellous bone in the main hindlimb bones of theropod dinosaurs for the first time. A comprehensive survey across many extant and extinct species is produced, identifying several patterns of similarity and contrast between groups. For instance, more stemward non-avian theropods (e.g. ceratosaurs and tyrannosaurids) exhibit cancellous bone architectures more comparable to that present in humans, whereas species more closely related to birds (e.g. paravians) exhibit architectural patterns bearing greater similarity to those of extant birds. Many of the observed patterns may be linked to particular aspects of locomotor biomechanics, such as the degree of hip or knee flexion during stance and gait. A further important observation is the abundance of markedly oblique trabeculae in the diaphyses of the femur and tibia of birds, which in large species produces spiralling patterns along the endosteal surface. Not only do these observations provide new insight into theropod anatomy and behaviour, they also provide the foundation for mechanistic testing of locomotor hypotheses via musculoskeletal biomechanical modelling.

Cancellous bone and theropod dinosaur locomotion. Part II-a new approach to inferring posture and locomotor biomechanics in extinct tetrapod vertebrates

This paper is the second of a three-part series that investigates the architecture of cancellous bone in the main hindlimb bones of theropod dinosaurs, and uses cancellous bone architectural patterns to infer locomotor biomechanics in extinct non-avian species. Cancellous bone is widely known to be highly sensitive to its mechanical environment, and therefore has the potential to provide insight into locomotor biomechanics in extinct tetrapod vertebrates such as dinosaurs. Here in Part II, a new biomechanical modelling approach is outlined, one which mechanistically links cancellous bone architectural patterns with three-dimensional musculoskeletal and finite element modelling of the hindlimb. In particular, the architecture of cancellous bone is used to derive a single 'characteristic posture' for a given species-one in which bone continuum-level principal stresses best align with cancellous bone fabric-and thereby clarify hindlimb locomotor biomechanics. The quasi-static approach was validated for an extant theropod, the chicken, and is shown to provide a good estimate of limb posture at around mid-stance. It also provides reasonable predictions of bone loading mechanics, especially for the proximal hindlimb, and also provides a broadly accurate assessment of muscle recruitment insofar as limb stabilization is concerned. In addition to being useful for better understanding locomotor biomechanics in extant species, the approach hence provides a new avenue by which to analyse, test and refine palaeobiomechanical hypotheses, not just for extinct theropods, but potentially many other extinct tetrapod groups as well.

Young's modulus of trabecular bone at the tissue level: A review

The tissue-level Young's modulus of trabecular bone is important for detailed mechanical analysis of bone and bone-implant mechanical interactions. However, the heterogeneity and small size of the trabecular struts complicate an accurate determination. Methods such as micro-mechanical testing of single trabeculae, ultrasonic testing, and nanoindentation have been used to estimate the trabecular Young's modulus. This review summarizes and classifies the trabecular Young's moduli reported in the literature. Information on species, anatomic site, and test condition of the samples has also been gathered. Advantages and disadvantages of the different methods together with recent developments are discussed, followed by some suggestions for potential improvement for future work. In summary, this review provides a thorough introduction to the approaches used for determining trabecular Young's modulus, highlights important considerations when applying these methods and summarizes the reported Young's modulus for follow-up studies on trabecular properties.

STATEMENT OF SIGNIFICANCE

The spongy trabecular bone provides mechanical support while maintaining a low weight. A correct measure of its mechanical properties at the tissue level, i.e. at a single-trabecula level, is crucial for analysis of interactions between bone and implants, necessary for understanding e.g. bone healing mechanisms. In this study, we comprehensively summarize the Young's moduli of trabecular bone estimated by currently available methods, and report their dependency on different factors. The critical review of different methods with recent updates is intended to inspire improvements in estimating trabecular Young's modulus. It is strongly suggested to report detailed information on the tested bone to enable statistical analysis in the future.

Trabecular architecture and joint loading of the proximal humerus in extant hominoids, Ateles, and Australopithecus africanus

OBJECTIVES

Several studies have investigated potential functional signals in the trabecular structure of the primate proximal humerus but with varied success. Here, we apply for the first time a "whole-epiphyses" approach to analysing trabecular bone in the humeral head with the aim of providing a more holistic interpretation of trabecular variation in relation to habitual locomotor or manipulative behaviors in several extant primates and Australopithecus africanus.

MATERIALS AND METHODS

We use a "whole-epiphysis" methodology in comparison to the traditional volume of interest (VOI) approach to investigate variation in trabecular structure and joint loading in the proximal humerus of extant hominoids, Ateles and A. africanus (StW 328).

RESULTS

There are important differences in the quantification of trabecular parameters using a "whole-epiphysis" versus a VOI-based approach. Variation in trabecular structure across knuckle-walking African apes, suspensory taxa, and modern humans was generally consistent with predictions of load magnitude and inferred joint posture during habitual behaviors. Higher relative trabecular bone volume and more isotropic trabeculae in StW 328 suggest A. africanus may have still used its forelimbs for arboreal locomotion.

DISCUSSION

A whole-epiphysis approach to analysing trabecular structure of the proximal humerus can help distinguish functional signals of joint loading across extant primates and can provide novel insight into habitual behaviors of fossil hominins.

Bone Mechanical Properties in Healthy and Diseased States

The mechanical properties of bone are fundamental to the ability of our skeletons to support movement and to provide protection to our vital organs. As such, deterioration in mechanical behavior with aging and/or diseases such as osteoporosis and diabetes can have profound consequences for individuals' quality of life. This article reviews current knowledge of the basic mechanical behavior of bone at length scales ranging from hundreds of nanometers to tens of centimeters. We present the basic tenets of bone mechanics and connect them to some of the arcs of research that have brought the field to recent advances. We also discuss cortical bone, trabecular bone, and whole bones, as well as multiple aspects of material behavior, including elasticity, yield, fracture, fatigue, and damage. We describe the roles of bone quantity (e.g., density, porosity) and bone quality (e.g., cross-linking, protein composition), along with several avenues of future research.

Dietary Advanced Glycation End Products Have Sex- and Age-Dependent Effects on Vertebral Bone Microstructure and Mechanical Function in Mice

Back pain is a leading cause of global disability that can arise from vertebral fracture and osteoporosis. Although poor general health and obesity are among the strongest risk factors for back pain, there is remarkably little known about how diet influences spinal diseases. Advanced glycation end-products (AGEs) are implicated in increased fracture risk, yet no studies investigated how dietary AGEs affect spinal health. We tested the hypothesis that high dietary AGE ingestion will diminish vertebral structure and function in a sex- and age-dependent manner. Female and male mice were fed low-AGE (L-AGE) or high-AGE (H-AGE) isocaloric diets for 6 and 18 months and multiple measurements of bone structure and function were taken. AGE levels in serum and cortical vertebrae were increased only for 6-month-old H-AGE female mice while blood glucose and body weight remained normal for all animals. When fed an H-AGE diet, 6-month-old female mice had inferior vertebral trabecular structure with decreased bone mineral density (BMD) and bone volume fraction. Biomechanical testing and analytical modeling further showed functional deterioration in 6-month-old H-AGE females with reduced shear and compression moduli, and maximum load to failure. At 18 months, H-AGE mice of both sexes had significant but small decreases in cortical BMD and thickness, yet functional biomechanical behaviors were not distinguishable from other aging changes. We conclude that an H-AGE diet, without diabetic or overweight conditions, diminished vertebral microstructure, mechanical behaviors, and fracture resistance in young female mice in a manner suggesting accelerated bone aging. © 2017 American Society for Bone and Mineral Research.

Calculation of cancellous bone elastic properties with the polarization-based FFT iterative scheme

The Fast Fourier Transform-based method, originally introduced by Moulinec and Suquet in 1994 has gained popularity for computing homogenized properties of composites. In this work, the method is used for the computational homogenization of the elastic properties of cancellous bone. To the authors' knowledge, this is the first study where the Fast Fourier Transform scheme is applied to bone mechanics. The performance of the method is analyzed for artificial and natural bone samples of 2 species: bovine femoral heads and implanted femurs of Hokkaido rats. Model geometries are constructed using data from X-ray tomographies, and the bone tissue elastic properties are measured using microindentation and nanoindentation tests. Computed results are in excellent agreement with those available in the literature. The study shows the suitability of the method to accurately estimate the fully anisotropic elastic response of cancellous bone. Guidelines are provided for the construction of the models and the setting of the algorithm.

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

INTRODUCTION

Previously, 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.

METHODS

Local 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%.

RESULTS

The 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%.

CONCLUSIONS

Of 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.

Mimetization of the elastic properties of cancellous bone via a parameterized cellular material

Bone tissue mechanical properties and trabecular microarchitecture are the main factors that determine the biomechanical properties of cancellous bone. Artificial cancellous microstructures, typically described by a reduced number of geometrical parameters, can be designed to obtain a mechanical behavior mimicking that of natural bone. In this work, we assess the ability of the parameterized microstructure introduced by Kowalczyk (Comput Methods Biomech Biomed Eng 9:135-147, 2006. doi: 10.1080/10255840600751473 ) to mimic the elastic response of cancellous bone. Artificial microstructures are compared with actual bone samples in terms of elasticity matrices and their symmetry classes. The capability of the parameterized microstructure to combine the dominant isotropic, hexagonal, tetragonal and orthorhombic symmetry classes in the proportions present in the cancellous bone is shown. Based on this finding, two optimization approaches are devised to find the geometrical parameters of the artificial microstructure that better mimics the elastic response of a target natural bone specimen: a Sequential Quadratic Programming algorithm that minimizes the norm of the difference between the elasticity matrices, and a Pattern Search algorithm that minimizes the difference between the symmetry class decompositions. The pattern search approach is found to produce the best results. The performance of the method is demonstrated via analyses for 146 bone samples.

Prediction of Mechanical Properties of the Cancellous Bone of the Mandibular Condyle

The mechanical properties of cancellous bone depend on the bone structure. The present study examined the extent to which the apparent stiffness of the cancellous bone of the human mandibular condyle can be predicted from its structure. Two models were compared. The first, a structure model, used structural parameters such as bone volume fraction and anisotropy to estimate the apparent stiffness. The second was a finite element model (FEM) of the cancellous bone. The bone structure was characterized by micro-computed tomography. The calculated stiffnesses of 24 bone samples were compared with measured stiffnesses. Both models could predict 89% of the variation in the measured stiffnesses. From the stiffness approximated by FEM in combination with the measured stiffness, the stiffness of the bone tissue was estimated to be 11.1 ± 3.2 GPa. It was concluded that both models could predict the stiffness of cancellous bone with adequate accuracy.

The effects of cracks on the quantification of the cancellous bone fabric tensor in fossil and archaeological specimens: a simulation study

Cancellous bone is very sensitive to its prevailing mechanical environment, and study of its architecture has previously aided interpretations of locomotor biomechanics in extinct animals or archaeological populations. However, quantification of architectural features may be compromised by poor preservation in fossil and archaeological specimens, such as post mortem cracking or fracturing. In this study, the effects of post mortem cracks on the quantification of cancellous bone fabric were investigated through the simulation of cracks in otherwise undamaged modern bone samples. The effect on both scalar (degree of fabric anisotropy, fabric elongation index) and vector (principal fabric directions) variables was assessed through comparing the results of architectural analyses of cracked vs. non-cracked samples. Error was found to decrease as the relative size of the crack decreased, and as the orientation of the crack approached the orientation of the primary fabric direction. However, even in the best-case scenario simulated, error remained substantial, with at least 18% of simulations showing a > 10% error when scalar variables were considered, and at least 6.7% of simulations showing a > 10° error when vector variables were considered. As a 10% (scalar) or 10° (vector) difference is probably too large for reliable interpretation of a fossil or archaeological specimen, these results suggest that cracks should be avoided if possible when analysing cancellous bone architecture in such specimens.

Disease and functional loading effect on the structural conformation and mechanical properties of the mandibular condyle in a transgenic rheumatoid arthritis murine model: an experimental study

AIM

The aim of the present study was to investigate the effect of rheumatoid arthritis (RA) and functional loading through diet modification on the structural conformation and the mechanical properties of the mandibular condyle in a transgenic mouse model and compare to healthy littermates.

MATERIALS AND METHODS

Four-week-old hybrid male mice from mixed background CBAxC57BL/6 were used. Four groups of animals were formed consisting of five animals each, either presenting RA (transgenic line hTNF 197), or wild-type (control), half receiving ordinary (hard) diet and half receiving soft diet within each category. Following sacrifice, resin-embedded and metallographically polished condylar specimens were evaluated employing scanning electron microscopy/ Energy dispersive x-ray spectroscopy and also tested for mechanical properties, through Vickers microhardness (HV100) measurements.

RESULTS

The multivariable analysis revealed significantly lower HV100 values for the RA groups after adjusting for diet (β = -10; 95% confidence interval: -16, -4; P = 0.001), while functional loading through diet modification did not appear as a significant predictor of the outcome.

CONCLUSIONS

There was evidence of compromised mechanical properties of the mandibular condylar bone for the diseased animals, whereas no association between functional loading and mechanical properties of the condyle could be established.

A review of trabecular bone functional adaptation: what have we learned from trabecular analyses in extant hominoids and what can we apply to fossils?

Many of the unresolved debates in palaeoanthropology regarding evolution of particular locomotor or manipulative behaviours are founded in differing opinions about the functional significance of the preserved external fossil morphology. However, the plasticity of internal bone morphology, and particularly trabecular bone, allowing it to respond to mechanical loading during life means that it can reveal greater insight into how a bone or joint was used during an individual's lifetime. Analyses of trabecular bone have been commonplace for several decades in a human clinical context. In contrast, the study of trabecular bone as a method for reconstructing joint position, joint loading and ultimately behaviour in extant and fossil non-human primates is comparatively new. Since the initial 2D studies in the late 1970s and 3D analyses in the 1990 s, the utility of trabecular bone to reconstruct behaviour in primates has grown to incorporate experimental studies, expanded taxonomic samples and skeletal elements, and improved methodologies. However, this work, in conjunction with research on humans and non-primate mammals, has also revealed the substantial complexity inherent in making functional inferences from variation in trabecular architecture. This review addresses the current understanding of trabecular bone functional adaptation, how it has been applied to hominoids, as well as other primates and, ultimately, how this can be used to better interpret fossil hominoid and hominin morphology. Because the fossil record constrains us to interpreting function largely from bony morphology alone, and typically from isolated bones, analyses of trabecular structure, ideally in conjunction with that of cortical structure and external morphology, can offer the best resource for reconstructing behaviour in the past.

Elastic Anisotropy of Trabecular Bone in the Elderly Human Vertebra

Knowledge of the nature of the elastic symmetry of trabecular bone is fundamental to the study of bone adaptation and failure. Previous studies have classified human vertebral trabecular bone as orthotropic or transversely isotropic but have typically obtained samples from only selected regions of the centrum. In this study, the elastic symmetry of human vertebral trabecular bone was characterized using microfinite element (μFE) analyses performed on 1019 cubic regions of side length equal to 5 mm, obtained via thorough sampling of the centrums of 18 human L1 vertebrae (age = 81.17 ± 7.7 yr; eight males and ten females). An optimization procedure was used to find the closest orthotropic representation of the resulting stiffness tensor for each cube. The orthotropic elastic constants and orientation of the principal elastic axes were then recorded for each cube and were compared to the constants predicted from Cowin's fabric-based constitutive model (Cowin, 1985, "The Relationship Between the Elasticity Tensor and the Fabric Tensor," Mech. Mater., 4(2), pp. 137-147.) and the orientation of the principal axes of the fabric tensor, respectively. Deviations from orthotropy were quantified by the "orthotropic error" (van Rietbergen et al., 1996, "Direct Mechanics Assessment of Elastic Symmetries and Properties of Trabecular Bone Architecture," J. Biomech., 29(12), pp. 1653-1657), and deviations from transverse isotropy were determined by statistical comparison of the secondary and tertiary elastic moduli. The orthotropic error was greater than 50% for nearly half of the cubes, and the secondary and tertiary moduli differed from one another (p < 0.0001). Both the orthotropic error and the difference between secondary and tertiary moduli decreased with increasing bone volume fraction (BV/TV; p ≤ 0.007). Considering only the cubes with an orthotropic error less than 50%, only moderate correlations were observed between the fabric-based and the μFE-computed elastic moduli (R2 ≥ 0.337; p < 0.0001). These results indicate that when using a criterion of 5 mm for a representative volume element (RVE), transverse isotropy or orthotropy cannot be assumed for elderly human vertebral trabecular bone. Particularly at low values of BV/TV, this criterion does not ensure applicability of theories of continuous media. In light of the very sparse and inhomogeneous microstructure found in the specimens analyzed in this study, further work is needed to establish guidelines for selecting a RVE within the aged vertebral centrum.

Clinical Use of Quantitative Computed Tomography-Based Finite Element Analysis of the Hip and Spine in the Management of Osteoporosis in Adults: the 2015 ISCD Official Positions-Part II

The International Society for Clinical Densitometry (ISCD) has developed new official positions for the clinical use of quantitative computed tomography (QCT)-based finite element analysis of the spine and hip. The ISCD task force for QCT reviewed the evidence for clinical applications and presented a report with recommendations at the 2015 ISCD Position Development Conference. Here we discuss the agreed upon ISCD official positions with supporting medical evidence, rationale, controversy, and suggestions for further study. Parts I and III address the clinical use of QCT of the hip, and the clinical feasibility of existing techniques for opportunistic screening of osteoporosis using CT scans obtained for other diagnosis such as colonography was addressed.

ULTRASONIC WAVE IS DETERMINED BY FABRIC TENSOR: AN APPLICATION TO CALCANEUS

The fabric tensor is a good measure for determining the mechanical properties of cancellous bone. Ultrasound is one method used to measure these mechanical properties. Ultrasound-generated speed of sound (SOS) measures the mechanical properties of cancellous bone. Thus, in this paper, we started with the fact that the fast wave in poroelastic theory is identical to the bulk wave velocity. We then formulate the equation for the fast wave in terms of fabric tensor for the calcaneus. The formulation in this paper is simpler than previously published results and will be easy to use in future experiments.

Experimentally-based multiscale model of the elastic moduli of bovine trabecular bone and its constituents

The elastic moduli of trabecular bone were modeled using an analytical multiscale approach. Trabecular bone was represented as a porous nanocomposite material with a hierarchical structure spanning from the collagen-mineral level to the trabecular architecture level. In parallel, compression testing was done on bovine femoral trabecular bone samples in two anatomical directions, parallel to the femoral neck axis and perpendicular to it, and the measured elastic moduli were compared with the corresponding theoretical results. To gain insights on the interaction of collagen and minerals at the nanoscale, bone samples were deproteinized or demineralized. After such processing, the treated samples remained as self-standing structures and were tested in compression. Micro-computed tomography was used to characterize the hierarchical structure of these three bone types and to quantify the amount of bone porosity. The obtained experimental data served as inputs to the multiscale model and guided us to represent bone as an interpenetrating composite material. Good agreement was found between the theory and experiments for the elastic moduli of the untreated, deproteinized, and demineralized trabecular bone.

Biomechanics and mechanobiology of trabecular bone: a review

Trabecular bone is a highly porous, heterogeneous, and anisotropic material which can be found at the epiphyses of long bones and in the vertebral bodies. Studying the mechanical properties of trabecular bone is important, since trabecular bone is the main load bearing bone in vertebral bodies and also transfers the load from joints to the compact bone of the cortex of long bones. This review article highlights the high dependency of the mechanical properties of trabecular bone on species, age, anatomic site, loading direction, and size of the sample under consideration. In recent years, high resolution micro finite element methods have been extensively used to specifically address the mechanical properties of the trabecular bone and provide unique tools to interpret and model the mechanical testing experiments. The aims of the current work are to first review the mechanobiology of trabecular bone and then present classical and new approaches for modeling and analyzing the trabecular bone microstructure and macrostructure and corresponding mechanical properties such as elastic properties and strength.

Gracility of the modern Homo sapiens skeleton is the result of decreased biomechanical loading

The postcranial skeleton of modern Homo sapiens is relatively gracile compared with other hominoids and earlier hominins. This gracility predisposes contemporary humans to osteoporosis and increased fracture risk. Explanations for this gracility include reduced levels of physical activity, the dissipation of load through enlarged joint surfaces, and selection for systemic physiological characteristics that differentiate modern humans from other primates. This study considered the skeletal remains of four behaviorally diverse recent human populations and a large sample of extant primates to assess variation in trabecular bone structure in the human hip joint. Proximal femur trabecular bone structure was quantified from microCT data for 229 individuals from 31 extant primate taxa and 59 individuals from four distinct archaeological human populations representing sedentary agriculturalists and mobile foragers. Analyses of mass-corrected trabecular bone variables reveal that the forager populations had significantly higher bone volume fraction, thicker trabeculae, and consequently lower relative bone surface area compared with the two agriculturalist groups. There were no significant differences between the agriculturalist and forager populations for trabecular spacing, number, or degree of anisotropy. These results reveal a correspondence between human behavior and bone structure in the proximal femur, indicating that more highly mobile human populations have trabecular bone structure similar to what would be expected for wild nonhuman primates of the same body mass. These results strongly emphasize the importance of physical activity and exercise for bone health and the attenuation of age-related bone loss.

Changes of elastic constants and anisotropy patterns in trabecular bone during disuse-induced bone loss assessed by poroelastic ultrasound

Currently, the approach most widely used to examine bone loss is the measurement of bone mineral density (BMD) using dual X-ray absorptiometry (DXA). However, bone loss due to immobilization creates changes in bone microarchitecture, which in turn are related to changes in bone mechanical function and competence to resist fracture.Unfortunately, the relationship between microarchitecture and mechanical function within the framework of immobilization and antiresorptive therapy has not being fully investigated. The goal of the present study was to investigate the structure–function relationship in trabecular bone in the real-world situations of a rapidly evolving osteoporosis(disuse), both with and without antiresorptive treatment. We evaluated the structure–function relationship in trabecular bone after bone loss (disuse-induced osteoporosis)and bisphosphonate treatment (antiresorptive therapy using risedronate) in canine trabecular bone using lCT and ultrasound wave propagation. Microstructure values determined from lCT images were used into the anisotropic poroelastic model of wave propagation in order to compute the apparent elastic constants (EC) and elastic anisotropy pattern of bone. Immobilization resulted in a significant reduction in trabecular thickness (Tb.Th) and bone volume fraction (BV/TV), while risedronate treatment combined with immobilization exhibited a lesser reduction in Tb.Th and BV/TV, suggesting that risedronate treatment decelerates bone loss, but it was unable to fully stop it. Risedronate treatment also increased the tissue mineral density (TMD), which when combined with the decrease in Tb.Th and BV/TV may explain the lack of significant differences invBMD in both immobilization and risedronate treated groups. Interestingly, changes inapparent EC were much stronger in the superior–inferior (SI) direction than in the medial–lateral (ML) and anterior–posterior (AP) anatomical directions, producing changes in elastic anisotropy patterns. When data were pooled together, vBMD was able to explain 58% of ultrasound measurements variability, a poroelastic wave propagation analytical model (i.e., BMD modulated by fabric directionality) was able to predict 81%of experimental wave velocity variability, and also explained 91% of apparent EC and changes in elastic anisotropy patterns. Overall, measurements of vBMD were unable to distinguish changes in apparent EC due to immobilization or risedronate treatment.However, anisotropic poroelastic ultrasound (PEUS) wave propagation was able to distinguish functional changes in apparent EC and elastic anisotropy patterns due to immobilization and antiresorptive therapy, providing an enhanced discrimination of anisotropic bone loss and the structure–function relationship in immobilized and risedronate-treated bone, beyond vBMD.

Effect of specimen-specific anisotropic material properties in quantitative computed tomography-based finite element analysis of the vertebra

Intra- and inter-specimen variations in trabecular anisotropy are often ignored in quantitative computed tomography (QCT)-based finite element (FE) models of the vertebra. The material properties are typically estimated solely from local variations in bone mineral density (BMD), and a fixed representation of elastic anisotropy ("generic anisotropy") is assumed. This study evaluated the effect of incorporating specimen-specific, trabecular anisotropy on QCT-based FE predictions of vertebral stiffness and deformation patterns. Orthotropic material properties estimated from microcomputed tomography data ("specimen-specific anisotropy"), were assigned to a large, columnar region of the L1 centrum (n = 12), and generic-anisotropic material properties were assigned to the remainder of the vertebral body. Results were compared to FE analyses in which generic-anisotropic properties were used throughout. FE analyses were also performed on only the columnar regions. For the columnar regions, the axial stiffnesses obtained from the two categories of material properties were uncorrelated with each other (p = 0.604), and the distributions of minimum principal strain were distinctly different (p ≤ 0.022). In contrast, for the whole vertebral bodies in both axial and flexural loading, the stiffnesses obtained using the two categories of material properties were highly correlated (R2 > 0.82, p < 0.001) with, and were no different (p > 0.359) from, each other. Only moderate variations in strain distributions were observed between the two categories of material properties. The contrasting results for the columns versus vertebrae indicate a large contribution of the peripheral regions of the vertebral body to the mechanical behavior of this bone. In companion analyses on the effect of the degree of anisotropy (DA), the axial stiffnesses of the trabecular column (p < 0.001) and vertebra (p = 0.007) increased with increasing DA. These findings demonstrate the need for accurate modeling of the peripheral regions of the vertebral body in analyses of the mechanical behavior of the vertebra.

μFEA successfully exhibits higher stresses and strains in microdamaged regions of whole vertebrae

Micro-finite element (μFE) modeling has shown promise in evaluating the structural integrity of trabecular bone. Histologic microcrack analyses have been compared to μFE models of trabecular bone cores to demonstrate the potential of this technique. To date this has not been achieved in whole bone structures, and comparisons of histologic microcrack and μFE results have been limited due to challenges in alignment of 2D sections with 3D data sets. The goal of this study was to ascertain if image registration can facilitate determination of a relationship between stresses and strains generated from μFE models of whole vertebrae and histologically identified microdamage. μFE models of three whole vertebrae, stained sequentially with calcein and fuchsin, were generated with accurate integration of element sets representing the histologic sections based on volumetric image registration. Displacement boundary conditions were applied to the μFE models based on registration of loaded and unloaded μCT images. Histologically labeled damaged regions were found to have significantly higher von Mises stresses and principle strains in the μFE models, as compared to undamaged regions. This work provides a new robust method for generating and histologically validating μFE models of whole bones that can represent trabecular damage resulting from complex physiologic loading.

A novel approach to estimate trabecular bone anisotropy using a database approach

Continuum finite element (FE) models of bones have become a standard pre-clinical tool to estimate bone strength. These models are usually based on clinical CT scans and material properties assigned are chosen as isotropic based only on the density distribution. It has been shown, however, that trabecular bone elastic behavior is best described as orthotropic. Unfortunately, the use of orthotropic models in FE analysis derived from CT scans is hampered by the fact that the measurement of a trabecular orientation (fabric) is not possible from clinical CT images due to the low resolution of such images. In this study, we explore the concept of using a database (DB) of high-resolution bone models to derive the fabric information that is missing in clinical images. The goal of this study was to investigate if models with fabric derived from a relatively small database can already produce more accurate results than isotropic models. A DB of 33 human proximal femurs was generated from micro-CT scans with a nominal isotropic resolution of 82 µm. Continuum FE models were generated from the images using a pre-defined mesh template in combination with an iso-anatomic mesh morphing tool. Each element within the mesh template is at a specific anatomical location. For each element within the cancellous bone, a spherical region around the element centroid with a radius of 2mm was defined. Bone volume fraction and the mean-intercept-length fabric tensor were analyzed for that region. Ten femurs were used as test cases. For each test femur, four different models were generated: (1) an orthotropic model based on micro-CT fabric measurements (gold standard), (2) an orthotropic model based on the fabric derived from the best-matched database model, (3) an isotropic-I model in which the fabric tensor was set to the identity tensor, and (4) a second isotropic-II model with its total bone stiffness fitted to the gold standard. An elastic-plastic damage model was used to simulate failure and post failure behavior during a fall to the side. The results show that all models produce a similar stress distribution. However, compared to the gold standard, both isotropic-I and II models underestimated the stress/damage distributions significantly. We found no significant difference between DB-derived and gold standard models. Compared to the gold standard, the isotropic-I models further underestimated whole bone stiffness by 26.3% and ultimate load by 14.5%, while these differences for the DB-derived orthotropic models were only 4.9% and 3.1% respectively. The results indicate that the concept of using a DB to estimate patient-specific anisotropic material properties can considerably improve the results. We expect that this approach can lead to more accurate results in particular for cases where bone anisotropy plays an important role, such as in osteoporotic patients and around implants.

Evaluating of bone healing around porous coated titanium implant and potential systematic bias on the traditional sampling method

INTRODUCTION

The mechanical properties of bone can largely be explained by bone density and the anisotropic orientation of the trabecular bone. The type of trabecular structure plays an important role in determining the mechanical properties of cancellous bone. Gap-healing and implant fixation could be affected by the various quality and quantity of bone in the local environment. Thus, implant fixation in one part might differ from the other part of the implant. This study aimed to investigate the influence of the sampling method on data evaluation.

MATERIAL AND METHODS

Titanium alloy implants (Biomet Inc.) of 10mm in length and 6mm in diameter were inserted bilaterally into the proximal humerus of 8 skeletally mature sheep. Thus two implants with a concentric gap of 2mm were implanted in each sheep. The gap was filled with allograft. Standardised surgical procedure was used. At sacrifice, 6 weeks after surgery, both proximal humeri were harvested. The specimens were randomized to superficial or profound groups. In the superficial group, mechanical testing or histological analysis was carried out on the superficial part of the implant. In the profound group, the mechanical testing or histological analysis was performed on the profound part of the implant.

RESULT

The mechanical fixation, bone volume and bone ongrowth showed no statistically significant differences. Mechanical test demonstrated a slight tendency to increased strength and failure energy were observed in the superficial group. Histomorphomety revealed bone ongrowth was slightly increased and volume fraction was decreased in the profound group.

CONCLUSION

Histological analysis and mechanical testing can be applied to the superficial or profound part of the implant.

Towards patient-specific material modeling of trabecular bone post-yield behavior

Bone diseases such as osteoporosis are one of the main causes of bone fracture and often result in hospitalization and long recovery periods. Researchers are aiming to develop new tools that consider the multiple determinants acting at the different scales of bone, and which can be used to clinically estimate patient-specific fracture risk and also assess the efficacy of new therapies. The main step towards this goal is a deep understanding of the bone organ, and is achieved by modeling the complexity of the structure and the high variability of the mechanical outcome. This review uses a hierarchical approach to evaluate bone mechanics at the macroscale, microscale, and nanoscale levels and the interactions between scales. The first section analyzes the experimental evidence of bone mechanics in the elastic and inelastic regions, microdamage generation, and post-yield toughening mechanisms from the organ level to the ultrastructural level. On the basis of these observations, the second section provides an overview of the constitutive models available to describe bone mechanics and predict patient-specific outcomes. Overall, the role of the hierarchical structure of bone and the interplay between each level is highlighted, and their effect is evaluated in terms of modeling biological variability and patient specificity.

A symmetry invariant formulation of the relationship between the elasticity tensor and the fabric tensor

The fabric tensor is employed as a quantitative stereological measure of the structural anisotropy in the pore architecture of a porous medium. Earlier work showed that the fabric tensor can be used additionally to the porosity to describe the anisotropy in the elastic constants of the porous medium. This contribution presents a reformulation of the relationship between fabric tensor and anisotropic elastic constants that is approximation free and symmetry-invariant. From specific data on the elastic constants and the fabric, the parameters in the reformulated relationship can be evaluated individually and efficiently using a simplified method that works independent of the material symmetry. The well-behavedness of the parameters and the accuracy of the method was analyzed using the Mori-Tanaka model for aligned ellipsoidal inclusions and using Buckminster Fuller's octet-truss lattice. Application of the method to a cancellous bone data set revealed that employing the fabric tensor allowed explaining 75-90% of the total variance. An implementation of the proposed methods was made publicly available.

Micro-architecture and mineralization of the human alveolar bone obtained with microCT

OBJECTIVES

The primary dental implant stability depends on the location of the implant in the jaw. This study analysed the architecture and mineralization of the trabecular bone at different jaw locations and thereby identified potential prognostic factors for implant failure. It has checked the hypotheses: (1) the mandible contains more compact and less mineralized trabecular bone than the maxilla and (2) within the mandible the trabecular bone is more compact and less mineralized in the anterior region.

METHODS

Alveolar bone specimens were produced from the cadavers of ten humans (7 males and 3 females; mean age: 73.7±12.5 years) and scanned with a high-resolution microCT system. Volumes of interest were chosen next to the roots of molars and incisors in both the maxilla and mandible. Several morphological parameters as well as the tissue mineral density were determined.

RESULTS

The alveolar bone specimens had a very high bone volume fraction (mean=0.31) with large differences (SD=0.17) between and within subjects. Yet several significant differences were found between the maxilla and the mandible. The bone volume fraction and trabecular thickness were significantly higher in the mandible than in the maxilla (p<0.01). But the tissue mineral density was not significantly different.

CONCLUSIONS

A higher primary implant stability coincides with a higher bone volume fraction and degree of anisotropy Although local differences remain more important for implantology, the results suggest that the micro-structure also affects the implant stability. The tissue mineral density seems to have no predictive value.

Predicting trabecular bone elastic properties from measures of bone volume fraction and fabric on the basis of micromagnetic resonance images

The relationship between fabric (a measure of structural anisotropy) and elastic properties of trabecular bone was examined by invoking morphology and homogenization theory on the basis of micromagnetic resonance images from the distal tibia in specimens (N = 30) and human subjects (N = 16) acquired at a 160 × 160 × 160 μm(3) voxel size. The fabric tensor was mapped in 7.5 × 7.5 × 7.5 mm(3) cubic subvolumes by a three-dimensional mean-intercept-length method. Elastic constants (three Young's and three shear moduli) were derived from linear microfinite element simulations of three-dimensional grayscale bone volume fraction-mapped images. In the specimen data, moduli fit power laws of bone volume fraction (bone volume/total volume) for all three test directions and subvolumes (R(2) = 0.92-0.98) with exponents ranging from 1.3 to 1.8. Weaker linear relationships were found for the in vivo data because of a narrower range in bone volume/total volume. When pooling the data for all test directions and subvolumes, bone volume/total volume predicted elastic moduli less well in the specimens (mean R(2) = 0.74) and not at all in vivo. A model of bone volume/total volume and fabric was highly predictive of microfinite element-derived Young's moduli: mean R(2) s of 0.98 and 0.82 (in vivo). The results show that fabric, an important predictor of bone mechanical properties, can be assessed in the limited resolution and signal-to-noise ratio regime of micromagnetic resonance images.