No abstract

With ageing and various diseases, the vascular pore volume fraction (porosity) in cortical bone increases, and the morphology of the pore network is altered. Cortical bone elasticity is known to decrease with increasing porosity, but the effect of the microstructure is largely unknown, while it has been thoroughly studied for trabecular bone. Also, popular micromechanical models have disregarded several micro-architectural features, idealizing pores as cylinders aligned with the axis of the diaphysis. The aim of this paper is to quantify the relative effects on cortical bone anisotropic elasticity of porosity and other descriptors of the pore network micro-architecture associated with pore number, size and shape. The five stiffness constants of bone assumed to be a transversely isotropic material were measured with resonant ultrasound spectroscopy in 55 specimens from the femoral diaphysis of 29 donors. The pore network, imaged with synchrotron radiation X-ray micro-computed tomography, was used to derive the pore descriptors and to build a homogenization model using the fast Fourier transform (FFT) method. The model was calibrated using experimental elasticity. A detailed analysis of the computed effective elasticity revealed in particular that porosity explains most of the variations of the five stiffness constants and that the effects of other micro-architectural features are small compared to usual experimental errors. We also have evidence that modelling the pore network as an ensemble of cylinders yields biased elasticity values compared to predictions based on the real micro-architecture. The FFT homogenization method is shown to be particularly efficient to model cortical bone.

Resonant ultrasound spectroscopy (RUS) is the state-of-the-art method used to investigate the elastic properties of anisotropic solids. Recently, RUS was applied to measure human cortical bone, an anisotropic material with low Q-factor (20), which is challenging due to the difficulty in retrieving resonant frequencies. Determining the precision of the estimated stiffness constants is not straightforward because RUS is an indirect method involving minimizing the distance between measured and calculated resonant frequencies using a model. This work was motivated by the need to quantify the errors on stiffness constants due to different error sources in RUS, including uncertainties on the resonant frequencies and specimen dimensions and imperfect rectangular parallelepiped (RP) specimen geometry. The errors were first investigated using Monte Carlo simulations with typical uncertainty values of experimentally measured resonant frequencies and dimensions assuming a perfect RP geometry. Second, the exact specimen geometry of a set of bone specimens were recorded by synchrotron radiation micro-computed tomography. Then, a "virtual" RUS experiment is proposed to quantify the errors induced by imperfect geometry. Results show that for a bone specimen of ∼1° perpendicularity and parallelism errors, an accuracy of a few percent ( <6.2%) for all the stiffness constants and engineering moduli is achievable.

The cortex of the femoral neck is a key structural element of the human body, yet there is not a reliable metric for predicting the mechanical properties of the bone in this critical region. This study explored the use of a range of non-destructive metrics to measure femoral neck cortical bone stiffness at the millimetre length scale. A range of testing methods and imaging techniques were assessed for their ability to measure or predict the mechanical properties of cortical bone samples obtained from the femoral neck of hip replacement patients. Techniques that can potentially be applied in vivo to measure bone stiffness, including computed tomography (CT), bulk wave ultrasound (BWUS) and indentation, were compared against in vitro techniques, including compression testing, density measurements and resonant ultrasound spectroscopy. Porosity, as measured by micro-CT, correlated with femoral neck cortical bone’s elastic modulus and ultimate compressive strength at the millimetre length scale. Large-tip spherical indentation also correlated with bone mechanical properties at this length scale but to a lesser extent. As the elastic mechanical properties of cortical bone correlated with porosity, we would recommend further development of technologies that can safely measure cortical porosity in vivo .

Bone extravascular matrix (EVM) elasticity at several tens micrometer scale plays a key role in the mechanical behavior of bone at different length scales with implications on bone biology through mechanotransduction.The elastic properties of cortical bone EVM have been evaluated by several experimental methods, including nanoindentation, scanning acoustic microscopy (SAM) and mechanical testing on µm sized bone specimens.Nevertheless, these methods hardly give access to elastic anisotropy. In this work, we propose a novel inverse homogenization method to evaluate the anisotropic elastic properties of cortical bone EVM based on the transverse isotropic elastic tensor of millimeter-sized bone specimens measured by using resonant ultrasound spectroscopy and Fast Fourier Transform homogenization method. With the inverse homogenization method, the anisotropic EVM stiffness constants were evaluated on 50 human femoral cortical bone specimens from an elderly group. To our knowledge, this is the first time that the whole set of the EVM stiffness tensor is evaluated on the same specimen and on a large number of samples. Further comparison with the results from SAM and the degree of mineralization of bone (DMB) showed the potential of this method. Empirical laws between DMB and EVM anisotropic stiffness constants were also provided for the first time. With the anisotropic elasticity evaluated by the proposed method, more accurate models can be developed to better understand bone mechanics and biology, such as mechanotransduction.

Elasticity assessment based on bulk wave velocity (BWV) measurements is the most popular technique to characterize the anisotropic stiffness tensor in cortical bone. Typically, a cuboid bone specimen is cut with its sides along the different anatomical directions. Then, the velocity of shear and longitudinal waves propagating along different directions are assessed, from which stiffness coefficients are calculated. Despite the importance of obtaining accurate elasticity values for bone research, there is no generally accepted protocol to measure BWV and the precision of the technique has been seldom investigated. The purpose of this work is to critically assess the method to measure BWV on cuboid specimens in terms of ultrasound frequency, specimen size and signal processing technique. In this study, we measured polycarbonate specimens of different dimensions and 55 human bone specimens with different transducers using frequencies ranging from 2.25 to 10 MHz and 1 to 5 MHz for longitudinal and shear waves, respectively. We compared four signal processing methods to detect the wave arrival time. The main results are that, (1) the measurement of shear waves is more complex than that of longitudinal wave, being less precise and more sensitive to sample size; (2) the estimated stiffness depends on the signal processing technique used (up to 10% variation for shear coefficients of bone); and (3) bone stiffness assessed from BWV using the first arrival of the signal to determine the time-of-flight is not different from stiffness assessed using resonant ultrasound spectroscopy (RUS). These results evidence that the measurement method can have an effect on the stiffness values estimates and hence, a well-defined protocol is needed to accurately measure bone stiffness coefficients based on BWV.

In the study of mechanical properties of human bone, specimens may be defatted before experiments to prevent contamination and the risk of infections. High energy synchrotron radiation micro-computed tomography (SR-μCT) is a popular technique to study bone microstructure. However, little is known about the effects of defatting or irradiation during SR-μCT imaging on different elastic coefficients including shear and longitudinal moduli in different anatomical directions. In this work, these effects are evaluated on a set of 24 samples using resonant ultrasound spectroscopy (RUS), which allows one to accurately measure the complete set of elastic coefficients of cortical bone non destructively. The results show that defatting with diethylether and methanol and irradiation up to 2.5kGy has no detectable effect on any of the elastic coefficients of human cortical bone.

Ultrasound imaging is a widely used diagnostic tool but has limitations in the imaging of deep lesions or obese patients where the large depth to aperture size ratio (f-number) reduces image quality. Reducing the f-number can improve image quality, and in this work, we combined three commercial arrays to create a large imaging aperture of 100 mm and 384 elements. To maintain the frame rate given the large number of elements, plane wave imaging was implemented with all three arrays transmitting a coherent wavefront. On wire targets at a depth of 100 mm, the lateral resolution is significantly improved; the lateral resolution was 1.27 mm with one array (1/3 of the aperture) and 0.37 mm with the full aperture. After creating virtual receiving elements to fill the inter-array gaps, an autoregressive filter reduced the grating lobes originating from the inter-array gaps by − 5.2 dB. On a calibrated commercial phantom, the extended field-of-view and improved spatial resolution were verified. The large aperture facilitates aberration correction using a singular value decomposition-based beamformer. Finally, after approval of the Stanford Institutional Review Board, the three-array configuration was applied in imaging the liver of a volunteer, validating the potential for enhanced resolution.