Search

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 [+ or -] 1.77 MPa, while the corresponding mean strain was 7.15 [+ or -] 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[+ or -]7.91 MPa, and the dissipated energy density was 5.67[+ or -]4.42 MPa. The latter variable was strongly related to the volume fraction (R2 =0.83) but the correlation improved only marginally with the inclusion of fabric (R2 =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. [DOI: 10.1115/1.4001361] Keywords: densification, fabric, large strain, mechanical properties, morphology, softening, trabecular bone

Classical finite element (FE) models can estimate vertebral stiffness and strength with much lower computational costs than [micro]FE analyses, but the accuracy of these models rely on calibrated material properties that are not necessarily consistent with experimental results, In general, trabecular bone material properties are scaled with computer tomography (CT) density alone, without accounting for local variations in anisotropy or micro-architecture. Moreover, the cortex is often omitted or assigned with a constant thickness. In this work, voxel FE models, as well as surface-based homogenized FE models with topologically-conformed geometry and assigned with experimentally validated properties for bone, were developed from a series of 12 specimens tested up to failure. The effects of changing from a digital mesh to a smooth mesh, including a cortex of variable thickness and/or including heterogeneous trabecular fabric, were investigated. In each case, FE predictions of vertebral stiffness and strength were compared with the experimental gold-standard, and changes in elastic strain energy density and damage distributions were reported. The results showed that a smooth mesh effectively removed zones of artificial damage locations occurring in the ragged edges of the digital mesh. Adding an explicit cortex stiffened and strengthened the models, unloading the trabecular centrum while increasing the correlations to experimental stiffness and strength. Further addition of heterogeneous fabric improved the correlations to stiffness ([R.sup.2] = 0.72) and strength ([R.sup.2] = 0.89) and moved the damage locations closer to the vertebral endplates, following the local trabecular orientations. B was furthermore demonstrated that predictions of vertebral stiffness and strength of homogenized FE models with topologically-conformed cortical shell and heterogeneous trabecular fabric correlated well with experimental measurements, after assigning purely experimental data for bone without further calibration of material laws at the macroscale of bone. This study successfully demonstrated the limitations of current classical FE methods and provided valuable insights into the damage mechanisms of vertebral bodies. [DOI: 10.1115/1.3212097]

The identification of anisotropic elastic properties of lamellar bone based on nanoindentation data is an open problem. Therefore, the purpose of this study was to develop a method to estimate the orthotropic elastic constants of human cortical bone secondary osteons using nanoindentation in two orthogonal directions. Since the indentation modulus depends on all elastic constants and, for anisotropic materials, also on the indentation direction, a theoretical model quantifying the indentation modulus from the stiffness tensor of a given material was implemented numerically (Swadener and Pharr, 2001, 'Indentation of Elastically Anisotropic Half-Spaces by Cones and Parabolae of Revolution,' Philos. Mag. A, 81(2), pp. 447-466). Nanoindentation was performed on 22 osteons of the distal femoral shaft: A new holding system was designed in order to indent the same osteon in two orthogonal directions. To interpret the experimental results and identify orthotropic elastic constants, an inverse procedure was developed by using a fabric-based elastic model for lamellar bone. The experimental indentation moduli were found to vary with the indentation direction and showed a marked anisotropy. The estimated elastic constants showed different degrees of anisotropy among secondary osteons of the same bone and these degrees of anisotropy were also found to be different than the one of cortical bone at the macroscopic level. Using the log-Euclidean norm, the relative distance between the compliance tensors of the estimated mean osteon and of cortical bone at the macroscopic level was 9.69%: Secondary osteons appeared stiffer in their axial and circumferential material directions, and with a greater bulk modulus than cortical bone, which is attributed to the absence of vascular porosity in osteonal properties. The proposed method is suitable for identification of elastic constants from nanoindentation experiments and could be adapted to other (bio)materials, for which it is possible to describe elastic properties using a fabric-based model. Keywords: anisotropy, cortical bone, elasticity, nanoindentation, secondary osteons

Measuring the microscopic mechanical properties of bone tissue is important in support of understanding the etiology and pathogenesis of many bone diseases. Knowledge about these properties provides a context for estimating the local mechanical environment of bone related cells that coordinate the adaptation to loads experienced at the whole organ level. The objective of this study was to determine the effects of experimental testing parameters on nanoindentation measures of lamellar-level bone mechanical properties. Specifically, we examined the effect of specimen preparation condition, indentation depth, repetitive loading, time delay, and displacement rate. The nanoindentation experiments produced measures of lamellar elastic moduli for human cortical bone (average value of 17.7 [+ or -] 4.0 GPa for osteons and 19.3 [+ or -] 4.7 GPa for interstitial bone tissue). In addition, the hardness measurements produced results consistent with data in the literature (average 0.52 [+ or -] 0.15 GPa for osteons and 0.59 [+ or -] 0.20 GPa for interstitial bone tissue). Consistent modulus values can be obtained from a 500-nm-deep indent. The results also indicated that the moduli and hardnesses of the dry specimens are significantly greater (22.6% and 56.9%, respectively) than those of the wet and wet and embedded specimens. The latter two groups were not different. The moduli obtained at a 5-nm/s loading rate were significantly lower than the values at the 10- and 20-nm/s loading rates while the 10- and 20-nm/s rates were not significantly different. The hardness measurements showed similar rate-dependent results. The preliminary results indicated that interstitial bone tissue has significantly higher modulus and hardness than osteonal bone tissue. In addition, a significant correlation between hardness and elastic modulus was observed.

In most finite element (FE) studies of vertebral bodies, axial compression is the loading mode of choice to investigate structural properties, but this might not adequately reflect the various loads to which the spine is subjected during daily activities or the increased fracture risk associated with shearing or bending loads. This work aims at proposing a patient-specific computer tomography (CT)-based methodology, using the currently most advanced, clinically applicable finite element approach to perform a structural investigation of the vertebral body by calculation of its full six dimensional (6D) stiffness matrix. FE models were created from voxel images after smoothing of the peripheral voxels and extrusion of a cortical shell, with material laws describing heterogeneous, anisotropic elasticity for trabecular bone, isotropic elasticity for the cortex based on experimental data. Validated against experimental axial stiffness, these models were loaded in the six canonical modes and their 6D stiffness matrix calculated. Results show that, on average, the major vertebral rigidities correlated well or excellently with the axial rigidity but that weaker correlations were observed for the minor coupling rigidities and for the image-based density measurements. This suggests that axial rigidity is representative of the overall stiffness of the vertebral body and that finite element analysis brings more insight in vertebral fragility than densitometric approaches. Finally, this extended patient-specific FE methodology provides a more complete quantification of structural properties for clinical studies at the spine. [ABSTRACT FROM AUTHOR]