phasevelocity.m - A Matlab script to calculate the frequency-dependent phase velocity and radial variation of elastic waves in cylindrical bars

Description

The relationships between normalised frequency and phase velocity, M₁ and M₂ can be used to account for first-mode dispersion effects in pressure bar measurements using programs such as dispersion.m. These relationships are dependent on the Poisson's ratio of the bar material, and so can be generated for any Poisson's ratio using this script, and saved in the file format used by dispersion.m.

The phasevelocity.m script finds the first root of Bancroft's (1941) equation using the bisection method, for a defined Poisson's ratio, and over a defined range of normalised wavelength (d/L). The result is the normalised phase velocity, c_p/c₀, which corresponds to the first mode of propagation for longitudinal waves in an elastic cylindrical bar. Normalised wavelengths are also converted to normalised frequency, fa/c₀.

Normalised phase velocities are then used to calucate Tyas and Wilson's (2001) factors M₁ and M₂, which account for wavelength-dependent radial variations in strain and Young's modulus, respectively. All values are recorded in .mat format.

References

• Bancroft, D. (1941) The Velocity of Longitudinal Waves in Cylindrical Bars. Physical Review, 59, 588-593.
• Barr, A. D. (2016) dispersion.m - A MatLab script for phase angle and amplitude correction of pressure bar signals. The University of Sheffield. Software. https://doi.org/10.15131/shef.data.3996876.v1 
• Tyas, A., Wilson, A. J. (2001) An investigation of frequency domain dispersion correction of pressure bar signals. International Journal of Impact Engineering, 25, 87-101.