Appearance of Supersonic Stoneley Waves in Auxetics

Abstract

It is shown that in auxetic materials (materials with negative Poisson’s ratio), supersonic Stoneley waves travelling without attenuation with a velocity equal to or exceeding maximum bulk wave velocity, may exist. Analytical expressions for the relation between negative Poisson’s ratio and Young’s moduli of the contacting isotropic media ensuring the condition of propagation for supersonic Stoneley waves, are derived by solving a secular equation for Stoneley waves.

1. Introduction

It is known [1,2,3,4,5,6,7,8,9] that genuine Stoneley waves travel along an interface between two dissimilar halfspaces with subsonic velocity, which is smaller than the minimal shear bulk wave velocity in the halfspaces. However, as it is shown theoretically and numerically, the supersonic Stoneley waves may also exist when both halfspaces have a common negative Poisson’s ratio (auxetics). Propagation of the interfacial Stoneley waves in linearly elastic homogeneous isotropic and auxetic halfspaces is analyzed by constructing analytical solutions of a secular equation for Stoneley waves [1]. According to [1], the interfacial Stoneley waves propagate along an interface between two acoustically different halfspaces under ideal contact conditions. Thus, the contacting halfspaces form an inhomogeneous space. It is known [1,2,3,4] that these waves propagate at a constant speed, independent of frequency; therefore, Stoneley waves do not exhibit dispersion; and in both halfspaces, Stoneley waves exponentially attenuate with depth.

Herein, the secular equation is written in a dimensionless form close to Scholte’s secular equation [2,3] for the velocity of Stoneley waves propagating along an interface between two isotropic elastic halfspaces. While Scholte’s secular equation has been previously used for deriving the velocity of Stoneley waves, it is being used for the first time for analyzing the velocity of supersonic Stoneley waves, when acoustical properties ensure Stoneley wave velocity exceeding bulk wave velocity in one of the halfspaces.

$$\begin{array}{l}P(c)\equiv c^4\left({\left({\mathsf{\rho }}_1-{\mathsf{\rho }}_2\right)}^2-\left({\mathsf{\rho }}_1{A}_2+{\mathsf{\rho }}_2{A}_1\right)\left({\mathsf{\rho }}_1{B}_2+{\mathsf{\rho }}_2{B}_1\right)\right)\\ +2K{c}^2\left({\mathsf{\rho }}_1{A}_2{B}_2-{\mathsf{\rho }}_2{A}_1{B}_1-{\mathsf{\rho }}_1+{\mathsf{\rho }}_2\right)+K^2\left(A_1{B}_1-1\right)\left(A_2{B}_2-1\right)=0\end{array}$$

(1)

where $c$ is the Stoneley wave velocity; ${\mathsf{\rho }}_k,k=1,2$ are material densities of the contacting halfspaces, and

$$K=2\left({\mathsf{\rho }}_1{\mathsf{\beta }}_1^2-{\mathsf{\rho }}_2{\mathsf{\beta }}_2^2\right), A_k=\sqrt{1-\frac{c^2}{{\mathsf{\alpha }}_k^2}},\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{B}_k=\sqrt{1-\frac{c^2}{{\mathsf{\beta }}_k^2}},\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}k=\hspace{0.17em}\hspace{0.17em}1,2$$

(2)

In Equation (2) ${\mathsf{\alpha }}_k,\hspace{0.17em}\hspace{0.17em}{\mathsf{\beta }}_k,\hspace{0.17em}\hspace{0.17em}k=1,2$ are longitudinal and shear bulk wave velocities:

$${\mathsf{\alpha }}_k=\sqrt{\frac{{\mathsf{\lambda }}_k+2{\mathsf{\mu }}_k}{{\mathsf{\rho }}_k}},\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{\mathsf{\beta }}_k=\sqrt{\frac{{\mathsf{\mu }}_k}{{\mathsf{\rho }}_k}}$$

(3)

where ${\mathsf{\lambda }}_k,\hspace{0.17em}\hspace{0.17em}{\mathsf{\mu }}_k,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}k=1,2$ are Lame’s constants.

The majority of theoretical works on Stoneley wave analyses [1,2,3,4,5,6,7,8,9] are concerned with the case of subsonic velocities, satisfying condition

$$c<\underset{k=1,2}{\min }\left({\mathsf{\beta }}_k\right)$$

(4)

Condition Equation (4) implies all the coefficients in the Stoneley equation to be real, and, as was asserted in [1], it is necessary for the existence of a positive root of the secular Equation (1) associated with the Stoneley wave velocity. Other methods for analyzing Stoneley waves propagating with subsonic velocities obeying condition Equation (4) relate to different variants of sextic formalisms, Stroh [10,11] and Cauchy [12], along with various asymptotic techniques [13,14,15]. The principal ability for leaky Stoneley waves propagating with velocities exceeding the limit in the right-hand side of Equation (4) is asserted in [10,11].

Herein, an explicit analytical condition between physical properties of the auxetic media is derived, ensuring the existence of supersonic leaky Stoneley waves propagating with velocity

$$c=\underset{k=1,2}{\max }\left({\mathsf{\alpha }}_k\right).$$

(5)

The derivation is based on constructing analytical solution for Stoneley secular Equation (1) at a special condition

$${\mathsf{\rho }}_1{\mathsf{\beta }}_1^2={\mathsf{\rho }}_2{\mathsf{\beta }}_2^2$$

(6)

Condition Equation (6) is in some cases less restrictive than the Wiechert condition [16,17] frequently used at Stoneley wave analyses [1,2,3,4,5,6,7,8]:

$$\frac{{\mathsf{\lambda }}_1}{{\mathsf{\lambda }}_2}=\frac{{\mathsf{\mu }}_1}{{\mathsf{\mu }}_2}=\frac{{\mathsf{\rho }}_1}{{\mathsf{\rho }}_2}$$

(7)

since in view of Equation (3), condition Equation (6) imposes restriction on ${\mathsf{\mu }}_k, k=1,2$ only, demanding

$${\mathsf{\mu }}_1={\mathsf{\mu }}_2.$$

(8)

Note that condition Equation (8) does not necessarily imply ${\mathsf{\lambda }}_1={\mathsf{\lambda }}_2$. Actually, condition Equation (6) resembles the condition for equal acoustic impedances [17]. However, as will be shown below, condition Equation (6) is much less restrictive than the condition of equal acoustic impedances.

2. Solving Secular Equation

2.1. Basic Relations

Consider stratified inhomogeneous space consisting of two isotropic and homogeneous halfspaces with ideal mechanical contact between them. Introducing dimensionless Lame’s constants $\tilde{\mathsf{\lambda }},\hspace{0.17em}\hspace{0.17em}\tilde{\mathsf{\mu }}$ and dimensionless density $\tilde{\mathsf{\rho }}$

$$\tilde{\mathsf{\lambda }}=\frac{{\mathsf{\lambda }}_1}{{\mathsf{\lambda }}_2},\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\tilde{\mathsf{\mu }}=\frac{{\mathsf{\mu }}_1}{{\mathsf{\mu }}_2},\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\tilde{\mathsf{\rho }}=\frac{{\mathsf{\rho }}_1}{{\mathsf{\rho }}_2}$$

(9)

the dimensionless bulk wave velocities in the halfspaces take the form

$$\tilde{\mathsf{\alpha }}=\frac{{\mathsf{\alpha }}_1}{{\mathsf{\alpha }}_2},\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\tilde{\mathsf{\beta }}=\frac{{\mathsf{\beta }}_1}{{\mathsf{\beta }}_2}$$

(10)

with notations Equations (9) and (10), secular Equation (1) can be rewritten in terms of the dimensionless variables

$$P(\tilde{c})\equiv L{\tilde{c}}^4+2M\hspace{0.17em}{\tilde{c}}^2+N=0$$

(11)

where $\tilde{c}$ is the dimensionless velocity

$$\tilde{c}=\frac{c}{{\mathsf{\beta }}_2}$$

(12)

and

$$\begin{array}{c}L=\left({\left(\tilde{\mathsf{\rho }}-1\right)}^2-\left(\tilde{\mathsf{\rho }}{\tilde{A}}_2+{\tilde{A}}_1\right)\left(\tilde{\mathsf{\rho }}{\tilde{B}}_2+{\tilde{B}}_1\right)\right)\\ M=\tilde{K}\left(\tilde{\mathsf{\rho }}{\tilde{A}}_2{\tilde{B}}_2-{\tilde{A}}_1{\tilde{B}}_1-\tilde{\mathsf{\rho }}+1\right)\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\\ N={\tilde{K}}^2\left({\tilde{A}}_1{\tilde{B}}_1-1\right)\left({\tilde{A}}_2{\tilde{B}}_2-1\right)\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\end{array}$$

(13)

In Equation (13) coefficient $\tilde{K}$ has the form

$$\tilde{K}=2\left(\tilde{\mathsf{\rho }}{\tilde{\mathsf{\beta }}}^2-1\right)$$

(14)

and

$${\tilde{A}}_k=\sqrt{1-{\tilde{c}}^2\frac{{\mathsf{\beta }}_2^2}{{\mathsf{\alpha }}_k^2}},\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{\tilde{B}}_k=\sqrt{1-{\tilde{c}}^2\frac{{\mathsf{\beta }}_2^2}{{\mathsf{\beta }}_k^2}},\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}k=\hspace{0.17em}\hspace{0.17em}1,2$$

(15)

2.2. Secular Equation at Condition

In terms of the dimensionless parameters condition Equation (6) reads as

$$\tilde{\mathsf{\rho }}{\tilde{\mathsf{\beta }}}^2=1$$

(16)

ensuring

$$\tilde{K}=0$$

(17)

The latter yields

$$M=N=0$$

(18)

Now, at Equation (18), Equation (11) becomes

$$L{\tilde{c}}^4=0$$

(19)

Equation (19) may have nontrivial solution(s) with respect to $\tilde{c}$, if

$$L=0$$

(20)

Suppose now that both contacting media have equal Poisson’s ratios then in view of condition Equation (16), relative bulk velocities become

$$\tilde{\mathsf{\alpha }}=\tilde{\mathsf{\beta }}=\frac{1}{\sqrt{\tilde{\mathsf{\rho }}}}$$

(21)

and

$$\begin{array}{c}{\tilde{A}}_1=\sqrt{1-\frac{{\tilde{c}}^2}{\tilde{\mathsf{\rho }}{\mathsf{\gamma }}^2}},\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{\tilde{A}}_2=\sqrt{1-\frac{{\tilde{c}}^2}{{\mathsf{\gamma }}^2}}\hspace{0.17em}\hspace{0.17em}\\ {\tilde{B}}_1=\sqrt{1-\frac{{\tilde{c}}^2}{\tilde{\mathsf{\rho }}}},\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{\tilde{B}}_2=\sqrt{1-{\tilde{c}}^2}\end{array}$$

(22)

where

$$\mathsf{\gamma }=\sqrt{2}\sqrt{\frac{1-\mathsf{\nu }}{1-2\mathsf{\nu }}}$$

(23)

In Equation (23) $\mathsf{\nu }$ is the common Poisson’s ratio of both media.

3. Supersonic Stoneley Waves

Consider now the case when velocity of Stoneley waves satisfies supersonic condition Equation (5), then Equation (20) with account of Equations (21) and (22), yields

$$\{\begin{array}{c}{(\tilde{\mathsf{\rho }}-1)}^2-\tilde{\mathsf{\rho }}\sqrt{1-\tilde{\mathsf{\rho }}}\left(\tilde{\mathsf{\rho }}\sqrt{1-{\mathsf{\gamma }}^2}+\sqrt{1-\tilde{\mathsf{\rho }}{\mathsf{\gamma }}^2}\right)=0,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\tilde{\mathsf{\rho }}>1\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\&\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{\tilde{c}}^2=\tilde{\mathsf{\rho }}{\mathsf{\gamma }}^2\\ {(\tilde{\mathsf{\rho }}-1)}^2-\sqrt{1-\frac{1}{\tilde{\mathsf{\rho }}}}\left(\tilde{\mathsf{\rho }}\sqrt{1-\frac{{\mathsf{\gamma }}^2}{\tilde{\mathsf{\rho }}}}+\sqrt{1-{\mathsf{\gamma }}^2}\right)=0,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\tilde{\mathsf{\rho }}<1\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\&\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{\tilde{c}}^2={\mathsf{\gamma }}^2\end{array}$$

(24)

Note that parameter $\mathsf{\gamma }\in \left(\frac{2}{\sqrt{3}};\hspace{0.17em}\hspace{0.17em}\infty \right)$ at $\mathsf{\nu }\in (-1,\hspace{0.17em}\hspace{0.17em}0.5)$, thus the left-hand sides of Equation (24) are real, despite all imaginary radicands.

3.1. Solutions for $\mathsf{\gamma }$

Equation (24) admit analytical closed form solutions in terms of relations between $\tilde{\mathsf{\rho }}$ and $\mathsf{\gamma }$:

$$\begin{array}{c}\mathsf{\gamma }=\{\begin{array}{l}{\tilde{\mathsf{\rho }}}^{-3/2}\sqrt{\left(\tilde{\mathsf{\rho }}+1\right)\left({\tilde{\mathsf{\rho }}}^2+\tilde{\mathsf{\rho }}-1\right)-2\sqrt{\tilde{\mathsf{\rho }}\left(\tilde{\mathsf{\rho }}-1\right)\left({\tilde{\mathsf{\rho }}}^2+\tilde{\mathsf{\rho }}-1\right)}}\\ {\tilde{\mathsf{\rho }}}^{-3/2}\sqrt{\left(\tilde{\mathsf{\rho }}+1\right)\left({\tilde{\mathsf{\rho }}}^2+\tilde{\mathsf{\rho }}-1\right)+2\sqrt{\tilde{\mathsf{\rho }}\left(\tilde{\mathsf{\rho }}-1\right)\left({\tilde{\mathsf{\rho }}}^2+\tilde{\mathsf{\rho }}-1\right)}}\end{array}\\ \tilde{\mathsf{\rho }}>1\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\&\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{\tilde{c}}^2=\tilde{\mathsf{\rho }}{\mathsf{\gamma }}^2\end{array}$$

(25)

and

$$\begin{array}{c}\mathsf{\gamma }=\{\begin{array}{l}\sqrt{\left(\tilde{\mathsf{\rho }}+1\right)\left(1-{\tilde{\mathsf{\rho }}}^2+\tilde{\mathsf{\rho }}\right)-2\tilde{\mathsf{\rho }}\sqrt{\left(1-\tilde{\mathsf{\rho }}\right)\left(1-{\tilde{\mathsf{\rho }}}^2+\tilde{\mathsf{\rho }}\right)}}\\ \sqrt{\left(\tilde{\mathsf{\rho }}+1\right)\left(1-{\tilde{\mathsf{\rho }}}^2+\tilde{\mathsf{\rho }}\right)+2\tilde{\mathsf{\rho }}\sqrt{\left(1-\tilde{\mathsf{\rho }}\right)\left(1-{\tilde{\mathsf{\rho }}}^2+\tilde{\mathsf{\rho }}\right)}}\end{array}\\ \hspace{0.17em}\tilde{\mathsf{\rho }}<1\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\&\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{\tilde{c}}^2={\mathsf{\gamma }}^2\end{array}$$

(26)

3.2. Solutions for $\mathsf{\nu }$

Relations Equation (25), can also be expressed in terms Poisson’s ratio and $\tilde{\mathsf{\rho }}$ by use of Equation (23):

$$\begin{array}{c}\mathsf{\nu }=\{\frac{1}{2}\frac{{\tilde{\mathsf{\rho }}}^3-2{\tilde{\mathsf{\rho }}}^2+1+2\sqrt{\tilde{\mathsf{\rho }}\left(\tilde{\mathsf{\rho }}-1\right)\left({\tilde{\mathsf{\rho }}}^2+\tilde{\mathsf{\rho }}-1\right)}}{1-2{\tilde{\mathsf{\rho }}}^2+2\sqrt{\tilde{\mathsf{\rho }}\left(\tilde{\mathsf{\rho }}-1\right)\left({\tilde{\mathsf{\rho }}}^2+\tilde{\mathsf{\rho }}-1\right)}}\\ \hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\tilde{\mathsf{\rho }}>1\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\&\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{\tilde{c}}^2=\tilde{\mathsf{\rho }}{\mathsf{\gamma }}^2\end{array}$$

(27a)

$$\begin{array}{c}\mathsf{\nu }=\{\frac{1}{2}\frac{{\tilde{\mathsf{\rho }}}^3-2{\tilde{\mathsf{\rho }}}^2+1-2\sqrt{\tilde{\mathsf{\rho }}\left(\tilde{\mathsf{\rho }}-1\right)\left({\tilde{\mathsf{\rho }}}^2+\tilde{\mathsf{\rho }}-1\right)}}{1-2{\tilde{\mathsf{\rho }}}^2-2\sqrt{\tilde{\mathsf{\rho }}\left(\tilde{\mathsf{\rho }}-1\right)\left({\tilde{\mathsf{\rho }}}^2+\tilde{\mathsf{\rho }}-1\right)}}\\ \hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\tilde{\mathsf{\rho }}>1\hspace{0.17em}\hspace{0.17em}\&\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{\tilde{c}}^2=\tilde{\mathsf{\rho }}{\mathsf{\gamma }}^2\end{array}$$

(27b)

and

$$\begin{array}{c}\mathsf{\nu }=\{\frac{1}{2}\frac{{\tilde{\mathsf{\rho }}}^3-2\tilde{\mathsf{\rho }}+1+2\tilde{\mathsf{\rho }}\sqrt{\left(1-\tilde{\mathsf{\rho }}\right)\left(1+\tilde{\mathsf{\rho }}-{\tilde{\mathsf{\rho }}}^2\right)}}{{\tilde{\mathsf{\rho }}}^2-2+2\sqrt{\left(1-\tilde{\mathsf{\rho }}\right)\left(1+\tilde{\mathsf{\rho }}-{\tilde{\mathsf{\rho }}}^2\right)}}\\ \hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\tilde{\mathsf{\rho }}<1\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\&\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{\tilde{c}}^2={\mathsf{\gamma }}^2\end{array}$$

(28a)

$$\begin{array}{c}\mathsf{\nu }=\{\frac{1}{2}\frac{{\tilde{\mathsf{\rho }}}^3-2\tilde{\mathsf{\rho }}+1-2\tilde{\mathsf{\rho }}\sqrt{\left(1-\tilde{\mathsf{\rho }}\right)\left(1+\tilde{\mathsf{\rho }}-{\tilde{\mathsf{\rho }}}^2\right)}}{{\tilde{\mathsf{\rho }}}^2-2-2\sqrt{\left(1-\tilde{\mathsf{\rho }}\right)\left(1+\tilde{\mathsf{\rho }}-{\tilde{\mathsf{\rho }}}^2\right)}}\\ \hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\tilde{\mathsf{\rho }}<1\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\&\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{\tilde{c}}^2={\mathsf{\gamma }}^2\end{array}$$

(28b)

3.3. Numerical Analysis

All four solutions in Equations (27) and (28) lead to auxetics (materials with negative Poisson’s ratio [18]), in particular, solution (27a) leads to negative Poisson’s ratios at $\tilde{\mathsf{\rho }}\in (1;\hspace{0.17em}\hspace{0.17em}~1.23)$; solution (27b) leads to negative Poisson’s ratios at $\tilde{\mathsf{\rho }}\in (~3.8;\hspace{0.17em}\hspace{0.17em}~12.0)$; solution (28a) leads to negative Poisson’s ratios at $\tilde{\mathsf{\rho }}\in (~0.82;\hspace{0.17em}\hspace{0.17em}1)$; and solution (28b) yields negative Poisson’s ratios at $\tilde{\mathsf{\rho }}\in (~0.09;\hspace{0.17em}\hspace{0.17em}0.25)$.

The plots revealing the variation of Poisson’s ratio vs relative material density $\tilde{\mathsf{\rho }}$ corresponding to solutions (27a) and (27b) are shown in Figure 1.

The plots obtained by Equations (27a) and (27b) are shown in Figure 2.

Thus, at both $\tilde{\mathsf{\rho }}<1$ and $\tilde{\mathsf{\rho }}>1$ there are intervals of relative material densities admitting propagation of non-attenuating leaky Stoneley waves for the considered auxetic media with common negative Poisson’s ratios; in this respect, see also [19,20,21]. It should also be noted that according to [22], auxetic materials need not be necessary anisotropic, but they may belong to the isotropic class. Another remark concerns some of the related works on acoustic wave propagation in media with more complicated properties; in this respect, see [23] for nonlinear waves and [24] for Stoneley waves in media satisfying the Wiechert condition.

4. Concluding Remarks

The obtained analytical expressions Equations (27) and (28) for the common Poisson’s ratios of the auxetic media and their densities, define supersonic leaky Stoneley waves propagating without attenuation in the direction of propagation along the interface plane, due to real values for the Stoneley wave velocity.

In contrast to the non-attenuating (in the direction of propagation) leaky Rayleigh waves that arise only at special kinds of elastic anisotropy [10,11], the analyzed supersonic and non-attenuating leaky Stoneley waves propagating with velocities

$$\{\begin{array}{ll}{\tilde{c}}^2=\tilde{\mathsf{\rho }}{\mathsf{\gamma }}^2& \hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}@\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\tilde{\mathsf{\rho }}>1\\ {\tilde{c}}^2={\mathsf{\gamma }}^2& \hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}@\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\tilde{\mathsf{\rho }}<1\end{array}$$

(29)

exist at the contacting isotropic auxetic halfspaces, satisfying condition Equation (16). Also note that at the intermediate velocity range

$$\{\begin{array}{ll}\tilde{\mathsf{\beta }}<{\tilde{c}}^2<\tilde{\mathsf{\rho }}{\mathsf{\gamma }}^2& \hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}@\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\tilde{\mathsf{\rho }}>1\\ 1<{\tilde{c}}^2<{\mathsf{\gamma }}^2& \hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}@\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\tilde{\mathsf{\rho }}<1\end{array}$$

(30)

non-attenuating in direction of propagation leaky Stoneley waves cannot exist because of appearing imaginary terms in secular Equation (24), leading to complex values for velocity $\tilde{c}$, and hence, attenuation in direction of propagation.