Anatomical variation of orthotropic elastic moduli of the proximal human tibia

The anatomical variation of orthotropic elastic moduli of the cancellous bone from three human proximal tibiae was investigated using an ultrasonic technique. With this technique, it was possible to measure three orthogonal elastic moduli and three shear moduli from cubic specimens of cancellous bone as small as 8 mm per side. Correlation with mechanical tensile testing has shown this technique to offer a precise measure of cancellous modulus (Eten = 0.94Eult + 144.6 MPa, r2 = 0.96, n = 34). The cancellous bone of the proximal tibia was found to be very inhomogeneous, with the axial modulus ranging between 340 and 3350 MPa. A course map is presented, showing measured Young's moduli as a function of anatomical position. The anisotropy of the cancellous bone, determined by the relative differences between the three orthogonal moduli, was shown to be relatively constant over the entire range of cancellous densities tested. The relationship between the axial elastic modulus and the apparent density was found to be approximately linear, as reported by others for proximal tibial cancellous bone.     

Validation of a voxel-based FE method for prediction of the uniaxial apparent modulus of human trabecular bone using macroscopic mechanical tests and nanoindentation

Assessment of the mechanical properties of trabecular bone is of major biological and clinical importance for the investigation of bone diseases, fractures and their treatments. Finite element (FE) methods are getting increasingly popular for quantifying the elastic and failure properties of trabecular bone. In particular, voxel-based FE methods have been previously used to calculate the effective elastic properties of trabecular microstructures. However, in most studies, bone tissue moduli were assumed or back-calculated to match the apparent elastic moduli from experiments, which often lead to surprisingly low values when compared to nanoindentation results. In this study, voxel-based FE analysis of trabecular bone is combined with physical measures of volume fraction, micro-CT (microCT) reconstructions, uniaxial mechanical tests and specimen-specific nanoindentation tests for proper validation of the method. Cylindrical specimens of cancellous bone were extracted from human femurs and their volume fraction determined with Archimede's method. Uniaxial apparent modulus of the specimens was measured with an improved tension-compression testing protocol that minimizes boundary artefacts. Their microCT reconstructions were segmented to match the measured bone volume fraction and used to create full-size voxel models with 30-45 microm element size. For each specimen, linear isotropic elastic material properties were defined based on specific nanoindentation measurements of its embedded bone tissue. Linear FE analyses were finally performed to simulate the uniaxial mechanical tests. Additional parametric analyses were performed to evaluate the potential errors on the predicted apparent modulus arising from variations in segmentation threshold, tissue modulus, and the use of 125-mm(3) cubic sub-regions. The results demonstrate an excellent correspondence between experimental measures and FE predictions of uniaxial apparent modulus. In conclusion, the adopted voxel-based FE approach is found to be a robust method to predict the linear elastic properties of human cancellous bone, provided segmentation of the microCT reconstructions is carefully calibrated, tissue modulus is known a priori and the entire region of interest is included in the analysis.     

Tensile and compressive stress yield criteria for cancellous bone

Experimentally based isotropic tensile and compressive strain yield criteria for cancellous bone are shown to imply tensile and compressive stress yield criteria, respectively. In particular, it is deduced that the directional variation of uniaxial yield stress is proportional to the one of Young's modulus. The obtained tensile and compressive stress yield criteria, in turn, imply information about the total stress yield criterion for cancellous bone.     

Mechanical properties of femoral trabecular bone in dogs

Background  

Studying mechanical properties of canine trabecular bone is important for a better understanding of fracture mechanics or bone disorders and is also needed for numerical simulation of canine femora. No detailed data about elastic moduli and degrees of anisotropy of canine femoral trabecular bone has been published so far, hence the purpose of this study was to measure the elastic modulus of trabecular bone in canine femoral heads by ultrasound testing and to assess whether assuming isotropy of the cancellous bone in femoral heads in dogs is a valid simplification.     

Finite element models predict cancellous apparent modulus when tissue modulus is scaled from specimen CT-attenuation

High-resolution architecture-based finite element models are commonly used for characterizing the mechanical behavior of cancellous bone. The vast majority of studies use homogeneous material properties to model trabecular tissue. The objectives of this study were to demonstrate that inhomogeneous finite element models that account for microcomputed tomography-measured tissue modulus variability more accurately predict the apparent stiffness of cancellous bone than homogeneous models, and to examine the sensitivity of an inhomogeneous model to the degree of tissue property variability. We tested five different material cases in finite element models of ten cancellous cubes in simulated uniaxial compression. Three of these cases were inhomogeneous and two were homogeneous. Four of these cases were unique to each specimen, and the remaining case had the same tissue modulus for all specimens. Results from all simulations were compared with measured elastic moduli from previous experiments. Tissue modulus variability for the most accurate of the three inhomogeneous models was then artificially increased to simulate the effects of non-linear CT-attenuation-modulus relationships. Uniqueness of individual models was more critical for model accuracy than level of inhomogeneity. Both homogeneous and inhomogeneous models that were unique to each specimen had at least 8% greater explanatory power for apparent modulus than models that applied the same material properties to all specimens. The explanatory power for apparent modulus of models with a tissue modulus coefficient of variation (COV) range of 21-31% was 13% greater than homogeneous models (COV=0). The results of this study indicate that inhomogenous finite element models that have tissue moduli unique to each specimen more accurately predict the elastic behavior of cancellous cubic specimens than models that have common tissue moduli between all specimens.     

The limitations of canine trabecular bone as a model for human: a biomechanical study

Distal canine femurs were sectioned into 8 mm cubic specimens. Orthogonal compression tests were performed to preyield in two or three directions and to failure in a third. Apparent density and ash weight density were measured for a subset of specimens. The results were compared to the human distal femur results of Ciarelli et al. (Transactions of the 32nd Annual Meeting of the Orthopaedic Research Society, Vol. 11, p. 42, 1986). Quantitative similarities existed in the fraction of components comprising the trabecular tissue of the two species. Qualitative similarities were seen in the positional and anisotropic variation of the mechanical properties, and also in the form and strength of the relationships between the mean modulus and bone density, ultimate stress and density, and ultimate stress and modulus. However, significantly different regression equations resulted for the mean modulus-density, and ultimate stress modulus relationships, indicating that for the same density, canine trabecular bone displays a lower modulus than human, and may achieve greater compressive strains before failure.     

Linear poroelastic cancellous bone anisotropy: trabecular solid elastic and fluid transport properties

The mechanical performance of cancellous bone is characterized using experiments which apply linear poroelasticity theory. It is hypothesized that the anisotropic organization of the solid and pore volumes of cancellous bone can be physically characterized separately (no deformable boundary interactive effects) within the same bone sample. Due to its spongy construction, the in vivo mechanical function of cancellous or trabecular bone is dependent upon fluid and solid materials which may interact in a hydraulic, convective fashion during functional loading. This project provides insight into the organization of the tissue, ie., the trabecular connectivity, by defining the separate nature of this biphasic performance. Previous fluid flow experiments [Kohles et al., 2001, Journal of Biomechanics, 34(11), pp. 1197-1202] describe the pore space via orthotropic permeability. Ultrasonic wave propagation through the trabecular network is used to describe the solid component via orthotropic elastic moduli and material stiffness coefficients. The linear poroelastic nature of the tissue is further described by relating transport (fluid flow) and elasticity (trabecular load transmission) during regression analysis. In addition, an empirical relationship between permeability and porosity is applied to the collected data. Mean parameters in the superior-inferior (SI) orientation of cubic samples (n=20) harvested from a single bovine distal femur were the largest (p<0.05) in comparison to medial-lateral (ML) and anterior-posterior (AP) orientations: Apparent elastic modulus (2,139 MPa), permeability (4.65x10(-10) m2), and material stiffness coefficient (13.6 GPa). A negative correlation between permeability as a predictor of structural elastic modulus supported a parametric relationship in the ML (R2=0.4793), AP (R2=0.3018), and SI (R2=0.6445) directions (p<0.05).     

Elastic modulus of trabecular bone material

An ultrasonic technique was used to measure both the elastic modulus (Young's modulus) of trabecular bone material and the elastic modulus of the cancellous structure. The average trabecular modulus, measured on specimens obtained from three human and one bovine distal femora, was 13.0 GPa (S.D. 1.47) and 10.9 GPa (S.D. 1.57), respectively. On human specimens the structural elastic modulus was found to be related to the structural (apparent) density raised to the 1.88 power. The elastic modulus from the bovine specimens showed a more linear relationship with the density of the cancellous structure (density raised to the 1.57 power).     

Tensile behavior of cortical bone: dependence of organic matrix material properties on bone mineral content

A porous composite model is developed to analyze the tensile mechanical properties of cortical bone. The effects of microporosity (volksman's canals, osteocyte lacunae) on the mechanical properties of bone tissue are taken into account. A simple shear lag theory, wherein tensile loads are transferred between overlapped mineral platelets by shearing of the organic matrix, is used to model the reinforcement provided by mineral platelets. It is assumed that the organic matrix is elastic in tension and elastic-perfectly plastic in shear until it fails. When organic matrix shear stresses at the ends of mineral platelets reach their yield values, the stress-strain curve of bone tissue starts to deviate from linear behavior. This is referred as the microscopic yield point. At the point where the stress-strain behavior of bone shows a sharp curvature, the organic phase reaches its shear yield stress value over the entire platelet. This is referred as the macroscopic yield point. It is assumed that after macroscopic yield, mineral platelets cannot contribute to the load bearing capacity of bone and that the mechanical behavior of cortical bone tissue is determined by the organic phase only. Bone fails when the principal stress of the organic matrix is reached. By assuming that mechanical properties of the organic matrix are dependent on bone mineral content below the macroscopic yield point, the model is used to predict the entire tensile mechanical behavior of cortical bone for different mineral contents. It is found that decreased shear yield stresses and organic matrix elastic moduli are required to explain the mechanical behavior of bones with lowered mineral contents. Under these conditions, the predicted values (elastic modulus, 0.002 yield stress and strain, and ultimate stress and strain) are within 15% of experimental data.     

Comparison of the elastic and yield properties of human femoral trabecular and cortical bone tissue

The ability to determine trabecular bone tissue elastic and failure properties has biological and clinical importance. To date, trabecular tissue yield strains remain unknown due to experimental difficulties, and elastic moduli studies have reported controversial results. We hypothesized that the elastic and tensile and compressive yield properties of trabecular tissue are similar to those of cortical tissue. Effective tissue modulus and yield strains were calibrated for cadaveric human femoral neck specimens taken from 11 donors, using a combination of apparent-level mechanical testing and specimen-specific, high-resolution, nonlinear finite element modeling. The trabecular tissue properties were then compared to measured elastic modulus and tensile yield strain of human femoral diaphyseal cortical bone specimens obtained from a similar cohort of 34 donors. Cortical tissue properties were obtained by statistically eliminating the effects of vascular porosity. Results indicated that mean elastic modulus was 10% lower (p<0.05) for the trabecular tissue (18.0+/-2.8 GPa) than for the cortical tissue (19.9+/-1.8 GPa), and the 0.2% offset tensile yield strain was 15% lower for the trabecular tissue (0.62+/-0.04% vs. 0.73+/-0.05%, p<0.001). The tensile-compressive yield strength asymmetry for the trabecular tissue, 0.62 on average, was similar to values reported in the literature for cortical bone. We conclude that while the elastic modulus and yield strains for trabecular tissue are just slightly lower than those of cortical tissue, because of the cumulative effect of these differences, tissue strength is about 25% greater for cortical bone.     

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

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

The bone tissue of the canine mandible is elastically isotropic

This paper reports experimental measurements which show that canine mandibular bone tissue is elastically isotropic. Earlier work has established that human, canine and bovine cortical bone tissue of the femur, tibia and skull are elastically anisotropic and therefore the reported isotropy of mandibular tissue was unexpected. The isotropic elastic moduli of the canine mandible are represented by a Young's modulus of 7.5 GPa and a Poisson's ratio of 0.4. Earlier work gave the three orthotropic Young's moduli of the cortical one of the canine femur as 12.8 GPa, 15.6 GPa and 20.1 GPa. The experimental technique employed is elastic wave propagation at ultrasonic frequencies.     

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.     

Alterations in damage processes in dense cancellous bone following gamma-radiation sterilization

Structurally intact cancellous bone allograft is an attractive tissue form because its high porosity can provide space for delivery of osteogenic factors and also allows for more rapid and complete in-growth of host tissues. Gamma radiation sterilization is commonly used in cancellous bone allograft to prevent disease transmission. Commonly used doses of gamma radiation sterilization (25–35 kGy) have been shown to modify cortical bone post-yield properties and crack propagation but have not been associated with changes in cancellous bone material properties. The purpose of this study was to determine the effects of irradiation on the elastic and yield properties and microscopic tissue damage processes in dense cancellous bone. Cancellous bone specimens (13 control, 14 irradiated to 30 kGy) from bovine proximal tibiae were tested in compression to 1.3% apparent strain and examined for microscopic tissue damage. The yield strain in irradiated specimens (0.93±0.11%, mean±SD) did not differ from that in control specimens (0.90±0.11%, p=0.44). No differences in elastic modulus were observed between groups after accounting for differences in bone volume fraction. However, irradiated specimens showed greater residual strain (p=0.01), increased number of microfractures (p=0.02), and reduced amounts of cross-hatching type damage (p<0.01). Although gamma radiation sterilization at commonly used dosing (30 kGy) does not modify elastic or yield properties of dense cancellous bone, it does cause modifications in damage processes, resulting in increased permanent deformation following isolated overloading.     

On the dependence of the elasticity and strength of cancellous bone on apparent density

This paper presents a statistical analysis of the pooled data from a number of previous experiments concerning the dependence of the Young's moduli and strength of cancellous bone tissue upon apparent density. The results show that both the Young's moduli and the strength are proportional to the square of apparent density of the tissue and are therefore proportional to one another. It is shown that the coefficient of proportionality is different for human and bovine tissue. It is concluded that the suggestion of Wolff (Das Gesetz der Transformation der Knochen, Hirschwald, Berlin, 1892) that compact bone tissue is simply more dense cancellous bone tissue is not an accurate statement when only the mechanical properties of these two tissues are considered. It is noted that estimates for the elastic modulus of the individual trabecula of human cancellous bone vary from 1 to 20 GPa and it is suggested that this question needs further study.     

The fabric dependence of the orthotropic elastic constants of cancellous bone

It has been proposed that the orthotropic elastic constants of cancellous bone depend upon a tensorial measure of anisotropy called fabric as well as the tissue's structural density. Cowin (1985, Mechanics Mater, 4, 137-147; 1986, J. biomech. Engng 108, 83-88) developed explicit relationships for the elastic constant, structural density and fabric relationship. In this study the orthotropic elastic moduli, structural density, and fabric components were measured for 11 cancellous bone specimens from five bovine femora and for 75 specimens from three human proximal tibiae and fitted to these relationships using a least squares analysis. The relationships explained between 72 and 94% of the variance in the elastic constants. The relationships between the elastic constants and squared or cubed power functions of structural density had better predictive value over the entire distribution of the data than did expressions with linear functions of structural density.     

Tensile damage and its effects on cortical bone

Plexiform bovine bone samples are repeatedly loaded in tension along their longitudinal axis. In order to induce damage in the bone tissue, bone samples are loaded past their yield point. Half of the bone samples from the damaged group were stored in saline to allow for viscoelastic recovery while the others were decalcified. Tensile tests were conducted on these samples to characterize the effects of damage on the mechanical behavior of the organic matrix (decalcified samples) as well as on bone tissue (stored in saline). The ultimate strain of the damaged decalcified bone is 29% higher compared to that of non-damaged decalcified (control) bone. The ultimate stresses as well as the elastic moduli are similar in both decalcified groups. This phenomenon is also observed in other collagenous tissue (tendon and ligament). This may suggest that damage in bone is caused by shear failure of the organic matrix; transverse separation of the collagen molecules or microfibrils from each other. In contrast, there is a trend towards lowered ultimate strains in damaged bone, which is soaked in saline, with respect to control bone samples (not damaged). The damaged bone tissue exhibits a bi-linear behavior in contrast to the mechanical behavior of non-damaged bone. The initial elastic modulus (below 55 MPa) and ultimate strength of damaged bone are similar to that in non-damaged bone.     

Changing the structurally effective mineral content of bone with in vitro fluoride treatment

Bovine femur cortical bone specimens were tested in tension after being treated in vitro for 3 days with sodium fluoride solutions of different molarity (0.145, 0.5, and 2.0M). The treatments alter the mechanical properties of the bone samples with different degrees as compared to control samples (untreated). The mechanical properties of the treated samples have lower elastic modulus, yield and ultimate stress, acoustic impedance and hardness, and higher ultimate strain and toughness as compared to control samples. The observed effects were intensified with the increasing molarity of the treatment solutions. This study shows that the fluoride treatment can be used to investigate the composite behavior of bone tissue by altering the structurally important bone mineral content in a controlled manner.     

The role of fabric in the large strain compressive behavior of human trabecular bone

Osteoporosis-related vertebral body fractures involve large compressive strains of trabecular bone. The small strain mechanical properties of the trabecular bone such as the elastic modulus or ultimate strength can be estimated using the volume fraction and a second order fabric tensor, but it remains unclear if similar estimations may be extended to large strain properties. Accordingly, the aim of this work is to identify the role of volume fraction and especially fabric in the large strain compressive behavior of human trabecular bone from various anatomical locations. Trabecular bone biopsies were extracted from human T12 vertebrae (n=31), distal radii (n=43), femoral head (n=44), and calcanei (n=30), scanned using microcomputed tomography to quantify bone volume fraction (BV/TV) and the fabric tensor (M), and tested either in unconfined or confined compression up to very large strains (～70%). The mechanical parameters of the resulting stress-strain curves were analyzed using regression models to examine the respective influence of BV/TV and fabric eigenvalues. The compressive stress-strain curves demonstrated linear elasticity, yielding with hardening up to an ultimate stress, softening toward a minimum stress, and a steady rehardening followed by a rapid densification. For the pooled experiments, the average minimum stress was 1.89 ± 1.77 MPa, while the corresponding mean strain was 7.15 ± 1.84%. The minimum stress showed a weaker dependence with fabric as the elastic modulus or ultimate strength. For the confined experiments, the stress at a logarithmic strain of 1.2 was 8.08 ± 7.91 MPa, and the dissipated energy density was 5.67 ± 4.42 MPa. The latter variable was strongly related to the volume fraction (R(2)=0.83) but the correlation improved only marginally with the inclusion of fabric (R(2)=0.84). The influence of fabric on the mechanical properties of human trabecular bone decreases with increasing strain, while the role of volume fraction remains important. In particular, the ratio of the minimum versus the maximum stress, i.e., the relative amount of softening, decreases strongly with fabric, while the dissipated energy density is dominated by the volume fraction. The collected results will prove to be useful for modeling the softening and densification of the trabecular bone using the finite element method.