Resonant Ultrasound Spectroscopy (RUS) & Elasticity Lab

1. The Principle of Free Mechanical Vibrations

When a solid body with stress-free boundary conditions is mechanically excited, it vibrates at a discrete set of natural resonance frequencies ($f_n$). These frequencies depend strictly upon three physical properties:

• Sample Geometry and Dimensions: Precisely oriented rectangular parallelepipeds, cylinders, or spheres ($L_x, L_y, L_z$ measured to $\pm 1\,\mu\text{m}$).
• Mass Density ($\rho$): Measured via precision Archimedes hydrostatic weighing or pycnometry.
• The Complete Elastic Stiffness Tensor ($C_{ijkl}$): The 4th-rank tensor relating stress $\sigma_{ij}$ to strain $\epsilon_{kl}$ via Generalized Hooke's Law: $\sigma_{ij} = C_{ijkl} \epsilon_{kl}$.

Because different vibration modes (torsional, flexural, dilatational, and shear) depend differently upon different components of $C_{ijkl}$, acquiring a single frequency spectrum containing 30 to 60 distinct resonances provides sufficient mathematical constraints to determine all independent elastic constants simultaneously.

2. The Rayleigh-Ritz Lagrangian Method

The elastodynamic equation of motion for a freely vibrating solid is governed by the action principle $\delta \int L \, dt = 0$, where the Lagrangian is:

In the Rayleigh-Ritz formulation, the displacement components $u_i(\mathbf{x})$ are expanded in terms of complete Cartesian polynomial basis functions $\psi_\lambda = x^p y^q z^r$ (where $p + q + r \le N$, typically $N = 10\text{ to }14$, generating several hundred basis functions). Requiring $\partial L / \partial a_{i\lambda} = 0$ yields the matrix eigenvalue equation:

where $\mathbf{\Gamma}$ is the elastic stiffness matrix and $\mathbf{E}$ is the kinetic energy overlap matrix.

3. Voigt-Reuss-Hill Bounds & Zener Anisotropy

From the single-crystal stiffness tensor, polycrystalline average engineering moduli are computed using Voigt (uniform strain upper bound) and Reuss (uniform stress lower bound) formulations:

• Bulk Modulus: $K_V = K_R = (C_{11} + 2C_{12}) / 3$ (for cubic systems).
• Shear Modulus: $G_V = (C_{11} - C_{12} + 3C_{44}) / 5$ and $G_R = 5 / [4/(C_{11} - C_{12}) + 3/C_{44}]$. The Voigt-Reuss-Hill average is $G_{\text{VRH}} = (G_V + G_R) / 2$.
• Zener Anisotropy Ratio: $A = \frac{2C_{44}}{C_{11} - C_{12}}$. For isotropic materials, $A = 1.0$. Deviations from unity generate intense directional variations in sound velocity and acoustic attenuation along different crystal axes.

4. Quality Factor ($Q$) & Phase Transition Diagnostics

The resonance peak linewidth $\Delta f$ yields the mechanical quality factor $Q = f_0 / \Delta f$ and the internal friction $Q^{-1} = \tan \delta$. Measuring $f_n(T)$ and $Q^{-1}(T)$ as a function of temperature provides unmatched sensitivity to:

• Ferroelastic and Martensitic Transitions: Acoustic softening where specific shear moduli (e.g., $C' = (C_{11} - C_{12})/2$) drop toward zero.
• Superconductivity: Discontinuities in sound velocity at $T_c$ reflecting coupling between the superconducting order parameter and the crystal lattice (thermodynamic Ehrenfest relations).
• Dislocation Dynamics: High-temperature anelastic relaxation governed by Granato-Lücke string-model damping.