Abstract
This study is devoted to the analysis of the band structures of the anti-plane transverse wave (SH wave) in nanoscale layered metamaterial structures. Attention is restricted to normal incidence of waves. The localization factor is introduced to characterize the band structures. The general transfer matrix method based on the nonlocal elastic continuum theory is employed to calculate the localization factor. Based on the analysis of band structures, the influences of random disorder of the internal characteristic length and the external thickness of each sub-layer, the aperiodic arrangements, the location of different material components, the ratio of mass density, the ratio of the transverse wave velocity, the ratio of the internal characteristic length or the external thickness of each sub-layer on the band structures, the cut-off frequency, the peak points and the dense band zones are investigated and discussed in detail, which can provide some new thoughts for the designs and applications of the nanoscale wave devices.
1 Introduction
The metamaterials, phononic crystals (PCs) (), have been studied intensely over the past 2 decades due to their potential capability of controlling and tuning the propagation of acoustic/elastic waves. These metamaterials have band gap characteristics, that is, waves in the band gap frequency range are prohibited from passing through these structures. The unusual effect of PCs with band gaps has a wide range of potential important applications such as sound detectors, transducers, filters, waveguides, sensors, etc. Compared with two-dimensional (2D) and three-dimensional (3D) PCs, one dimensional (1D) layered PCs have simpler structure and can fully show the characteristics of wave propagation, thus, many experimental and theoretical researches on the band structures of 1D macroscale layered PCs are witnessed over the past decades (; ; ; ; ; ; ). The PCs are generally periodic. However, the random disorder (; ; ) and quasi-periodic arrangement (; ; ; ; ), may exhibit unique characteristics of a mixture of acoustic/elastic wave propagation and localization, which are of significant interest in both basic and applied sciences (). Although the macroscale quasi-periodic or aperiodic phononic crystals (APNCs) have been extensively investigated and reported in literature (; ; ; ; ), (; ; ; ; ; ), very little theoretical study on the band structures of nanoscale APNCs has been performed.
In recent years, owing to the wide potential applications in new thermo-elecrical, acousto-optical, nanoscale electro-mechanical devices and computer chips (; ; ), more and more researchers have conducted extensive research on nanoscale structures. As we all know, when the structure size is several namometers, the size effect must be considered. In this case, the constitutive relationships cannot be described accurately by the conventional elastic continuum theory (; ). Therefore, many methods have been developed to study the mechanical behaviors of nanoscale materials and structures (; ; ; ; ; ; ; ; ; ), in which the nonlocal elastic (NLE) continuum theory proposed by Eringen (; ) can describe the long-range inter-atomic interactions and can account for the nanoscale size effect inside the structures. By utilizing the NLE continuum theory, Artan et al. () studied the effect of nonlocality on the dynamic behavior of laminated composites by means of dispersion of SV waves propagating in the direction parallel to layering. studied the sound wave propagation in single-walled carbon nanotubes using NLE continuum theory, and revealed the significance of the small-scale effect on wave propagation in single-walled carbon nanotubes. presented the paradoxes in the existing solutions of the nonlocal field equation by introducing the high-order boundary conditions. investigated the nonlinear vibration of the piezoelectric nanobeams based on the NLE continuum theory and Timoshenko beam theory. And the influences of the nonlocal parameter, temperature change and external electric voltage on the size-dependent nonlinear vibration characteristics of the piezoelectric nanobeams are conducted. studied the propagation of Love waves in an isotropic homogeneous elastic medium in the frame of the NLE continuum theory, and determined the nonlocal modulus by comparing the dispersion equation of the plane transverse waves with the corresponding equation given by the atomic lattice dynamics. analyzed the vibration of a nano-plate based on the NLE continuum theory. In addition, the nonlocality also plays an important role in electronic and magnetic materials (; ; ; ; ; ). For example, presented an exact closed-form solution for the three-dimensional free vibrational response of a simply-supported and multilayered magneto-electro-elastic plate considering the nonlocal effect. derived the analytical solutions for propagation of time-harmonic waves in three-dimensional magneto-electro-elastic multilayered plates with nonlocal effect, and investigated the influences of the nonlocal parameter on the dispersion curves. It should be noted that by developing the transfer matrix method based on the NLE continuum theory (; ; ; ; ; ), a series of extensive studies on wave propagation in nanoscale periodic structures have been carried out. The results showed that a cut-off frequency was found, beyond which the waves are prohibited from passing through the structure. Besides, the dense band zones (DBZs) appeared in the band structures when the nanoscale size-effect is taken into account. However, the PCs in the above studies are all perfect periodic. For nearly periodic nanoscale layered PCs, studied the size effect on the band structures of randomly disordered, quasi-periodic and defected nanoscale PCs. Therein, only the disorder of the external thickness of the first sub-layer is considered for simplicity. Besides, only Fibonacci sequence is studied. However, the influences of random disorder of the internal characteristic length and the external thickness of each sub-layer, the aperiodic arrangements, the location of different material components, the ratio of the mass density, the ratio of the transverse wave velocity, the ratio of the internal characteristic length or the external thickness of each sub-layer on the band structures, the cut-off frequency, the peak points and the DBZs have not been investigated, which requires a detailed study of these problems.
In this paper, we attempt to address these questions and the band structures of the SH wave in the nanoscale layered structures are studied in detail. The general transfer matrix method based on the NLE continuum theory is used to calculate the localization factor describing the band structures. A detailed parametric study is conducted to investigate the influences of random disorder of the internal characteristic length and the external thickness of each sub-layer, the aperiodic arrangements, the location of different components, the ratio of the mass density, the ratio of the transverse wave velocity, the ratio of the internal characteristic length or the external thickness of each sub-layer on the band structures, the cut-off frequency, the peak points and the DBZs.
The paper is structured as follows: Section 2 introduces the nonlocal elastic continuum theory. And the theoretical models and the general transfer matrix method are given in Section 3. Section 4 is devoted to the illustration and discussion of the results based on the calculations of the localization factor, where different influence factors are taken into account. Finally, some conclusions and future perspectives are presented in Section 5.
2 The Nonlocal Elastic Continuum Theory
In nonlocal elastic theory, owing to the long-range interaction between atoms or molecules in nanoscale materials and structures, the stresses at a point are related not only to the strains at the same point, but also to the strains at other points of the whole body. The nonlocal elastic continuum model proposed by Eringen (; ) well explains that the physical phenomenon represented by one point in the continuum is affected by all other points in the whole domain, and the results are consistent with the experimental observations of lattice atomic dynamics and phonon scattering. For homogeneous, isotropic and elastic solids, the nonlocal and classical stress tensor has the following relationship, which includes an integral involving the whole region, i.e., ().where represents the nonlocal stress tensor at point , the kernel function is the nonlocal modulus, is the Euclidean distance, is the internal characteristic length, is the classical stress tensor whose components can be defined bywith the strain componentswhere and are classical Lamé constants, is the Kronecker-delta, and are the displacement components, respectively.
The kernel function in Eq. 1 which depends on the internal characteristic length can be determined by matching the dispersion curves with those obtained from atomic lattice dynamics, first principle method and experiments. Because the structures considered in this paper are infinite along the direction, the kernel function along the direction is supposed to be a Delta function, and thenwhere is the kernel function along the direction. Therefore, Eq. 1 can be written as the following component form
It is well known that the kernel function has different forms (). Considering the time-harmonic elastic waves in this paper, it is more suitable to choose as (), i.e., exponential kernel function
In Eringen’s NLE theory, the integral form of Eq. 1 can be approximated in the following differential form ():where is the Laplace operator. The equations of wave motion without body forces are expressed aswith being the mass density. Here, the repeated indices denote the conventional summation rule. Substituting Eqs 2, 3, and 7 into Eq. 8, the wave motion equation based on the NLE theory can be written as the following differential ones
3 Theoretical Model and the General Transfer Matrix Method
The SH wave propagating normally in the nanoscale periodic structures are presented in Figure 1A. Figures 1B–D shows the schematic diagrams of the nanoscale layered PCs arranged as Thue-Morse sequence, Rudin-Shapiro sequence and Fibonacci sequence, respectively, and the random disorder, quasi-periodic and aperiodic structures considered in this paper can be obtained by the following theoretical model designs.
FIGURE 1
3.1 Theoretical Model
3.1.1 Nanoscale Random Disordered Structure
Here, we consider the elastic SH waves propagating in normally distributed randomly disordered PCs. let denote the internal characteristic length and the external thickness of the each sub-layer, respectively. For the normally distributed randomly disordered PC, the characteristic length can be written aswhere is the mean of (corresponding to the perfect periodic distribution), and is the variance of the internal characteristic length or the external thickness representing the disorder degree of this system, corresponds to a perfect periodic system. , in which are standard uniformly distributed random variables, .
3.1.2 Nanoscale Quasi-Periodic Structure
Here, we consider the nanoscale quasi-periodic layered structures arranged in the Fibonacci sequence () as shown in Figure 1D. The Fibonacci sequence can be obtained by repeating operations of the concurrent substitution rules: and (). The generation of the Fibonacci sequence is denoted as with and . Then the Fibonacci sequence can be written as , for example, where and are sub-layers made up of different materials.
3.1.3 Nanoscale Aperiodic Structures
Being a bridge of linking periodic models with quasi-periodic systems in a geometrical structure, Thue-Morse system () and Rudin-Shapiro systems illustrated in Figures 1B,C are thought to be more random than the quasi-periodic Fibonacci lattices.
The Thue-Morse sequence is based on the two letter alphabet , and can be generated by the inflation rules, as follows: , . The lower-order Thue-Morse lattices are the strings , etc.
The Rudin-Shapiro sequence is an infinite sequence and can be generated by a four state automaton as follows: , , , , respectively. It should be noticed that for the two aperiodic systems mentioned above, the letters and denote two different material sub-layers.
3.2 General Transfer Matrix Method
In this section, we start from the periodic system, i.e., the normal propagation of the time-harmonic SH elastic waves in a nanoscale periodic multilayered phononic crystal is considered. This structure depicted in Figure 1A is composed of materials and layers, with the thicknesses and , respectively, and is the thickness of one unit-cell. The local coordinates of the monolayers are also given in the figure. We assume that the layered composite consists of unit-cells. Each unit-cell includes, unless otherwise stated, two sub-layers which are denoted by the subscript . For this problem, the displacement components in the - and -directions, i.e., and are zero, and the only non-zero displacement is along the -direction, which is perpendicular to the -plane. Then Eq. 9 can be simplified into the following form for the sub-layer
By introducing the dimensionless local coordinate , , Eq. 11 can be rewritten into the following dimensionless formwhere is the non-dimensional frequency with being the velocity of the transverse elastic wave, and is the ratio of the internal characteristic length and the external thickness of the unit-cell.
Then, the general harmonic solution for the sub-layer can be obtained, which has the following form:where , and are the unknown coefficients to be determined. According to Eqs 13, 5, the nonlocal stresses can be obtained as
Then, the state vector with the bar denoting the dimensionless parameters is chosen in order to obtain the transfer matrix of the unit-cell. Based on the continuity conditions on the left (subscript ) and right (subscript ) sides of the interface of two adjacent sub-layers in the unit cell, the state vectors can be defined aswhere the matrices and can be obtained from Eqs 13, 5, whose elements are given by
Obviously, the two state vectors in Eq. 15 have the following relation by eliminating the common vector, i.e.,where is the transfer matrix of the sub-layer.
The displacements and nonlocal stresses are continuous at the interface of two adjacent sub-layers in the same unit-cell and between the and the unit-cells, that is,
From Eqs 17, 18, the following relation can be obtainedwhich shows the relationship between the state vectors of the and the unit-cells, where is the transfer matrix between the two consecutive unit-cells, i.e., the transfer matrix of the unit-cell. For perfect periodic two-component PCs, for all are the same and denoted as . It should be noticed that the above derivation is applicable for not only the ordered periodic PNCs but also the disordered, quasi-periodic and aperiodic ones. However the transfer matrices of the “unit-cells” of the disordered, quasi-periodic and APNCs are different from those of the perfect periodic ones. For quasi-periodic or aperiodic structures, are not all the same and the Bloch theory is not applicable, for example, Using the above method, for the aperiodic structure composed of unit-cells, the total transfer matrix can be obtained, that is,
The detailed mathematical derivation is not given here for the sake of brevity. Interested readers may refer to many publications for more details.
In this paper, we use the well-defined localization factor to characterize the band structures and localization phenomenon of 1D nanoscale layered PCs. The localization factor is defined as the minimum positive Lyapunov exponent which describes the average exponential rate of growth or attenuation of the wave amplitude (). And it can be calculated by using the Wolf’s method () once the transfer matrix is obtained. If the dimension of the transfer matrices is , then the smallest positive Lyapunov exponent is the localization factor. The expression for the localization factor of the system with unit-cells is given as follows:where the vector in Eq. 21 is given byin which are orthogonal unit vectors, denotes the dot-product, is the vector norm, and represents the number of the unit-cells. The orthogonal unit state vector is obtained through the iteration using the Gram–Schmidt orthonormalization procedures (). If the localization factor is equal to zero, the corresponding frequency intervals are known as pass-bands. Otherwise if the localization factor is positive, the frequency intervals are known as stop-bands or band-gaps. In this paper, only the normal incidence of SH wave is considered, thus, the dimension of the transfer matrix is and is the localization factor, which is denoted as in the following analysis.
4 Numerical Results and Discussions
In this section, the band structures and localization properties of the anti-plane elastic waves propagating normally in nanoscale layered structures are studied by the general transfer matrix method. Different factors affecting the band structures are considered. Numerical results are presented and discussed. During the calculations, we refer to Ref. () for the material constants and list the values in Table 1. For convenience, the frequency is normalized as with .
TABLE 1
| Component materials | HfO2 | ZrO2 | Al | Cu |
|---|---|---|---|---|
| Mass density (Kg/m3) | ρ = 10,873 | ρ = 6,488 | ρ = 2,730 | ρ = 8,950 |
| Shear modulus (Pa) | µ = 6.60 × 1010 | µ = 6.88 × 1010 | µ = 2.87 × 1010 | µ = 7.53 × 1010 |
Material constants.
4.1 Nanoscale Random Disordered Layered PCs
Firstly, In order to check the correctness of the present method, the 1D nanoscale periodic layered structures arranged alternately by HfO2 () and ZrO2 (), as shown in Figure 1A, are studied. The same material constants as the Ref. () are selected, and the results are shown in Figure 2A. We find that the present results (black solid lines) are in good agreement with the results (red solid lines) of Figure 2C in Ref. (), which verifies the correctness and effectiveness of the current method. In the following, In order to reveal the influence of random disorder of the sub-layers on the band structures by considering the internal characteristic lengths and the external thicknesses, the localization factors in the cases of (i.e., disorder is considered only for the internal characteristic length of the second sub-layer ) and , (i.e., disorder is considered for both internal characteristic lengths and external thickness of sub-layers and ), are calculated and the results are illustrated in Figures 2B,C. The disorder degree denoted by is assumed to be and , respectively. It can be seen that there is one peak point in Figure 2B and two peak points in Figure 2C whose positions are determined by the cut-off frequency of a specific material sub-layer. The second peak point in Figure 2C is a cut-off frequency beyond which the localization factor becomes stable and positive, which means that the elastic waves cannot propagate through the structure over the cut-off frequency. In addition, the dense band zones (DBZs) as defined in Ref () appears, where the localization factors are very big with multiple, dense, flat and narrow band-gaps in the frequency range in Figure 2B and in Figure 2C, meaning a very strong wave localization phenomenon. The disorder degree has little effect on the DBZ, the cut-off frequency and the localization factor whose frequencies are larger than the peak points and . However, when the disorder caused by the sub-layer’s length or thickness is introduced to the periodic phononic crystals, the disorder degree still has some influences on the band structures in both pass-bands and band gaps, which are demonstrated by the enlarged sub-figure in Figures 2B,C. For example, in the pass-band for shown in the partly enlarged regions in Figure 2B, the localization factors become positive with the increase of , indicating a wave localization phenomenon, and the localization degree increases with the increase of in the pass-bands while the localization degree decreases with the increase of in the band gaps. Next, the disorder of the internal characteristic length and the external thickness of the same sub-layer, the disorder of different sub-layers and the number of disorder are considered to show more detailed results in Figures 3A–C, respectively. The disorder degree is . It can be seen from Figure 3A that the influences of the disorder caused by the internal characteristic lengths are a little bigger than those caused by the external thicknesses. In addition, in Figures 3B,C the localization factors change a lot and have a strong dependence on the disorder of different sub-layers and the number of disorder. Moreover, with the increase of the number of disorder, the degree of wave localization in the pass-bands increases, and the localization phenomenon becomes more and more obvious. From Figure 3, we can conclude that the disorder degree and the localization factor depend on different types of disorders of each sub-layer.
FIGURE 2
FIGURE 3
4.2 Nanoscale Aperiodic Layered PCs
Due to the lack of periodicity, a finite but sufficiently large number of unit-cell is very important to calculate the localization factor. After trial calculation, is taken in the following computations. Next, the influences of the material combinations on the band structures are analyzed by considering changing the order of material components, only the single sub-layer material and two sub-layers are both changed, respectively. Specifically, the order of material components is selected as , and , . Only the change of the second sub-layer material is taken as , and . The examples where the materials of two sub-layers are both changed are , and . The band structures for the nanoscale Thue-Morse laminate with the change in order of material components are plotted in Figure 4. The ratio of internal characteristic length and the thickness of the unit-cell is . It can be seen that with the exchange of material components, the distance between the two peak points becomes narrower, and the two peak points and the first distinct band gap move to the low frequency zone. Furthermore, Compared with the two DBZs (5.42, 5.56) and (6.53, 7.34) for , the first and the second DBZs move left to (4.1, 4.21) and (4.98, 5.56) for , respectively. While the two DBZs for move left to (4.61, 4.72) and (5.3, 5.56) compared with the corresponding DBZs (5.31, 5.56) and (5.9, 6.54) for . It is worth noted that compared with the first peak points for and , the second peak points, i.e., the cut-off frequencies, stay in the almost same position for and , that is to say, the cut-off frequencies for and are almost which are exactly the right edges of the corresponding DBZs, i.e., the cut-off frequency of the material . In addition, when the frequency is higher than the cut-off frequency, the value of the localization factor will become stable at approximately 3.23 (as illustrated by the dashed lines in Figure 4).
FIGURE 4
In the following, we only change the material of the second sub-layer and the localization factors for the nanoscale aperiodic Thue-Morse sequences with different material combinations are presented in Figure 5. Here, the material combinations is taken as , and , respectively. It can be observed that the peak points, i.e., the DBZs with the frequency ranges of (5.3, 5.56), (6.48, 6.54), (7.03, 7.35) and (7, 7.32) correspond to materials , , and , respectively. The band structures between the two DBZs have almost the same distributions, however, when the frequency is lower than the first peak point, the localization factors change a lot. More specifically, the first distinct band gap becomes wider and moves to the low frequency zone. Compared with Figures 5B,C shows the two DBZs are located at almost the same position, because the transverse wave velocities of the two materials and , have almost similar values. In addition, the examples where the materials of two sub-layers are both changed are chosen as , and , respectively. It can be seen from Figure 6 that the position of the peak point shown by the dashed line corresponds to materials and because the shear wave velocities of the two materials have the similar values. The larger the difference in transverse wave velocity, the farther the two DBZs are from each other.
FIGURE 5
FIGURE 6
Next, the band structures of the systems with different aperiodic arrangements are calculated. As shown in Figures 1B–D, three different aperiodic arrangements are chosen, i.e., the Thue-Morse sequence, the Rudin-Shaprio sequence and the Fibonacci sequence, respectively. Figure 7 shows the influences of aperiodic arrangements on localization factors of the SH wave propagating normally in the nanoscale aperiodic systems consisting of for . It can be seen that the tendencies of all curves are coincident for the three aperiodic systems, which implies the aperiodic arrangements have little effect on the peak points, the DBZs and the cut-off frequency. For example, the localization factors oscillate quickly with big values between the two peak points. The cut-off frequencies are all around about . However, the band structures before the first peak point becomes different for the three aperiodic arrangements. i.e., the main band gaps of the Fibonacci structure change a lot with the gradual disappearing and narrowing of the band gaps. Specifically, in the dot rectangle, compared with the Fibonacci structure, the localization factor is more like a defect state for the Thue-Morse and Rudin-Shapiro structures.
FIGURE 7
Additionally, the influences of the ratio of the mass density and the ratio of the transverse wave velocity on the localization factors are examined. Here, the nanoscale Fibonacci aperiodic laminate is selected as an example. It can be observed from Figure 8 that when the ratio of the mass density is not equal to , the first distinct band gap emerges, while for , it disappears (enlarged regions shown in Figure 8). The positions of the two DBZs, the two peak points, the cut-off frequency and the localization factors whose frequencies are larger than those of the first peak point have no changes. For example, after the first peak point, the tendencies of all curves are coincident, i.e., the localization factors oscillate quickly with big values after the first peak point, and then the second peak point appears. More detailed results are illustrated in Figure 9 for the localization factors varying with the normalized frequency and . Compared with the results in Figure 8, the first distinct band gap disappears for , when is less than 1 or is bigger than 1, the first distinct band gap appears and becomes wider with decreasing or increasing. Combining Figures 8, 9, we can conclude that the first distinct band gap will emerge for . and the ratio of the mass density has no influences on the peak points, the cut-off frequency, the DBZs and the localization factors whose frequencies are larger than those of the first peak point. From Figure 10, we can see that there is only one peak point when the ratio of transverse wave velocity is , i.e., the cut-off frequency, and when the velocity ratio changes, the position of the peak point remains unchanged. However, when the ratio is not equal to , there are two peak points. Among them, when the ratio of velocity is less than , one peak point tends to the high-frequency region, and when the ratio is larger than , a peak point tends to the low-frequency region. In addition, the low-frequency band gap is more likely to emerge with the increase of the ratio. This shows that the ratio of transverse wave velocity has significant effect on the cut-off frequency, the DBZ, the peak point and the low-frequency band gap. Furthermore, the localization factors varying with normalized frequency and are presented in Figure 11. There is only one peak point when the ratio of transverse wave velocity is , i.e., the cut-off frequency of material Al, and when the velocity ratio changes, the position of the peak point remains unchanged since is a fixed value. On the contrary, the peak point other than the cut-off frequency moves the high-frequency region when is less than 1, and moves the low-frequency zone when is larger than 1. Combining Figures 10, 11, we can conclude that the position and the number of the peak point, i.e., the DBZ have strong dependence on . There exists a wide pass-band when is less than 1, and the pass-band becomes wider with the decrease of .
FIGURE 8
FIGURE 9
FIGURE 10
FIGURE 11
Finally, the influences of the structural parameter, i.e., the ratio of external characteristic thickness and the ratio of internal characteristic length on the localization factors are investigated. The localization factors for the nanoscale Fibonacci laminates with different values of are presented in Figure 12, From Figure 12, we can see that the pass-bands and band-gaps are affected by varying the ratio of external characteristic thickness, while the two peak points, the two DBZs and the cut-off frequency remain in the same position. Detailed results are shown in Figure 13, the position of the cut-off frequency remains unchanged, but the localization degree enhances with the increase of . From Figure 14, we can see that the bigger the difference of internal characteristic length between the two materials, the farther the distance between the two peak points. The first distinct band gap disappears with the increase of . More detailed numerical results are shown in Figure 15, It can be seen that the cut-off frequency does not appear when , However, when , the cut-off frequency appears and decreases with the increasing. When , the cut-off frequency remains unchanged and the first peak point tends to the low-frequency zone with the increasing. In addition, the localization degree of the cut-off frequency, the pass-bands and the band gaps varies with the . And the localization degree whose frequencies are larger than the cut-off frequency decreases with increasing.
FIGURE 12
FIGURE 13
FIGURE 14
FIGURE 15
5 Conclusion
The results presented in this study are obtained by the numerical calculations of the wave localization properties in the nanoscale layered structures by using the general transfer matrix method based on the nonlocal elastic continuum theory. The key conclusions from this analysis can be summarized as follows:
1) No matter what kind of disorders, the disorder degree has little effect on the DBZ, the peak point, the cut-off frequency and the localization factor whose frequencies are larger than the cut-off frequency. However, the influences of the disorder caused by the internal characteristic lengths are a little bigger than those caused by the external thicknesses. The localization factors have a strong dependence on the disorder of different sub-layers and the number of disorder.
2) The first distinct band gap, the peak point, the cut-off frequency, the DBZ and the localization factor have strong dependence on the material combinations. With the exchange of material components, the distance between the two peak points becomes narrower, and the two peak points and the first distinct band gap move to the low frequency zone. The position of the peak point and the DBZ depends on the sub-layer materials under consideration.
3) The aperiodic arrangements have little effect on the peak points, the DBZs and the cut-off frequency. However, the band structures before the first peak point becomes different for the three aperiodic arrangements.
4) The first distinct band gap will emerge for . and the ratio of the mass density has no influences on the peak points, the cut-off frequency, the DBZs and the localization factors whose frequencies are larger than those of the first peak point.
5) There is only one peak point when the ratio of transverse wave velocity is . However, when the ratio of velocity is less than , one peak point tends to the high-frequency region, and when the ratio is larger than , a peak point tends to the low-frequency region. In addition, the low-frequency band gap is more likely to emerge with the increase of the ratio. There exists a wide pass-band when is less than 1, and the pass-band becomes wider with the decrease of .
6) The pass-bands and band gaps are affected by varying the ratio of external characteristic thickness, while the two peak points, the two DBZs and the cut-off frequency remain in the same position. In addition, the localization degree enhances with the increase of .
7) When , the cut-off frequency appears and decreases with the increasing. When , the cut-off frequency remains unchanged and the first peak point tends to the low-frequency zone with the increasing, And the localization degree whose frequencies are larger than the cut-off frequency decreases with increasing.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
The authors greatly acknowledge the financial support from the National Natural Science Foundation of China (No. 11002026, 11372039), Beijing Natural Science Foundation (No. 3133039), and the Scientific Research Foundation for the Returned Overseas Chinese Scholars (No. 20121832001).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AdhikariS.GilchristD.MurmuT.McCarthyM. A. (2015). Nonlocal normal Modes in Nanoscale Dynamical Systems. Mech. Syst. Signal Process.60-61, 583–603. 10.1016/j.ymssp.2014.12.004
2
AifantisE. C. (1999). Strain Gradient Interpretation of Size Effects. Int. J. Fract95 (1-4), 299–314. 10.1007/978-94-011-4659-3_16
3
AlibeiglooA. (2011). Free Vibration Analysis of Nano-Plate Using Three-Dimensional Theory of Elasticity. Acta Mech.222, 149–159. 10.1007/s00707-011-0518-7
4
AndersonP. W. (1958). Absence of Diffusion in Certain Random Lattices. Phys. Rev.109, 1492–1505. 10.1103/physrev.109.1492
5
ArtanR.AltanB. S. (2002). Propagation of SV Waves in a Periodically Layered media in Nonlocal Elasticity. Int. J. Sol. Structures39, 5927–5944. 10.1016/s0020-7683(02)00476-6
6
AynaouH.BoudoutiE. H. E.Djafari-RouhaniB.AkjoujA.VelascoV. R.VelascoV.R (2005). Propagation and localization of acoustic waves in Fibonacci phononic circuits. J. Phys.: Condens. Matter17, 4245–4262. 10.1088/0953-8984/17/27/002
7
BarcoOOrtunoM (2012). Localization length of nearly periodic layered metamaterials. Phys Rev A86 (2), 023846. 10.1103/physreva.86.023846
8
BovierA.GhezJ.-M. (1995). Remarks on the Spectral Properties of Tight-Binding and Kronig-Penney Models with Substitution Sequences. J. Phys. A: Math. Gen.28, 2313–2324. 10.1088/0305-4470/28/8/022
9
ChenA.-L.TianL.-Z.WangY.-S. (2017). Band Structure Properties of Elastic Waves Propagating in the Nanoscaled Nearly Periodic Layered Phononic Crystals. Acta Mechanica Solida Sinica30, 113–122. 10.1016/j.camss.2017.03.005
10
ChenA.-L.WangY.-S.KeL.-L.GuoY.-F.WangZ.-D. (2013). Wave Propagation in Nanoscaled Periodic Layered Structures. Jnl Comp. Theo Nano10, 2427–2437. 10.1166/jctn.2013.3225
11
ChenA.-L.WangY.-S. (2011). Size-effect on Band Structures of Nanoscale Phononic Crystals. Physica E: Low-dimensional Syst. Nanostructures44, 317–321. 10.1016/j.physe.2011.08.032
12
ChenA.-L.WangY.-S. (2007). Study on Band Gaps of Elastic Waves Propagating in One-Dimensional Disordered Phononic Crystals. Physica B: Condensed Matter392, 369–378. 10.1016/j.physb.2006.12.004
13
ChenA.-L.WangY.-S.ZhangC. (2012). Wave Propagation in One-Dimensional Solid-Fluid Quasi-Periodic and Aperiodic Phononic Crystals. Physica B: Condensed Matter407, 324–329. 10.1016/j.physb.2011.10.041
14
ChenA.-L.YanD.-J.WangY.-S.ZhangC. (2016). Anti-plane Transverse Waves Propagation in Nanoscale Periodic Layered Piezoelectric Structures. Ultrasonics65, 154–164. 10.1016/j.ultras.2015.10.006
15
ChenA.-L.YanD.-J.WangY.-S.ZhangC. (2019). In-plane Elastic Wave Propagation in Nanoscale Periodic Piezoelectric/piezomagnetic Laminates. Int. J. Mech. Sci.153-154, 416–429. 10.1016/j.ijmecsci.2019.02.017
16
ChenJ.-J.WangQ.HanX. (2010). Lamb Wave Transmission through One-Dimensional Three-Component Fibonacci Composite Plates. Mod. Phys. Lett. B24, 161–167. 10.1142/s0217984910022305
17
ChenJ.GuoJ.PanE. (2017). Wave Propagation in Magneto-Electro-Elastic Multilayered Plates with Nonlocal Effect. J. Sound Vibration400, 550–563. 10.1016/j.jsv.2017.04.001
18
DuS. G.ShiD. M.DengH. (2000). Special effects and applications of nanostructured materials. Ziran Zazhi22 (2), 101–106.
19
EconomouE. N.SigalasM. (1994). Stop bands for elastic waves in periodic composite materials. The Journal of the Acoustical Society of America95, 1734–1740. 10.1121/1.408692
20
El-NabulsiR. A. (2018). Nonlocal Approach to Energy Bands in Periodic Lattices and Emergence of Electron Mass Enhancement. J. Phys. Chem. Sol.122, 167–173. 10.1016/j.jpcs.2018.06.028
21
El-NabulsiR. A. (2018). On Nonlocal Complex Maxwell Equations and Wave Motion in Electrodynamics and Dielectric media. Opt. Quant Electron.50, 170. 10.1007/s11082-018-1436-x
22
EringenA. C. (2006). Nonlocal Continuum Mechanics Based on Distribution. Int. J. Eng. Sci.44 (3-4), 141–147. 10.1016/j.ijengsci.2005.11.002
23
EringenA. C. (1972). Nonlocal Polar Elastic Continua. Int. J. Eng. Sci.10 (1), 1–16. 10.1016/0020-7225(72)90070-5
24
EringenA. C. (1983). On Differential Equations of Nonlocal Elasticity and Solutions of Screw Dislocation and Surface Waves. J. Appl. Phys.54 (9), 4703–4710. 10.1063/1.332803
25
Fernández-AlvarezL.VelascoV. R. (1998). Sagittal Elastic Waves in Fibonacci Superlattices. Phys. Rev. B57, 14141–14147. 10.1103/physrevb.57.14141
26
GastanierM. P.PierreC. (1997). Predicting Localization via Lyapunov Exponents Statistics. J. Sound Vib203, 151–157.
27
GaziN. A.BernhardG. (2014). Quasi-periodic Fibonacci and Periodic One-Dimensional Hypersonic Phononic Crystals of Porous Silicon: Experiment and Simulation. J. Appl. Phys.116, 094903.
28
GolubM. V.ZhangC.WangY.-S. (2012). SH-wave Propagation and Scattering in Periodically Layered Composites with a Damaged Layer. J. Sound Vibration331, 1829–1843. 10.1016/j.jsv.2011.12.008
29
GurtinM. E.WeissmüllerJ.LarchéF. (1998). A General Theory of Curved Deformable Interfaces in Solids at Equilibrium. Philosophical Mag. A78 (5), 1093–1109. 10.1080/01418619808239977
30
HashemiS. H.SamaeiA. T. (2011). Buckling Analysis of Micro/nanoscale Plates via Nonlocal Elasticity Theory. Physica E: Low-dimensional Syst. Nanostructures43, 1400–1404. 10.1016/j.physe.2011.03.012
31
HeirecheH.TounsiA.BenzairA.MaachouM.Adda BediaE. A. (2008). Sound Wave Propagation in Single-Walled Carbon Nanotubes Using Nonlocal Elasticity. Physica E: Low-dimensional Syst. Nanostructures40, 2791–2799. 10.1016/j.physe.2007.12.021
32
HepplestoneS. P.SrivastavaG. P. (2008). Hypersonic Modes in Nanophononic Semiconductors. Phys. Rev. Lett.101 (10), 105502. 10.1103/physrevlett.101.105502
33
Hu anA.Jiang Shu-shengS. S.Peng Ru-WenR. W.Zhang Chun-shengC. S.Feng DuanD. (1992). Extended One-Dimensional Fibonacci Structures. Acta Phys. Sin.41 (1), 62–68. 10.7498/aps.41.62
34
HuH.LiD.-J.ZhuJ.-L.XiongJ.-J. (2000). Size Effects on Excitons in Nano-Rings. J. Phys. Condens. Matter12, 9145–9151. 10.1088/0953-8984/12/43/304
35
HuangG. L.SunC. T. (2007). Continuum Modelling of Solids with Micro/nanostructures. Philosophical Mag.87, 3689–3707. 10.1080/14786430701400137
36
KanaN.KhamlichS.Kana KanaJ. B.MaazaM. (2013). Peculiar Surface Size-Effects in Nacl Nano-Crystals. Surf. Rev. Lett.20, 1350001. 10.1142/s0218625x13500017
37
KeL.-L.WangY.-S.WangZ.-D. (2012). Nonlinear Vibration of the Piezoelectric Nanobeams Based on the Nonlocal Theory. Compos. Structures94, 2038–2047. 10.1016/j.compstruct.2012.01.023
38
KingP. D. C.CoxT. J. (2007). Acoustic Band Gaps in Periodically and Quasiperiodically Modulated Waveguides. J. Appl. Phys.102, 014902. 10.1063/1.2749483
39
KisselG. J. (1991). Localization Factor for Multichannel Disordered Systems. Phys. Rev. A.44, 1008–1014. 10.1103/physreva.44.1008
40
KushwahaM. S.HaleviP.DobrzynskiL.Djafari-RouhaniB. (1993). Acoustic Band Structure of Periodic Elastic Composites. Phys. Rev. Lett.71, 2022–2025. 10.1103/physrevlett.71.2022
41
LuntiaovM.RogersonG. A. (2010). Long Wave Motion in Layered Elastic media. Int. J. Eng. Sci.48 (21), 1856–1871.
42
MerlinR.BajemaK.ClarkeR.JuangF.-Y.BhattacharyaP. K. (1985). Quasiperiodic GaAs-AlAs Heterostructures. Phys. Rev. Lett.55, 1768–1770. 10.1103/physrevlett.55.1768
43
MindlinR. D. (1965). Second Gradient of Strain and Surface-Tension in Linear Elasticity. Int. J. Sol. Structures1 (4), 417–438. 10.1016/0020-7683(65)90006-5
44
NguyenT. T.AbdelmoulaR.LiJ.RoussigneY.StashkevichA. (2016). Wave Propagating in Multilayers Composed of Piezo Electric and Piezo Magnetic Layers. Composites B: Eng.93, 289–301. 10.1016/j.compositesb.2016.03.019
45
NougaouiA.RouhaniB. D. (1987). Elastic Waves in Periodically Layered Infinite and Semi-infinite Anisotropic media. Surf. Sci.185, 125–153. 10.1016/s0039-6028(87)80618-0
46
NowinskiJ. L. (1984). On the Nonlocal Aspects of the Propagation of Love Waves. Int. J. Eng. Sci.22, 383–392. 10.1016/0020-7225(84)90073-9
47
NowinskiJ. L. (1991). On Wave Propagation in an Elastic Micropolar Medium with Nonlocal Cohesive Forces. J. Franklin Inst.328 (4), 459–470. 10.1016/0016-0032(91)90020-4
48
ParsonsL. C.AndrewsG. T. (2009). Observation of Hypersonic Phononic crystal Effects in Porous Silicon Superlattices. Appl. Phys. Lett.95, 241909. 10.1063/1.3275742
49
RamprasadR.ShiN. (2005). Scalability of Phononic crystal Heterostructures. Appl. Phys. Lett.87 (11), 111101. 10.1063/1.2043242
50
SesionP. D.AlbuquerqueE. L.ChesmanC.FreireV. N. (2007). Acoustic Phonon Transmission Spectra in Piezoelectric AlN/GaN Fibonacci Phononic Crystals. Eur. Phys. J. B58, 379–387. 10.1140/epjb/e2007-00249-5
51
ShaatM. (2017). A General Nonlocal Theory and its Approximations for Slowly Varying Acoustic Waves. Int. J. Mech. Sci.130, 52–63. 10.1016/j.ijmecsci.2017.05.038
52
SigalasM. M.SoukoulisC. M. (1995). Elastic-wave Propagation through Disordered And/or Absorptive Layered Systems. Phys. Rev. B Condens Matter51, 2780–2789. 10.1103/physrevb.51.2780
53
ToupinR. A. (1962). Elastic Materials with Couple-Stresses. Arch. Rational Mech. Anal.11 (1), 385–414. 10.1007/bf00253945
54
WaksmanskiN.PanE. (2017). An Analytical Three-Dimensional Solution for Free Vibration of a Magneto-Electro-Elastic Plate Considering the Nonlocal Effect. J. Intell. Mater. Syst. Structures28 (11), 1501–1513. 10.1177/1045389x16672734
55
WolfA.SwiftJ. B.SwinneyH. L.VastanoJ. A. (1985). Determining Lyapunov Exponents from a Time Series. Physica D: Nonlinear Phenomena16, 285–317. 10.1016/0167-2789(85)90011-9
56
YanD.-J.ChenA.-L.WangY.-S.ZhangC.GolubM. (2018). In-plane Elastic Wave Propagation in Nanoscale Periodic Layered Piezoelectric Structures. Int. J. Mech. Sci.142-143, 276–288. 10.1016/j.ijmecsci.2018.04.054
57
YanD.-J.ChenA.-L.WangY.-S.ZhangC. (2020). Size-effect on the Band Structures of the Transverse Elastic Wave Propagating in Nanoscale Periodic Laminates. Int. J. Mech. Sci.180, 105669. 10.1016/j.ijmecsci.2020.105669
58
YanZ.-Z.ZhangC. (2012). Band Structures and Localization Properties of Aperiodic Layered Phononic Crystals. Physica B: Condensed Matter407, 1014–1019. 10.1016/j.physb.2012.01.001
59
YanZ.-Z.ZhangC.WangY.-S. (2009). Analysis of Wave Propagation and Localization in Periodic/disordered Layered Composite Structures by a Mass-spring Model. Appl. Phys. Lett.94, 161909–161913. 10.1063/1.3119206
60
YanZ.-Z.ZhangC.WangY.-S. (2010). Wave Propagation and Localization in Randomly Disordered Layered Composites with Local Resonances. Wave Motion47, 409–420. 10.1016/j.wavemoti.2010.02.002
61
YangF.ChongA. C. M.LamD. C. C.TongP. (2002). Couple Stress Based Strain Gradient Theory for Elasticity. Int. J. Sol. Structures39 (10), 2731–2743. 10.1016/s0020-7683(02)00152-x
62
YuD.WenJ.ShenH.XiaoY.WenX. (2012). Propagation of Flexural Wave in Periodic Beam on Elastic Foundations. Phys. Lett. A376, 626–630. 10.1016/j.physleta.2011.11.056
63
ZárateJ. E.Fernández-AlvarezL.VelascoV. R. (1999). Transverse Elastic Waves in Fibonacci Superlattices. Superlattices and Microstructures25, 519–526. 10.1006/spmi.1998.0682
Summary
Keywords
nanoscale layered structures, disorder, band structures, nonlocal elastic continuum theory, aperiodicity
Citation
Yan Z and Yang X (2022) Band Structure Analysis of SH Wave Propagating in Nanoscale Layered Metamaterial Structures. Front. Mater. 9:842073. doi: 10.3389/fmats.2022.842073
Received
23 December 2021
Accepted
20 January 2022
Published
01 March 2022
Volume
9 - 2022
Edited by
Yan-Feng Wang, Tianjin University, China
Updates
Copyright
© 2022 Yan and Yang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhizhong Yan, zzyan@bit.edu.cn
This article was submitted to Metamaterials, a section of the journal Frontiers in Materials
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.