Orientation of bone mineral and its role in the anisotropic mechanical properties of bone--transverse anisotropy

An orientation of hydroxyapatite (HAP) crystals in bovine femur mineral was investigated by means of X-ray pole figure analysis (XPFA). It was found that the c-axis of HAP generally orients parallel to the longitudinal axis of bone (bone axis) and a significant amount of c-axis was oriented in other directions, in particular, perpendicular to the bone axis. Comparing these results with those of the small angle X-ray scattering (SAXS) investigation by Matsushima et al. (Jap. J. appl. Phys. 21, 186-189, 1982) at least two types of morphology of bone mineral were found; rod like bone mineral having the c-axis of HAP crystal parallel to the longitudinal axis of the rod and that having the c-axis not parallel, in a particular case, perpendicular to its longitudinal axis. Transverse anisotropy in mechanical properties of bone was reproduced by the estimation of Young's moduli by using the structural results mainly from XPFA. It is concluded that the anisotropy in mechanical properties of bone is well explained by taking account of the non-longitudinal (off-bone) axial distribution of orientation of bone mineral.     

Tensile yield in compact bone is determined by strain, post-yield behaviour by mineral content

Compact bone specimens from many species were examined to determine the relationships, in tension, between mineral content, Young's modulus, yield stress, yield strain, post-yield stress, post-yield strain, ultimate stress, ultimate strain and work under the stress-strain curve. Yield strain varied much less than the post-yield strain, and yield stress was strongly dependent on Young's modulus. Mineral content was a rather poor predictor of yield stress. However, post-yield events were predicted better by mineral (calcium) content than by Young's modulus. The greater the mineral content the less the post-yield work under the curve and the less the increase in post-yield stress and strain. The findings are compared with those of Les et al. who compressed specimens from equine metacarpals. Where they can be compared, the results are consistent with each other.     

Spatial orientation in bone samples and Young's modulus

Bone mass is the most important determinant of the mechanical strength of bones, and spatial structure is the second. In general, the spatial structure and mechanical properties of bones such as the breaking strength are direction dependent. The mean intercept length (MIL) and line frequency deviation (LFD) are two methods for quantifying directional aspects of the spatial structure of bone. Young's modulus is commonly used to describe the stiffness of bone, which is also a direction-dependent mechanical property. The aim of this article is to investigate the relation between MIL and LFD on one hand and Young's modulus on the other. From 11 human mandibular condyles, 44 samples were taken and scanned with high-resolution computer tomography equipment (micro-CT). For each sample the MIL and LFD were determined in 72602 directions distributed evenly in 3D space. In the same directions Young's modulus was determined by means of the stiffness tensor that had been determined for each sample by finite element analysis. To investigate the relation between the MIL and LFD on one hand and Young's modulus on the other, multiple regression was used. On average the MIL accounted for 69% of the variance in Young's modulus in the 44 samples and the LFD accounted for 72%. The average percentage of variance accounted for increased to 80% when the MIL was combined with the LFD to predict Young's modulus. Obviously MIL and LFD to some extent are complementary with respect to predicting Young's modulus. It is known that directional plots of the MIL tend to be ellipses or ellipsoids. It is speculated that ellipsoids are not always sufficient to describe Young's modulus of a bone sample and that the LFD partly compensates for this.     

The underestimation of Young's modulus in compressive testing of cancellous bone specimens

In order to determine the accuracy of measurements of Young's modulus of cancellous bone by conventional compression testing, two independent strain measurements were made simultaneously during non-destructive uniaxial compression to 0.8% strain of rectangular specimens (n = 18). Strain was measured by an extensometer attached to the compression anvils close to the specimen and by an optical system covering the central half of the specimens. Mean Young's modulus determined by the extensometer technique was 689 MPa, but was 871 MPa when determined by the optical technique (mean difference = 182 MPa, SED = 50 MPa, p less than 0.002). Uneven strain distribution due to lack of support of cut vertical trabeculae at the anvil-specimen interface is believed to be causing the underestimation of Young's modulus measured by the extensometer technique. The influence of friction at the specimen-anvil interface was studied by performing a finite element analysis. It is concluded that Young's modulus of specimens of the chosen geometry on average is underestimated by about 20% by conventional compressing testing. The underestimation seems not to be dependent upon specimen density.     

The elastic properties of trabecular and cortical bone tissues are similar: results from two microscopic measurement techniques

Acoustic microscopy (30-60 microm resolution) and nanoindentation (1-5 microm resolution) are techniques that can be used to evaluate the elastic properties of human bone at a microstructural level. The goals of the current study were (1) to measure and compare the Young's moduli of trabecular and cortical bone tissues from a common human donor, and (2) to compare the Young's moduli of bone tissue measured using acoustic microscopy to those measured using nanoindentation. The Young's modulus of cortical bone in the longitudinal direction was about 40% greater than (p<0.01) the Young's modulus in the transverse direction. The Young's modulus of trabecular bone tissue was slightly higher than the transverse Young's modulus of cortical bone, but substantially lower than the longitudinal Young's modulus of cortical bone. These findings were consistent for both measurement methods and suggest that elasticity of trabecular tissue is within the range of that of cortical bone tissue. The calculation of Young's modulus using nanoindentation assumes that the material is elastically isotropic. The current results, i.e., the average anisotropy ratio (E(L)/E(T)) for cortical bone determined by nanoindentation was similar to that determined by the acoustic microscope, suggest that this assumption does not limit nanoindentation as a technique for measurement of Young's modulus in anisotropic bone.     

On the relationship between the microstructure of bone and its mechanical stiffness.

A recent study of bone structure shows that the plate-shaped carbonate apatite crystals in individual lamellae are arranged in layers across the lamellae, and that the orientation of these layers are different in alternate lamellae. Based on these findings, a new micromechanical model for the Young's modulus of bone is proposed, which accounts for the anisotropy and geometrical characteristics of the material. The model incorporates the platelet-like geometry of the basic reinforcing unit, the presence of alternating thin and thick lamellae, and the orientations of the crystal platelets in the lamellae. The thin and thick lamellae are modeled as orthotropic composite layers made up of thin rectangular apatite platelets within a collagen matrix, and classical orthotropic elasticity theory is used to calculate the Young's modulus of the lamellae. Bone is viewed as an assembly of such orthotropic lamellae bent into cylindrical structures, and having a constant, alternating angle between successive lamellae. The micromechanical model employs a modified rule-of-mixtures to account for the two types of lamellae. The model provides a curve similar to the published experimental data on the angular dependence of Young's modulus, including a local maximum at an angle between 0 and 90 degrees. A rigorous testing of the model awaits additional experimental data.     

Effect of bone mineral content on the tensile properties of cortical bone: experiments and theory

The effect of mineral volume fraction on the tensile mechanical properties of cortical bone tissue is investigated by theoretical and experimental means. The mineral content of plexiform, bovine bone was lowered by 18% and 29% by immersion in fluoride solutions for 3 days and 12 days, respectively. The elastic modulus, yield strength and ultimate strength of bone tissue decreased, while the ultimate strain increased with a decrease in mineral content. The mechanical behavior of bone tissue was modeled by using a micromechanical shear lag theory consisting of overlapped mineral platelets reinforcing the organic matrix. The decrease in yield stress, by the 0.002 offset method, of the fluoride treated bones were matched in the theoretical curves by lowering the shear yield stress of the organic matrix. The failure criterion used was based on failure stresses determined from a failure envelope (Mohr's circle), which was constructed using experimental data. It was found that the model predictions of elastic modulus got worse with a decrease in mineral content (being 7.9%, 17.2% and 33.0% higher for the control, 3-day and 12-day fluoride-treated bones). As a result, the developed theory could not fully predict the yield strain of bones with lowered mineral content, being 12.9% and 21.7% lower than the experimental values. The predicted ultimate stresses of the bone tissues with lower mineral contents were within +/- 10% of the experimental values while the ultimate strains were 12.7% and 26.3% lower than the experimental values. Although the model developed in this study did not take into account the presence of hierarchical structures, voids, orientation of collagen molecules and micro cracks, it still indicated that the mechanical properties of the organic matrix depend on bone mineral content.     

Incompatible mechanical properties in compact bone

The covariation of a number of mechanical of properties, and some physical characteristics, of compact bones from a wide range of bones were examined. Young's modulus was well predicted by a combination of mineral content and porosity. Increasing Young's modulus was associated with: increasing stress at yield, increasing bending strength, and a somewhat higher resilience, tensile strength and fatigue strength. Contrarily, in the post-yield region a higher Young's modulus (and more clearly, a higher mineral content) was associated with: a reduced work to fracture in tension, a reduced impact strength and an increased notch sensitivity in impact. Increasing porosity is associated with deleterious effects in the pre-yield region, but has little effect in the post-yield region. Bone, like many other materials, is unable to have good qualities in both the pre- and post-yield regions. Since an increase in mineral or Young's modulus is more potent, that is deleterious, in the post-yield than it is advantageous in the pre-yield region, it is likely that mineral content will be selected to be slightly lower than would be the case if it were equally potent in both regions. As is usual in biology, different adaptive extremes are incompatible.     

Micromechanical response of mineral and collagen phases in bone

We report the first simultaneous quantification of Young's modulus in the separate material phases of bone: collagen and carbonated hydroxyapatite. High-energy X-ray scattering and in situ loading revealed macroscopic, mineral, and collagen Young's moduli (90% confidence limit) for a canine fibula equaled 24.7(0.2) GPa, 38.2(0.5) GPa {for 00.4 and 43.6(1.4) GPa for 22.2}, and 18(1.2) GPa, respectively. The mineral contained compressive residual stresses on the order of -60 to -80 MPa before loading and had a stress enhancement (ratio of internal to applied stress) between 2.0 and 2.3. The diffraction peak width increased with increasing applied stress, mainly along the bone's longitudinal direction, and peak widths returned to pre-deformation values when load was removed. In a second fibula section from the same animal, the mineral's internal stress changed from -50 MPa (22.2 reflection) to -75 MPa (00.4) just after removal from formalin to -10 MPa after eight hours immersion in phosphate-buffered saline; the corresponding change in collagen D-spacing DeltaD/D equaled 4.2x10(-3).     

Evaluation of bone plate with low-stiffness material in terms of stress distribution

The aim of this study is to evaluate a newly developed bone plate with low-stiffness material in terms of stress distribution. In this numerical study, 3D finite element models of the bone plate with low-stiffness material and traditional bone plates made of stainless steel and Ti alloy have been developed by using the ANSYS software. Stress analyses have been carried out for all three models under the same loading and boundary conditions. Compressive stresses occurring in the intact portion of the bone (tibia) and at the fractured interface at different stages of bone healing have been investigated for all three types of bone-plate systems. The results obtained have been compared and presented in graphs. It has been seen that the bone plate with low-stiffness material offers less stress-shielding to the bone, providing a higher compressive stress at the fractured interface to induce accelerated healing in comparison with Ti alloy and stainless-steel bone plate. In addition, the effects of low-stiffness materials with different Young's modulus on stress distribution at the fractured interface have been investigated in the newly developed bone-plate system. The results showed that when a certain value of Young's modulus of low-stiffness material is exceeded, increase in stiffness of the bone plate does not occur to a large extent and stress distributions and micro-motions at the fractured interface do not change considerably.     

Nanostructure and elastic modulus of single trabecula in bovine cancellous bone

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

The dependence of transversely isotropic elasticity of human femoral cortical bone on porosity

The objective of this study was to examine the dependence of the elastic properties of cortical bone as a transversely isotropic material on its porosity. The longitudinal Young's modulus, transverse Young's modulus, longitudinal shear modulus, transverse shear modulus, and longitudinal Poisson's ratio of cortical bone were determined from eighteen groups of longitudinal and transverse specimens using tensile and torsional tests on a servo-hydraulic material testing system. These cylindrical waisted specimens of cortical bone were harvested from the middle diaphysis of three pairs of human femora. The porosity of these specimens was assessed by means of histology. Our study demonstrated that the longitudinal Young's and shear moduli of human femoral cortical bone were significantly (p<0.01) negatively correlated with the porosity of cortical bone. Conversely, the elastic properties in the transverse direction did not have statistically significant correlations with the porosity of cortical bone. As a result, the transverse elastic properties of cortical bone were less sensitive to changes in porosity than those in the longitudinal direction. Additionally, the anisotropic ratios of cortical bone elasticity were found to be significantly (p<0.01) negatively correlated with its porosity, indicating that cortical bone tended to become more isotropic when its porosity increased. These results may help a number of researchers develop more accurate micromechanics models of cortical bone.     

Elastic properties of woven bone: effect of mineral content and collagen fibrils orientation

Woven bone is a type of tissue that forms mainly during fracture healing or fetal bone development. Its microstructure can be modeled as a composite with a matrix of mineral (hydroxyapatite) and inclusions of collagen fibrils with a more or less random orientation. In the present study, its elastic properties were estimated as a function of composition (degree of mineralization) and fibril orientation. A self-consistent homogenization scheme considering randomness of inclusions’ orientation was used for this purpose. Lacuno-canalicular porosity in the form of periodically distributed void inclusions was also considered. Assuming collagen fibrils to be uniformly oriented in all directions led to an isotropic tissue with a Young’s modulus \(E = 1.90\) GPa, which is of the same order of magnitude as that of woven bone in fracture calluses. By contrast, assuming fibrils to have a preferential orientation resulted in a Young’s modulus in the preferential direction of 9–16 GPa depending on the mineral content of the tissue. These results are consistent with experimental evidence for woven bone in foetuses, where collagen fibrils are aligned to a certain extent.     

The effect of porosity and mineral content on the Young's modulus of elasticity of compact bone

The Young's modulus of elasticity, the calcium content and the volume fraction (1-porosity) of 23 tension specimens and 80 bending specimens, taken from compact bone of 18 species of mammal, bird and reptile, were determined. There was a strong positive relationship between Young's modulus and both calcium content and volume fraction. A power law model fits the data better than a linear model. Young's modulus has a roughly cubic relationship with both calcium content and volume fraction. Over 80% of the total variation in Young's modulus in this data set is explained by these two variables.     

Apparent Young's modulus of vertebral cortico-cancellous bone specimens

Up to now, due to cortical thickness and imaging resolution, it is not possible to derive subject-specific mechanical properties on the 'vertebral shell' from imaging modalities applicable in vivo. As a first step, the goal of this study was to assess the apparent Young's modulus of vertebral cortico-cancellous bone specimens using an inverse method. A total of 22 cortico-cancellous specimens were harvested from 22 vertebral bodies. All specimens were tested in compression until failure. To compute the apparent Young's modulus of the specimen from the inverse method, the boundary conditions of the biomechanical experiments were faithfully reproduced in a finite element model (FEM), and an optimisation routine was used. The results showed a mean of the apparent Young's modulus of 374?±?208?MPa, ranging from 87 to 791?MPa. By computing an apparent Young's modulus of a cortico-cancellous medium, this study gives mechanical data for an FEM of an entire vertebra including an external shell combining both bone tissues.     

Determination of mechanical properties of human femoral cortical bone by the Hopkinson bar stress technique

The Hopkinson bar stress technique and a universal testing machine (Instron 1125) have been used to investigate the dynamic and static mechanical properties of cortical bone taken from a human femur respectively. We found that the average dynamic Young's modulus value (Ed = 19.9 GPa) to be 23% higher than the average static Young's modulus value (Ed = 16.2 GPa). Furthermore, the Poisson's ratio did not exhibit any significant variation for the two different types of loading. No difference was observed between the values of the dynamic Young's modulus in tension and those found in compression. A comparison was made of the results of this study with those found by other researchers using different techniques, such as ultrasonics, and it was found that they agree well with most of the results of previous studies. Finally, the viscosity for cortical bone found in this study correlates with viscosity reported by Tennyson et al. [Expl Mech. 12, 502-507 (1972)] for ten days post mortem age specimens.     

Topographical distribution of trabecular bone strength in the human os calcanei

Mechanical properties of twenty human os calcanei were determined by uniaxial compression testing of bone specimens from facies articularis talaris posterior, facies articularis cuboidea, and tuber calcanei. Specimens were taken oriented perpendicular to the planes of the facies articularis, and in tuber along the presumed loading axis throughout the gait cycle. Young's modulus and strength at facies articularis cuboidea and facies articularis talaris posterior were about three times those at the tuber calcanei. The variation of the relationship between Young's modulus and apparent density indicated differences in the orientation of the trabecula, in relation to the direction of evaluation between these locations. A more detailed analysis of the topographical variation of strength within each location was made using penetration testing of a further nineteen specimens. The results of both types of measurements indicated that the major part of the load during walking is carried by facies articularis talaris posterior and facies articularis cuboidea.     

Trabecular bone's mechanical properties are affected by its non-uniform mineral distribution

The bone remodeling process takes place at the surface of trabeculae and results in a non-uniform mineral distribution. This will affect the mechanical properties of cancellous bone, because the properties of bone tissue depend on its mineral content. We investigated how large this effect is by simulating several non-uniform mineral distributions in 3D finite element models of human trabecular bone and calculating the apparent stiffness of these models. In the ‘linear model’ we assumed a linear relation between mineral content and Young's modulus of the tissue. In the ‘exponential model’ we included an empirical exponential relation in the model. When the linear model was used the mineral distribution slightly changed the apparent stiffness, the difference varied between an 8% decrease and a 4% increase compared to the uniform model with the same BMD. The exponential model resulted in up to 20% increased apparent stiffness in the main load-bearing direction. A thin less mineralized surface layer (28 μm) and highly mineralized interstitial bone (mimicking mineralization resulting from anti-resorptive treatment) resulted in the highest stiffness. This could explain large reductions in fracture risk resulting from small increases in BMD. The non-uniform mineral distribution could also explain why bone tissue stiffness determined using nano-indentation is usually higher than finite element (FE)-determined stiffness. We conclude that the non-uniform mineral distribution in trabeculae does affect the mechanical properties of cancellous bone and that the tissue stiffness determined using FE-modeling could be improved by including detailed information about mineral distribution in trabeculae in the models.     

Anisotropic viscoelastic properties of cortical bone

Relaxation Young's modulus of cortical bone was investigated for two different directions with respect to the longitudinal axis of bone (bone axis, BA): the modulus parallel (P) and normal (N) to the BA. The relaxation modulus was analyzed by fitting to the empirical equation previously proposed for cortical bones, i.e., a linear combination of two Kohlraush-Williams-Watts (KWW) functions (Iyo et al., 2003. Biorheology, submitted): E(t)=E0 (A1 exp[-(t/tau1)beta]+(1-A1) exp[-(t/tau2)gamma]), [0 < A1, beta, gamma < 1], where E0 is the initial modulus value E0. Tau1 and tau2(>tau1) are characteristic times of the relaxation, A1 is the fractional contribution of the fast relaxation (KWW1 process) to the whole relaxation process, and beta and gamma are parameters describing the shape of the relaxation modulus. In both P and N samples, the relaxation modulus was described well by the empirical equation. The KWW1 process of a P sample almost completely coincided with that of an N sample. In the slow process (KWW2 process), there was a difference between the relaxation modulus of a P sample and that of an N sample. The results indicate that the KWW1 process in the empirical equation represents the relaxation in the collagen matrix in bone and that the KWW2 process is related to a higher-order structure of bone that is responsible for the anisotropic mechanical properties of bone.     

Nanostructural alteration in bone quantified in terms of orientation distribution of mineral crystals: a possible tool for fracture risk assessment

There may be different causes of failures in bone; however, their origin generally lies at the lowest level of structural hierarchy, i.e., at the mineral-collagen composite. Any change in the nanostructure affects the affinity or bonding effectiveness between and within the phases at this level, and hence determines the overall strength and quality of bone. In this study, we propose a novel concept to assess change in the nanostructure and thereby change in the bonding status at this level by revealing change in the orientation distribution characteristics of mineral crystals. Using X-ray diffraction method, a parameter called Degree of Orientation (DO) has been quantified. The DO accounts for the azimuthal distribution of mineral crystals and represents their effective amount along any direction. Changes in the DOs in cortical bone samples from bovine femur with different preferential orientations of mineral crystals were estimated under external loads. Depending on the applied loads, change in the azimuthal distribution of the DOs and the degree of reversibility of the crystals was observed to vary. The characteristics of nanostructural change and thereby possible affect on the strength of bone was then predicted from the reversible or irreversible characteristics of distributed mineral crystals. Significant changes in the organization of mineral crystals were observed; however, variations in the applied stresses and elastic moduli were not evinced at the macroscale level. A novel concept to assess the alteration in nanostructure on the basis of mineral crystals orientation distribution has been proposed. The importance of nanoscale level information obtained noninvasively has been emphasized, which acts as a precise tool to estimate the strength and predict the possible fracture risks in bone.