Nanoindentation of osteonal bone lamellae

Microcrack behavior in bone: Stress field analysis at osteon cement line tips

Bone microstructure governs microcrack propagation complexity. Current research, relying on linear elastic fracture mechanics, inadequately considers authentic multi-level structures, like cement lines and osteons, impacting stress intensity at cracks. This study, by constructing models encompassing single or multiple osteons, delves into the influence of factors like crack length, osteon radius, and modulus ratio on the stress intensity factor at the crack tip. Employing a fracture mechanics phase-field approach to simulate crack propagation paths, it particularly explores the role of cement lines as weak interfaces in crack extension. The aim is to comprehensively and systematically elucidate the critical factors of bone microstructure in the context of crack propagation.

Comprehensively characterizing heterogeneous and transversely isotropic properties of femur cortical bones

Comprehensive characterization of the transversely isotropic mechanical properties of long bones along both the longitudinal and circumferential gradients is crucial for developing accurate mathematical models and studying bone biomechanics. In addition, mechanical testing to derive elastic, plastic, and failure properties of bones is essential for modeling plastic deformation and failure of bones. To achieve these, we machined a total of 336 cortical specimens, including 168 transverse and 168 longitudinal specimens, from four different quadrants of seven different sections of 3 bovine femurs. We conducted three-point bending tests of these specimens at a loading rate of 0.02 mm/s. Young's modulus, yield stress, tangential modulus, and effective plastic strain for each specimen were derived from correction equations based on classical beam theory. Our statistical analysis reveals that the longitudinal gradient has a significant effect on the Young's modulus, yield stress, and tangential modulus of both longitudinal and transverse specimens, whereas the circumferential gradient significantly influences the Young's modulus, yield stress, and tangential modulus of transverse specimens only. The differences in Young's modulus and yield stress between longitudinal specimens from different sections are greater than 40%, whereas those between transverse specimens are approximately 30%. The Young's modulus and yield stress of transverse specimens in the anterior quadrant were 18.81%/15.46% and 18.34%/14.88% higher than those in the posterior and lateral quadrants, respectively. There is no significant interaction between the longitudinal gradient and the circumferential gradient. Considering the transverse isotropy, it is crucial to consider loading direction when investigating the impact of circumferential gradients in the anterior, lateral, medial, and posterior directions. Our findings indicate that the conventional assumption of homogeneity in simulating the cortical bone of long bones may have limitations, and researchers should consider the anatomical position and loading direction of femur specimens for precise prediction of mechanical responses.

Elasticity and material anisotropy of lamellar cortical bone in adult bovine tibia characterized via AFM nanoindentation

The research focuses on the evaluation of the mechanical properties of osteonal cortical bone at the lamellar level. Elastic properties of the mid-diaphysis region of the bovine tibia are investigated via cantilever-based nanoindentation at the submicron length scale utilizing Atomic Force Microscopy, where the force-displacement curves are used for the elastic assessment using the Derjaguin-Muller-Toropov model to calculate indentation modulus. Variations of the modulus and the directional mechanical response of the osteonal bone at different distances from the Haversian canal are investigated. Additionally, the effects of demineralization on the indentation modulus are discussed. It was found that in the axial direction, the first and last untreated thick lamella layers show a significant indentation modulus difference compared to all other layers (4.26 ± 0.4 and 4.6 ± 0.3 GPa vs ∼3.5 GPa). On the other hand, the indentation modulus of transverse thick lamella layers shows a periodic variation between ∼3 ± 0.7 GPa and ∼4 ± 0.3 GPa from near the Haversian canal to near the interstitial bone. A periodic variation in the anisotropy ratio was found. Mineral content was quantified via energy-dispersive X-ray microanalysis at different levels of mineralization and shows a positive correlation with the indentation modulus.

The influence of different shaped osteocyte lacunae on microcrack initiation and propagation

BACKGROUND

The morphology of osteocyte lacunae varies in bones of different ages and bone pathologies. Osteocyte lacunae can cause stress concentration and initiate microcracks. However, the influence of changes in osteocyte lacunar shape on microcrack is unknown. Therefore, the aim of this study was to determine the effects of osteocyte lacunae with different shapes on microcrack initiation and propagation.

METHODS

Osteon models containing osteocyte lacunae with different shapes were established. The progressive damage analysis method, based on computer simulations, was used to study the evolution of microdamage within the osteon, including the processes of intralaminar and interlaminar microdamage.

FINDINGS

Models with larger DoE values can effectively delay or prevent the formation of linear microcracks, which ensures high fracture toughness of cortical bone. It is subjected to stronger mechanical stimulation, making it more sensitive to loads. Models with smaller DoE values increase the load threshold for microdamage generation and reduces its impact on bone mechanical performance, making it less susceptible to microdamage than models with larger DoE values.

INTERPRETATION

These findings enhance the limited knowledge of the influence of the lacunar shape on microdamage and contribute to a better understanding of bone biomechanics.

Cortical Bone Modifications after Radiotherapy: Cortex Porosity and Osteonal Changes Evaluated Over Time

Aiming to evaluate cortical bone microarchitecture and osteonal morphology after irradiation, twelve male New Zealand rabbits were used. The animals were divided: control group (no radiation-NIr); and 3 irradiated groups, sacrificed after: 7 (Ir7d); 14 (Ir14d) and 21 (Ir21d) days. A single radiation dose of 30 Gy was used. Computed microtomography analyzed the cortical microarchitecture: cortical thickness (CtTh), bone volume (BV), total porosity (Ct.Po), intracortical porosity (CtPo-cl), channel/pore number (Po.N), fractal dimension (FD) and degree of anisotropy (Ct.DA). After scan, osteonal morphology was histologically assessed by means: area and perimeter of the osteons (O.Ar; O.p) and of the Haversian canals (C.Ar; C.p). Microtomographic analysis were performed by ANOVA, followed by Tukey and Dunnet tests. Osteon morphology analyses were performed by Kruskal-Wallis, and test Dunn's. Cortical thickness was significant difference (p<0.010) between the NIr and irradiated groups, with thicker cortex at Ir7d (1.15±0.09). The intracortical porosity revealed significant difference (p<0.001) between irradiated groups and NIr, with lower value for Ir7d (0.29±0.09). Bone volume was lower in Ir14d compared to control. Area and perimeter of the osteons were statistically different (p<0.0001) between NIr and Ir7d. Haversian canals also revealed lower values (p<0.0001) in Ir7d (80.57±9.3; 31.63±6.5) compared to NIr and irradiated groups. Cortical microarchitecture was affected by radiation, and the effects appear to be time-dependent, mostly regarding the osteons morphology at the initial days. Cortex structure in Ir21d revealed similarities to control suggesting that microarchitecture resembles normal condition after a period.

Finite element analysis on multi-toughening mechanism of microstructure of osteon

The toughening mechanism of cortical bone is closely related to its hierarchical microstructure. Osteon is the most important microstructure of cortical bone. Therefore, it is very important to study the toughening mechanism of the microstructure of osteon. There are three main kinds of cracks in cortical bone: external crack of osteon, internal radial cracks of osteon and microporous damage cracks. Numerical models for these three kinds of cracks are established by XFEM and the progressive damage approach, respectively. The multi-toughening mechanisms of microstructure of osteon are found. The cement line on the outside of osteon is its first toughening mechanism, which can make the crack deflection and improve the fracture resistance of osteon. The resistance of cement line to fracture increases with the decrease of the strength and the increase of the thickness. The second toughening mechanism is elliptical osteocyte lacunae, which can attract the crack into the elliptical lacunae and cause stress redistribution to prevent the crack propagation. The annularly elliptical lacuna structure is an optimized arrangement and shape of microstructure, which is the third toughening mechanism of osteon. This microstructure can determine the location of the crack initiation and make the microcracks propagate along the annular direction rather than penetrating into the haversian cannal to protect the integrity of the osteon. The study of these toughening mechanisms provides new ideas for the research and design of synthetic composite structures.

Study on the biomechanical responses of the loaded bone in macroscale and mesoscale by multiscale poroelastic FE analysis

BACKGROUND

Bone is a hierarchically structured composite material, and different hierarchical levels exhibit diverse material properties and functions. The stress and strain distribution and fluid flow in bone play an important role in the realization of mechanotransduction and bone remodeling.

METHODS

To investigate the mechanotransduction and fluid behaviors in loaded bone, a multiscale method was developed. Based on poroelastic theory, we established the theoretical and FE model of a segment bone to provide basis for researching more complex bone model. The COMSOL Multiphysics software was used to establish different scales of bone models, and the properties of mechanical and fluid behaviors in each scale were investigated.

RESULTS

FE results correlated very well with analytical in macroscopic scale, and the results for the mesoscopic models were about less than 2% different compared to that in the macro-mesoscale models, verifying the correctness of the modeling. In macro-mesoscale, results demonstrated that variations in fluid pressure (FP), fluid velocity (FV), von Mises stress (VMS), and maximum principal strain (MPS) in the position of endosteum, periosteum, osteon, and interstitial bone and these variations can be considerable (up to 10, 8, 4 and 3.5 times difference in maximum FP, FV, VMS, and MPS between the highest and the lowest regions, respectively). With the changing of Young's modulus (E) in each osteon lamella, the strain and stress concentration occurred in different positions and given rise to microscale spatial variations in the fluid pressure field. The heterogeneous distribution of lacunar-canalicular permeability (klcp) in each osteon lamella had various influence on the FP and FV, but had little effect on VMS and MPS.

CONCLUSION

Based on the idealized model presented in this article, the presence of endosteum and periosteum has an important influence on the fluid flow in bone. With the hypothetical parameter values in osteon lamellae, the bone material parameters have effect on the propagation of stress and fluid flow in bone. The model can also incorporate alternative material parameters obtained from different individuals. The suggested method is expected to provide dependable biological information for better understanding the bone mechanotransduction and signal transduction.

Functional disuse initiates medullary endosteal micro-architectural impairment in cortical bone characterized by nanoindentation

In this study, we evaluated the effect of functional disuse-induced bone remodeling on its mechanical properties, individually at periosteum and medullary endosteum regions of the cortical bone. Left middle tibiae were obtained from 5-month-old female Sprague-Dawley rats for the baseline control as well as hindlimb suspended (disuse) groups. Micro-nano-mechanical elastic moduli (at lateral region) was evaluated along axial (Z), circumferential (C) and radial (R) orientations using nanoindentation. Results indicated an anisotropic microstructure with axial orientation having the highest and radial orientation with the lowest moduli at periosteum and medullary endosteum for both baseline control as well as disuse groups. Between the groups: at periosteum, an insignificant difference was evaluated for each of the orientations (p > 0.05) and at endosteum, a significant decrease of elastic moduli in the radial (p < 0.0001), circumferential (p < 0.001) and statistically insignificant difference in axial (p > 0.05) orientation. For the moduli ratios between groups: at periosteum, only significant difference in the Z/R (p < 0.05) anisotropy ratio, whereas at endosteum, a statistically significant difference in Z/C (p < 0.001), and Z/R (p < 0.001), as well as C/R (p < 0.05) anisotropy ratios, was evaluated. The results suggested initial bone remodeling impaired bone micro-architecture predominantly at the medullary endosteum with possible alterations in the geometric orientations of collagen and mineral phases inside the bone. The findings could be significant for studying the mechanotransduction pathways involved in maintaining the bone micro-architecture and possibly have high clinical significance for drug use against impairment from functional disuse.

Application of the Johnson-Cook plasticity model in the finite element simulations of the nanoindentation of the cortical bone

The mechanical behavior of the cortical bone in nanoindentation is a complicated mechanical problem. The finite element analysis has commonly been assumed to be the most appropriate approach to this issue. One significant problem in nanoindentation modeling of the elastic-plastic materials is pile-up deformation, which is not observed in cortical bone nanoindentation testing. This phenomenon depends on the work-hardening of materials; it doesn't occur for work-hardening materials, which suggests that the cortical bone could be considered as a work-hardening material. Furthermore, in a recent study [59], a plastic hardening until failure was observed on the micro-scale of a dry ovine osteonal bone samples subjected to micropillar compression. The purpose of the current study was to apply an isotropic hardening model in the finite element simulations of the nanoindentation of the cortical bone to predict its mechanical behavior. The Johnson-Cook (JC) model was chosen as the constitutive model. The finite element modeling in combination with numerical optimization was used to identify the unknown material constants and then the finite element solutions were compared to the experimental results. A good agreement of the numerical curves with the target loading curves was found and no pile-up was predicted. A Design Of Experiments (DOE) approach was performed to evaluate the linear effects of the material constants on the mechanical response of the material. The strain hardening modulus and the strain hardening exponent were the most influential parameters. While a positive effect was noticed with the Young's modulus, the initial yield stress and the strain hardening modulus, an opposite effect was found with the Poisson's ratio and the strain hardening exponent. Finally, the JC model showed a good capability to describe the elastoplastic behavior of the cortical bone.

Role of porosity and matrix behavior on compressive fracture of Haversian bone using random spring network model

Haversian remodeling is known to result in improved resistance to compressive fracture in healthy cortical bone. Here, we examine the individual roles of the mean porosity, structure of the network of pores and remodeled bone matrix properties in the fracture behavior of Haversian bone. The detailed structure of porosity network is obtained both pre- and post-testing of dry cubical bone samples using micro-Computed Tomography. Based on the periodicity in the features of porosity along tangential direction, we develop a two dimensional porosity-based random spring network model for Haversian bone. The model is shown to capture well the macroscopic response and reproduce the avalanche statistics similar to recently reported experiments on porcine bone. The predictions suggest that at the millimeter scale, the remodeled bone matrix of Haversian bone is less stiff but tougher than that of plexiform/primary bone.

Splitting fracture in bovine bone using a porosity-based spring network model

We examine the specific role of the structure of the network of pores in plexiform bone in its fracture behaviour under compression. Computed tomography scan images of the sample pre- and post-compressive failure show the existence of weak planes formed by aligned thin long pores extending through the length. We show that the physics of the fracture process is captured by a two-dimensional random spring network model that reproduces well the macroscopic response and qualitative features of fracture paths obtained experimentally, as well as avalanche statistics seen in recent experiments on porcine bone.

Haversian microstructure in bovine femoral cortices: An adaptation for improved compressive strength

Microstructural variations in bovine femoral cortices and its possible implications for the bone's mechanical behavior are characterized for a mature and a young bovine femur. Histological examination at several locations shows the presence of Haversian systems to be largely confined to the posterior region of any cross-section. Haversian bone is shown to have higher compressive strength than the non-Haversian primary bone present in the corresponding anterior regions. The anatomical variation in the compressive strength along diaphysis is found to correlate strongly with the Haversian density. Based on the differences in the failure surfaces observed from compressive failure, it is argued that the presence of Haversian systems plays a role in deflection of crack path, leading to non-prismatic failure surfaces. As biomaterials, such as bone cement and implants, closely interact with bone material, the structure-property relation established here can provide a basis for better design of future biomaterials.

Elastic Properties of Human Osteon and Osteonal Lamella Computed by a Bidirectional Micromechanical Model and Validated by Nanoindentation

Knowledge of the anisotropic elastic properties of osteon and osteonal lamellae provides a better understanding of various pathophysiological conditions, such as aging, osteoporosis, osteoarthritis, and other degenerative diseases. For this reason, it is important to investigate and understand the elasticity of cortical bone. We created a bidirectional micromechanical model based on inverse homogenization for predicting the elastic properties of osteon and osteonal lamellae of cortical bone. The shape, the dimensions, and the curvature of osteon and osteonal lamellae are described by appropriately chosen curvilinear coordinate systems, so that the model operates close to the real morphology of these bone components. The model was used to calculate nine orthotropic elastic constants of osteonal lamellae. The input values have the elastic properties of a single osteon. We also expressed the dependence of the elastic properties of the lamellae on the angle of orientation. To validate the model, we performed nanoindentation tests on several osteonal lamellae. We compared the experimental results with the calculated results, and there was good agreement between them. The inverted model was used to calculate the elastic properties of a single osteon, where the input values are the elastic constants of osteonal lamellae. These calculations reveal that the model can be used in both directions of homogenization, i.e., direct homogenization and also inverse homogenization. The model described here can provide either the unknown elastic properties of a single lamella from the known elastic properties at the level of a single osteon, or the unknown elastic properties of a single osteon from the known elastic properties at the level of a single lamella.

Nanomechanical mapping of bone tissue regenerated by magnetic scaffolds

Nanoindentation can provide new insights on the maturity stage of regenerating bone. The aim of the present study was the evaluation of the nanomechanical properties of newly-formed bone tissue at 4 weeks from the implantation of permanent magnets and magnetic scaffolds in the trabecular bone of rabbit femoral condyles. Three different groups have been investigated: MAG-A (NdFeB magnet + apatite/collagen scaffold with magnetic nanoparticles directly nucleated on the collagen fibers during scaffold synthesis); MAG-B (NdFeB magnet + apatite/collagen scaffold later infiltrated with magnetic nanoparticles) and MAG (NdFeB magnet). The mechanical properties of different-maturity bone tissues, i.e. newly-formed immature, newly-formed mature and native trabecular bone have been evaluated for the three groups. Contingent correlations between elastic modulus and hardness of immature, mature and native bone have been examined and discussed, as well as the efficacy of the adopted regeneration method in terms of "mechanical gap" between newly-formed and native bone tissue. The results showed that MAG-B group provided regenerated bone tissue with mechanical properties closer to that of native bone compared to MAG-A or MAG groups after 4 weeks from implantation. Further, whereas the mechanical properties of newly-formed immature and mature bone were found to be fairly good correlated, no correlation was detected between immature or mature bone and native bone. The reported results evidence the efficacy of nanoindentation tests for the investigation of the maturity of newly-formed bone not accessible through conventional analyses.

Cyclic cryopreservation affects the nanoscale material properties of trabecular bone

Tissues such as bone are often stored via freezing, or cryopreservation. During an experimental protocol, bone may be frozen and thawed a number of times. For whole bone, the mechanical properties (strength and modulus) do not significantly change throughout five freeze-thaw cycles. Material properties at the trabecular and lamellar scales are distinct from whole bone properties, thus the impact of freeze-thaw cycling at this scale is unknown. To address this, the effect of repeated freezing on viscoelastic material properties of trabecular bone was quantified via dynamic nanoindentation. Vertebrae from five cervine spines (1.5-year-old, male) were semi-randomly assigned, three-to-a-cycle, to 0-10 freeze-thaw cycles. After freeze-thaw cycling, the vertebrae were dissected, prepared and tested. ANOVA (factors cycle, frequency, and donor) on storage modulus, loss modulus, and loss tangent, were conducted. Results revealed significant changes between cycles for all material properties for most cycles, no significant difference across most of the dynamic range, and significant differences between some donors. Regression analysis showed a moderate positive correlation between cycles and material property for loss modulus and loss tangent, and weak negative correlation for storage modulus, all correlations were significant. These results indicate that not only is elasticity unpredictably altered, but also that damping and viscoelasticity tend to increase with additional freeze-thaw cycling.

Micro-morphological properties of osteons reveal changes in cortical bone stability during aging, osteoporosis, and bisphosphonate treatment in women

SUMMARY

We analyzed morphological characteristics of osteons along with the geometrical indices of individual osteonal mechanical stability in young, healthy aged, untreated osteoporotic, and bisphosphonate-treated osteoporotic women. Our study revealed significant intergroup differences in osteonal morphology and osteocyte lacunae indicating different remodeling patterns with implications for fracture susceptibility.

INTRODUCTION

Bone remodeling is the key process in bone structural reorganization, and its alterations lead to changes in bone mechanical strength. Since osteons reflect different bone remodeling patterns, we hypothesize that the femoral cortices of females under miscellaneous age, disease and treatment conditions will display distinct osteonal morphology and osteocyte lacunar numbers along with different mechanical properties.

METHODS

The specimens used in this study were collected at autopsy from 35 female donors (young group, n = 6, age 32 ± 8 years; aged group, n = 10, age 79 ± 9 years; osteoporosis group, n = 10, age 81 ± 9 years; and bisphosphonate group, n = 9, age 81 ± 7 years). Von Kossa-modified stained femoral proximal diaphyseal sections were evaluated for osteonal morphometric parameters and osteocyte lacunar data. Geometrical indices of osteonal cross-sections were calculated to assess the mechanical stability of individual osteons, in terms of their resistance to compression, bending, and buckling.

RESULTS

The morphological assessment of osteons and quantification of their osteocyte lacunae revealed significant differences between the young, aged, osteoporosis and bisphosphonate-treated groups. Calculated osteonal geometric indices provided estimates of the individual osteons' resistance to compression, bending and buckling based on their size. In particular, the osteons in the bisphosphonate-treated group presented improved osteonal geometry along with increased numbers of osteocyte lacunae that had been formerly impaired due to aging and osteoporosis.

CONCLUSIONS

The data derived from osteons (as the basic structural units of the cortical bone) in different skeletal conditions can be employed to highlight structural factors contributing to the fracture susceptibility of various groups of individuals.

Investigation of the three-dimensional orientation of mineralized collagen fibrils in human lamellar bone using synchrotron X-ray phase nano-tomography

We investigate the three-dimensional (3-D) organization of mineralized collagen fibrils in human cortical bone based on synchrotron X-ray phase nano-tomography images. In lamellar bone the collagen fibrils are assumed to have a plywood-like arrangement, but due to experimental limitations the 3-D fibril structure has only been deduced from section surfaces so far and the findings have been controversial. Breakthroughs in synchrotron tomographic imaging have given access to direct 3-D information on the bone structure at the nanoscale level. Using an autocorrelation-based orientation measure we confirm that the fibrils are unidirectional in quasi-planes of sub-lamellae and find two specific dominant patterns, oscillating and twisted plywoods coexisting in a single osteon. Both patterns exhibit smooth orientation changes between adjacent quasi-planes. Moreover, we find that the periodic changes in collagen fibril orientation are independent of fluctuations in local mass density. These data improve our understanding of the lamellar arrangement in bone and allow more detailed investigations of structure-function relationships at this scale, providing templates for bio-inspired materials. The presented methodology can be applied to non-destructive 3-D characterization of the sub-micron scale structure of other natural and artificial mineralized biomaterials.

Osteonal lamellae elementary units: lamellar microstructure, curvature and mechanical properties

The mechanical and structural properties of the sublayers of osteonal lamellae were studied. Young's modulus (E) of adjacent individual lamellae was measured by nanoindentation of parallel slices every 1-3 μm, in planes parallel and perpendicular to the osteon axis (OA). In planes parallel to the OA, the modulus of a lamella could vary significantly between sequential slices. Significant modulus variations were also sometimes found on opposing sides of the osteonal canal for the same lamella. These results are rationalized by considerations involving the microstructural organization of the collagen fibrils in the lamellae. Scanning electron microscope imaging of freeze fractured surfaces revealed that the substructure of a single lamella can vary significantly on the opposing sides of the osteonal axis. Using a serial surface view method, parallel planes were exposed every 8-10 nm using a dual-beam microscope. Analysis of the orientations of fibrils revealed that the structure is rotated plywood like, consisting of unidirectional sublayers of fibrils of several orientations, with occasional randomly oriented sublayers. The dependence of the measured mechanical properties of the lamellae on the indentation location may be explained by the observed structure, as well as by the curvature of the osteonal lamellae through simple geometrical-structural considerations. Mechanical advantages arising from the curved laminate structure are discussed.

Micro-structure and mechanical properties of the turtle carapace as a biological composite shield

Turtle shell is a multi-scale bio-composite in which the components are arranged in various spatial patterns, leading to an unusually strong and durable structure. The keratin-coated dorsal shell, termed the carapace, exhibits a flat bone, sandwich-like structure made up of two exterior cortices enclosing a cancellous interior. This unique structure was developed by nature to protect the reptile from predator attacks by sustaining impact loads and dissipating energy. In the present study we attempt to correlate the micro-scale architecture with the mechanical properties of the carapace sub-regions of the red-eared slider turtle. The microscopic structural features were examined by scanning electron microscopy and micro-computed tomography. Nanoindentation tests were performed under dry and wet conditions on orthogonal anatomical planes to evaluate the elastic modulus and hardness of the various carapace sub-regions. The mineral content was also measured in the different regions of the carapace. Consequently, we discuss the influence of hydration on the carapace sub-regions and the contribution of each sub-region to the overall mechanical resistance of the assemblage.

Microfibril orientation dominates the microelastic properties of human bone tissue at the lamellar length scale

The elastic properties of bone tissue determine the biomechanical behavior of bone at the organ level. It is now widely accepted that the nanoscale structure of bone plays an important role to determine the elastic properties at the tissue level. Hence, in addition to the mineral density, the structure and organization of the mineral nanoparticles and of the collagen microfibrils appear as potential key factors governing the elasticity. Many studies exist on the role of the organization of collagen microfibril and mineral nanocrystals in strongly remodeled bone. However, there is no direct experimental proof to support the theoretical calculations. Here, we provide such evidence through a novel approach combining several high resolution imaging techniques: scanning acoustic microscopy, quantitative scanning small-Angle X-ray scattering imaging and synchrotron radiation computed microtomography. We find that the periodic modulations of elasticity across osteonal bone are essentially determined by the orientation of the mineral nanoparticles and to a lesser extent only by the particle size and density. Based on the strong correlation between the orientation of the mineral nanoparticles and the collagen molecules, we conclude that the microfibril orientation is the main determinant of the observed undulations of microelastic properties in regions of constant mineralization in osteonal lamellar bone. This multimodal approach could be applied to a much broader range of fibrous biological materials for the purpose of biomimetic technologies.