Measurement of Shear Elastic Moduli in Quasi-Incompressible Soft Solids

Abstract

Recently a nonlinear equation describing the plane shear wave propagation in isotropic quasi‐incompressible media has been developed using a new expression of the strain energy density, as a function of the second, third and fourth order shear elastic constants (respectively μ, A, D) [1]. In such a case, the shear nonlinearity parameter βs depends only from these last coefficients. To date, no measurement of the parameter D have been carried out in soft solids. Using a set of two experiments, acoustoelasticity and finite amplitude shear waves, the shear elastic moduli up to the fourth order of soft solids are measured. Firstly, this theoretical background is applied to the acoustoelasticity theory, giving the variations of the shear wave speed as a function of the stress applied to the medium. From such variations, both linear (μ) and third order shear modulus (A) are deduced in agar‐gelatin phantoms. Experimentally the radiation force induced by a focused ultrasound beam is used to generate quasi‐plane linear shear waves within the medium. Then the shear wave propagation is imaged with an ultrafast ultrasound scanner. Secondly, in order to give rise to finite amplitude plane shear waves, the radiation force generation technique is replaced by a vibrating plate applied at the surface of the phantoms. The propagation is also imaged using the same ultrafast scanner. From the assessment of the third harmonic amplitude, the nonlinearity parameter βS is deduced. Finally, combining these results with the acoustoelasticity experiment, the fourth order modulus (D) is deduced. This set of experiments provides the characterization, up to the fourth order, of the nonlinear shear elastic moduli in quasi‐incompressible soft media. Measurements of the A moduli reveal that while the behaviors of both soft solids are close from a linear point of view, the corresponding nonlinear moduli A are quite different. In a 5% agar‐gelatin phantom, the fourth order elastic constant D is found to be 30±10 kPa.