Quantification of stiffness measurement errors in resonant ultrasound spectroscopy of human cortical bone

Application of differential evolution on elasticity measurement of low quality factor materials using FEM-based resonant ultrasound spectroscopy

Finite element method based resonant ultrasound spectroscopy (FEM-based RUS) allows elasticity measurement for a material with high quality factor (Q) and arbitrary geometry by minimizing the differences between its theoretically calculated resonant frequencies and the corresponding experimentally measured ones. As Q decreases, some experimental frequencies remain undetermined, which makes it difficult to pair the calculated and experimental frequencies and to correctly identify the elastic constants. Additional difficulty need be tackled for irregularly-shaped low-Q materials due to the adoption of time-consuming FEM, thus efficiency of the identification method needs to be focused on. To apply FEM-based RUS to low-Q materials, a new elastic constant identification method is proposed based on a differential evolution algorithm in this paper. This method can perform a global search combining with local optimizations in the elastic constant space, and improve the overall efficiency by limiting the number of the frequency calculations. By using numerical experiments, the effectiveness of the proposed method under different frequency missing situations was verified and its efficiency was measured from the required frequency calculation numbers, showing an approximate two third reduction compared with an existing method. Finally, the elastic constants of an actual irregular cortical bone-mimicking material (Q ≈ 25) were measured using the two methods, yielding consistent Young's moduli (calculated from the identified constants) with the data provided by the manufacturer and a similar improvement in computational efficiency of the proposed method.

Maximum effect of the heterogeneity of tissue mineralization on the effective cortical bone elastic properties

The mineralization level is heterogeneous in cortical bone extracellular matrix as a consequence of remodeling. Models of the effective elastic properties at the millimeter scale have been developed based on idealizations of the vascular pore network and matrix properties. Some popular models do not take into account the heterogeneity of the matrix. However, the errors on the predicted elasticity when the difference in elastic properties between osteonal and interstitial tissues is not modeled have not been quantified. This work provides an estimation of the maximum error. We compare the effective elasticity of a representative volume element (RVE) assuming (1) different elastic properties in osteonal and interstitial tissues vs. (2) average matrix properties. In order to account for the variability of bone microstructure, we use a collection of high resolution images of the pore network to build RVEs. In each RVE we assumed a constant osteonal wall thickness and we artificially varied this thickness between 35 and 140 [Formula: see text]m to create RVEs with different amounts of osteonal tissue. The homogenization problem was solved with a fast Fourier transform (FFT)-based numerical scheme. We found that the error depends on pore volume fraction and varies on average from 1 to [Formula: see text] depending on the assumed diameter of the osteons. The results suggest that matrix heterogeneity may be disregarded in cortical bone models in most practical cases.

Bulk Wave Velocities in Cortical Bone Reflect Porosity and Compression Strength

The goal of this study was to evaluate whether ultrasonic velocities in cortical bone can be considered as a proxy for mechanical quality of cortical bone tissue reflected by porosity and compression strength. Micro-computed tomography, compression mechanical testing and resonant ultrasound spectroscopy were used to assess, respectively, porosity, strength and velocity of bulk waves of both shear and longitudinal polarisations propagating along and perpendicular to osteons, in 92 cortical bone specimens from tibia and femur of elderly human donors. All velocities were significantly associated with strength (r = 0.65-0.83) and porosity (r = -0.64 to -0.77). Roughly, according to linear regression models, a decrease in velocity of 100 m/s corresponded to a loss of 20 MPa in strength (which is approximately 10% of the largest strength value) and to an increase in porosity of 5%. These results provide a rationale for the in vivo measurement of one or several velocities for the diagnosis of bone fragility.

Frequency-dependent analysis of ultrasound apparent absorption coefficient in multiple scattering porous media: application to cortical bone

The effect of matrix viscoelastic absorption on frequency-dependent attenuation in porous structures mimicking simplified cortical bone is addressed in this numerical study. An apparent absorption is defined to quantify the difference between total attenuation (resulting from both absorption and scattering) and attenuation exclusively due to scattering. A power-law model is then used to describe the frequency-dependent apparent absorption as a function of pore diameter and density. The frequency response of the porous structures to a Gaussian pulse is studied to determine the frequency range over which the system can be considered linear. The results show that for low scattering regimes (normalized frequency [Formula: see text]0.80), the system and its apparent absorption can be considered linear. Hence, the total attenuation coefficient results from the summation of scattering and absorption coefficients. However, for highly scattering regimes, the system can no longer be considered linear, as the apparent absorption vs. frequency deviates from a linear trend. As the pore density increases, the apparent absorption coefficient increases as well.

Homogenization of cortical bone reveals that the organization and shape of pores marginally affect elasticity

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.

Artificial neural network to estimate micro-architectural properties of cortical bone using ultrasonic attenuation: A 2-D numerical study

The goal of this study is to estimate micro-architectural parameters of cortical porosity such as pore diameter (φ), pore density (ρ) and porosity (ν) of cortical bone from ultrasound frequency dependent attenuation using an artificial neural network (ANN). First, heterogeneous structures with controlled pore diameters and pore densities (mono-disperse) were generated, to mimic simplified structure of cortical bone. Then, more realistic structures were obtained from high resolution CT scans of human cortical bone. 2-D finite-difference time-domain simulations were conducted to calculate the frequency-dependent attenuation in the 1-8 MHz range. An ANN was then trained with the ultrasonic attenuation at different frequencies as the input feature vectors while the output was set as the micro-architectural parameters (pore diameter, pore density and porosity). The ANN is composed of three fully connected dense layers with 24, 12 and 6 neurons, connected to the output layer. The dataset was trained over 6000 epochs with a batch size of 16. The trained ANN exhibits the ability to predict the micro-architectural parameters with high accuracy and low losses. ANN approaches could potentially be used as a tool to help inform physics-based modelling of ultrasound propagation in complex media such as cortical bone. This will lead to the solution of inverse-problems to retrieve bone micro-architectural parameters from ultrasound measurements for the non-invasive diagnosis and monitoring osteoporosis.

Anisotropic elastic properties of human femoral cortical bone and relationships with composition and microstructure in elderly

The strong dependence between cortical bone elasticity at the millimetre-scale (mesoscale) and cortical porosity has been evidenced by previous studies. However, bone is an anisotropic composite material made by mineral, proteins and water assembled in a hierarchical structure. Whether the variations of structural and compositional properties of bone affect the different elastic coefficients at the mesoscale is not clear. Aiming to understand the relationships between bone elastic properties and compositions and microstructure, we applied state-of-the-art experimental modalities to assess these aspects of bone characteristics. All elastic coefficients (stiffness tensor of the transverse isotropic bone material), structure of the vascular pore network, collagen and mineral properties were measured in 52 specimens from the femoral diaphysis of 26 elderly donors. Statistical analyses and micromechanical modeling showed that vascular pore volume fraction and the degree of mineralization of bone are the most important determinants of cortical bone anisotropic mesoscopic elasticity. Though significant correlations were observed between collagen properties and elasticity, their effects in bone mesoscopic elasticity were minor in our data. This work also provides a unique set of data exhibiting a range of variations of compositional and microstructural cortical bone properties in the elderly and gives strong experimental evidence and basis for further development of biomechanical models for human cortical bone. STATEMENT OF SIGNIFICANCE: This study reports the relationships between microstructure, composition and the mesoscale anisotropic elastic properties of human femoral cortical bone in elderly. For the first time, we provide data covering the complete anisotropic elastic tensor, the microstructure of cortical vascular porosity, mineral and collagen characteristics obtained from the same or adjacent samples in each donor. The results revealed that cortical vascular porosity and degree of mineralization of bone are the most important determinants of bone anisotropic stiffness at the mesoscale. The presented data gives strong experimental evidence and basis for further development of biomechanical models for human cortical bone.

Assessment of elastic coefficients of child cortical bone using resonant ultrasound spectroscopy

The assessment of the anisotropic elastic properties of non-pathological child cortical bone remains a challenge for the biomechanical engineering community and an important clinical issue. Resonant ultrasound spectroscopy (RUS) can be used to determine bone stiffness coefficients from the mechanical resonances of bone specimens. Here, a RUS protocol was used on 7 fibula specimens from children (mean age 14 ± 3 years) to estimate the whole elastic stiffness tensor of non-pathological child cortical bone considered as orthotropic. Despite a small number of sample, results are consistent with this hypothesis, even if a trend towards transverse isotropy is discussed. Indeed, the average values of the 9 independent stiffness coefficients obtained in this study for child bone are: C₁₁ = 16.73 ± 0.19 GPa, C₂₂ = 16.19 ± 0.12 GPa, C₃₃ = 24.47 ± 0.30 GPa, C₄₄ = 4.14 ± 0.08 GPa, C₅₅ = 4.16 ± 0.07 GPa, C₆₆ = 3.13 ± 0.05 GPa, C₁₂ = 10.14 ± 0.20 GPa, C₁₃ = 10.67 ± 0.27 GPa, C₂₃ = 10.25 ± 0.14 GPa.

Elastic properties measurement of human enamel based on resonant ultrasound spectroscopy

OBJECTIVES

To investigate the elastic properties of human enamel using resonant ultrasound spectroscopy (RUS).

METHODS

Six rectangular parallelepiped specimens were prepared from six human third molars. For all specimens, the theoretical resonant frequencies were calculated using the Rayleigh-Ritz method, knowing the specimen mass density and dimensions, and using a priori stiffness constants. The experimental resonant frequencies were measured and extracted by RUS. Then, the optimal stiffness constants were retrieved by adjustment of the theoretical resonant frequencies to the measured ones based on the Levenberg-Marquardt method. The engineering elastic moduli, including Young's moduli, shear moduli, and Poisson's ratios, were also calculated based on the optimal stiffness constants.

RESULTS

The five independent stiffness constants C₁₁, C₁₂, C₁₃, C₃₃, and C₄₄ were 90.2 ± 6.65 GPa, 34.7 ± 6.90 GPa, 29.5 ± 4.82 GPa, 83.5 ± 8.93 GPa, and 37.0 ± 10.9 GPa, respectively. Young's moduli E₁₁ and E₃₃, shear moduli G₁₃ and G₁₂, and Poisson's ratios υ₁₂ and υ₁₃ were 71.7 ± 7.34 GPa, 69.2 ± 7.32 GPa, 37.0 ± 10.9 GPa, 28.1 ± 4.35 GPa, 0.303 ± 0.098, and 0.248 ± 0.060, respectively.

SIGNIFICANCE

Elastic properties are critical for developing dental materials and designing dental prostheses. The RUS method may provide more precise measurement of elastic properties of dental materials.