Analysis of Floquet Waves in Periodic Multilayered Isotropic Media with the Method of Reverberation-Ray Matrix

Abstract

The in-plane elastic waves in periodically multilayered isotropic structures, which are decoupled from the out-of-plane waves, are represented mainly by the frequency–wavenumber spectra and occasionally by the frequency–phase velocity spectra as well as being studied predominantly for periodic bi-layered media along and perpendicular to the thickness direction in the existing research. This paper investigates their comprehensive dispersion characteristics along arbitrary in-plane directions and in entire (low and high) frequency ranges, including the frequency–wavelength, wavenumber–phase velocity, wavelength–phase velocity spectra, the dispersion surfaces and the slowness curves with fixed frequencies, as well as the frequency–wavenumber and frequency–phase velocity spectra. Specially, the dispersion surfaces and the slowness curves completely reflect the propagation characteristics of in-plane waves along all directions. First, the method of reverberation-ray matrix (MRRM) combined with the Floquet theorem is extended to derive the dispersion equation of in-plane elastic waves in general periodic multilayered isotropic structures by means of the elastodynamic theory of isotropic materials and the state space formalism of layers. The correctness of the derivation and the numerical stability of the method in both low and high frequency ranges, particularly its superiority over the method of the transfer matrix (MTM) within the ranges near the cutoff frequencies, are verified by several numerical examples. From these demonstrations for periodic octal- and bi-layered media, the comprehensive dispersion curves are provided and their general characteristics are summarized. It is found that although the frequencies associated with the dimensionless wavenumber along thickness ql=nπ (n is an integer) are always the demarcation between pass and stop bands in the case of perpendicular incident wave, but this is not always exist in the case of the oblique incident wave due to the coupling between the two modes of in-plane elastic waves. The slowness curves with fixed frequencies of Floquet waves in periodically multilayered isotropic structures, as compared to their counterpart of body waves in infinite isotropic media obtained from the Christoffel equation now have periodicity along the thickness direction, which is consistent to the configuration of the structures. The slowness curves associated with higher frequencies have a smaller minimum positive period and have more propagation modes due to the cutoff properties of these additional modes.