# Research on the Band Gap Characteristics of Two-Dimensional Phononic Crystals Microcavity with Local Resonant Structure

Authors: Mao Liu; Pei Li; Yongteng Zhong; Jiawei Xiang
Journal: Shock and Vibration
 (2015)

---

## Abstract

A new two-dimensional locally resonant phononic crystal with microcavity structure is proposed. The acoustic wave band gap characteristics of this new structure are studied using finite element method. At the same time, the corresponding displacement eigenmodes of the band edges of the lowest band gap and the transmission spectrum are calculated. The results proved that phononic crystals with microcavity structure exhibited complete band gaps in low-frequency range. The eigenfrequency of the lower edge of the first gap is lower than no microcavity structure. However, for no microcavity structure type of quadrilateral phononic crystal plate, the second band gap disappeared and the frequency range of the first band gap is relatively narrow. The main reason for appearing low-frequency band gaps is that the proposed phononic crystal introduced the local resonant microcavity structure. This study provides a good support for engineering application such as low-frequency vibration attenuation and noise control.


---

## Body

## 1. Introduction

In recent years, a growing interest has been focused on the study of the propagation of elastic waves in the periodic phononic crystals (PCs) [1–4]. PCs composed of two or more materials with different mechanical properties will generate band gaps (BGs) within which the propagation of elastic or acoustic waves is prohibited. The great attention on BGs is mainly due to their rich physical properties and potential applications in the design of acoustic devices such as noise reduction, waveguides, and acoustic filters [5, 6].In order to promote the application of phononic crystals in the fields of vibration control and noise isolation, obtaining tunable BGs with large bandwidth in the low-frequency domain is significant. In the 2DPCs, a lot of studies on the research of large band gap have been carried out to obtain PC structures with excellent BGs. For example, Lai et al. have studied the PC structure composed of the square lattice of steel cylinders in air background and found that the BGs could be tunable with various microstructures [7]. Wang et al. have introduced the narrow slit structures into the inclusions of the two-dimensional PC structure composed of steel tubes in air and obtained large BGs in audible frequency range. Wang et al. discussed the absolute band gap (ABG) structures in PCs with cross-like holes [8]. Liu et al. designed a tree-component PC structure that exhibits BGs two orders of magnitude smaller than the relevant wavelength which proposed the localized resonance (LR) mechanism [9]. As the LRPCs allowed the acquisition of low-frequency band gaps without considerably increasing the size of the structures, they encouraged more extensive application prospects and attracted more and more attention.Recently, some novel periodic structures with large band gaps are reported. Trabelsi et al. investigated the band properties of a phononic crystal composed of alternating fluid and fluid-saturated porous layers [10]. They observed low-frequency band gaps of underwater structure. Manktelow et al. studied the nonlinear periodic materials [11] and obtained large band gap shifts in one and two dimensions via optimal topologies. Salehi et al. researched the propagation of acoustic waves in the phononic crystal of 3D with rhombohedral lattice [12]. Sainidou and Stefanou analyse guided and quasiguided elastic waves in a thin glass slab coated on one side with a periodic monolayer of polymer spheres, immersed in water [13]. Oudich et al. researched the band gaps in the PC structure constructed by periodically depositing single-layer or two-layer stubs on the surface of a thin plate and proposed a “spring-mass” model to explain the opening of the LR band gaps [14]. Yu et al. studied 2DPCs with neck structures and obtained the large band gaps in the low frequency [15]. Wang et al. also explained the ternary locally resonant PCs with a comb-like coating using “spring-mass” model [8].In this paper, new band gaps structure composed of square scatter with microcavity structure embedded in a homogeneous matrix is considered. We adopted the finite element method to investigate the variable of the band gaps in low frequency. Different numbers of connection bridge plates and connection forms will be studied to confirm the locally resonant mechanism. Furthermore, microcavity structures with different geometric dimension are studied in this work. Finally, some conclusions are given.


## 2. Model and Method of the Calculation

PC structure considered here is a square lattice unit cell which contains a square inclusion and junction plates embedded in a homogeneous matrix. Figures1(a) and 1(b) illustrate the schematic view of the cross-sections of the PC structure and one of the unit cells. Figures 1A, 1B, and 1C denote the scatter, the matrix, and the microcavity, respectively. The extended direction is selected along thex- andy-axes. As shown in Figure 1(b), the microcavity structure falls in between scatter and matrix. Square scatter is connected to the matrix via 4 connection bridge plates (in dashed box shown in Figure 1(b)) which are the same materials with scatter. The length and width of the connection bridge plates are L and b, respectively. The lattice spacing of the square constant is a.Figure 1
The cross-section structure of the phononic crystal (a) and the unit cell structure of the phononic crystal (b).
(a)
(b)In the present study, to research the band gap characteristics of the proposed 2DPC structure, a lot of computes on the dispersion relations and transmission spectra are performed with the FEM [16]. For the calculations of the band structures of the proposed PCs, the considered structure refers to an infinite system. The governing field equations for elastic wave propagation in solids are given by (1)
∑
j
=
1
3
∂
∂
x
j
∑
l
=
1
3
∑
k
=
1
3
c
i
j
k
l
∂
u
k
∂
x
l
=
ρ
∂
2
u
i
∂
t
2
i
=
1,2
,
3
,where ρ is the mass density, c
i
j
k
l are the elastic constants, u
i is the displacement, t is the time, and x
j (j
=
1,2
,
3) denotes the coordinate variables x, y, and z, respectively. Because the infinite system is periodic along the x- and y-directions simultaneously, only the unit cell shown in Figure 1(b) is taken into account, based on the Bloch theorem. With the FEM, the Bloch theorem is adopted through the implementation of the Bloch periodic boundary conditions on the mentioned boundaries of the unit cell. The periodic boundary conditions on the two group boundaries of the unit cell yield (2)
U
r
+
a
=
e
i
k
·
a
U
r
,where r is located at the boundary nodes, k is wave vectors, U is the displacement at the nodes, anda is the vector that generates the point lattice associated with the phononic crystals. As the angular frequency ω is a periodical function of the wave vector, the problem can be reduced to the first Brillouin zone. The dispersion curves are eventually built by varying the wave vector k
=
(
k
x
,
k
y
) on the first Brillouin zone for a given propagation direction. The full band structure as well as the eigenmodes of the structure is then deduced using symmetries.For the analysis of the transmission spectra, a finite system must be defined. We consider here the structure being finite in thex-direction that contains N units. In the y-direction are M units; the periodic boundary conditions are still utilized to represent the infinite units. In this case, a finite array structure composed of N
×
M units is modeled for the calculations. The plane waves with single frequency, supplied by the acceleration excitation source, are incident from the left side of the finite array and propagate along the x-direction, and the corresponding transmitted acceleration value is recorded on the right side. The transmission is defined as (3)
T
L
=
20
l
o
g
10
v
o
v
i
,where v
o and v
i are the value of the transmitted and the incident displacement, respectively. Then, by varying the excitation value of the incident displacement, the transmission spectra can be obtained.


## 3. Numerical Examples and Discussion

In this section, the proposed 2DPCs with microcavity structures will be studied using the finite element software COMSOL. The detailed process can be found in [17]. The wave vectors are restricted within the x-y plane and point to the boundary ΓХMΓ of the irreducible first BZ of the square lattice structure. Based on the Bloch theorem, this is sufficient to determine the complete frequency band structures. The material parameter of two kinds of materials is shown in Table 1. The following structure parameters are used: L
=
0.001 m and b
=
0.0002 m, and the square lattice constant is a
=
0.003 m.Table 1
Materials properties.
Materials
Density (ρ)
Young’s modulus (E)
Poisson’s ratio (ν)
(Kg/m3)
Pa
A
4.08
11600
0.369
B
0.435
1180
0.368We observe that bands are contained in the frequency range of 0–500 kHz, in Figure2(a), where three complete band gaps are involved. The lowest one ranges from 14.93 to 26.92 kHz (between the second and third band), while the second and third band gaps are extended from 28.26 to 31.76 kHz (between the third and fourth bands) and 42.29 to 49.02 kHz (between the tenth and eleventh bands), respectively. In order to verify the calculations of the band gaps, the transmission spectra through the PC structure composed of a finite array of 5 × 5 units are computed. As shown in Figure 2(b), clearly existing a frequency range extending from 14.93 to 32 kHz in which the attenuation is so large, thus it can be treated as a complete band gap. The starting and terminating frequencies of the band structure coincided well with the result of band gap compute using FEM in Figure 2(a), validating the numerical results of the band structure calculations. The appearance of the two improving frequency ranges in the band structure, which correspond precisely to the third and fourth bands in the band gap, is mainly owed to the existence of the localized modes.Figure 2
The dispersion curves of the phononic crystal with microcavity structure (a) and the transmission spectra of the phononic crystal with microcavity (b).
(a)
(b)As a comparison, we also compute the gap characteristics of the traditional PC composed of scatter inclusions embedded in the epoxy matrix with full contact, namely, no microcavity structure. Obviously, it is a classic 2DPC structure based on the Bragg scattering mechanism. During the computations, the lattice constant and filling ratio are selected to be the same as those used in Figure2. The results are shown in Figure 3. It can be revealed that there are few bands in the frequency range of 0–500 kHz. The only band gap locates upon the second band, with starting frequency 20.34 kHz. At the same time, the transmission spectra of the finite array of 5 × 5 units are calculated to confirm the analysis of the band structure. The results are shown in Figure 5; the attenuation frequency range of the acceleration responses coincides well with the frequency range that parallels the band gap. Comparing the structures with microcavity, we can discover that the frequency range of the band structure in PC with microcavity structures is lower at the lower edge of band gap than that of the PCs having no microcavity structures; at the same time, the width of the gap is extremely large compared to classic PC structures and meanwhile the whole weight of the microcavity structures is lighter via some material will removed. All these excellent characteristics will promote the applications of these new microcavity PC structures in the low-frequency noise control in engineering field and the design of the low-frequency reductions.Figure 3
The dispersion curves of the phononic crystal with no microcavity structure (a) and the transmission spectra of the phononic crystal with no microcavity (b).
(a)
(b)In order to confirm the physical mechanism for the occurrence of the low-frequency band structures in the proposed PC structure, we compute the eigenmode shapes and displacement vector fields at the edges of the first complete band gap. The results are shown in Figure4 and the behavior is discussed in the following. The color maps of Figure 4 are the magnitude of the total displacement vector field. The modes on the lower and upper edges of microcavity structure are discussed first. It can be obviously observed that the vibration of mode A, which corresponds to the lower edge of the first band structure, is mainly concentrated at the square scatter and four connection bridge plates, while the matrix nearly remains motionless. The connection bridge plates extrude the matrix with the result that localized large deformation. The vibration process can be seen as a mass-spring system. The square scatter works as a resonator and the connection bridge plates play the role of the spring. With different mode A, the vibration of mode B is mainly the torsion of the matrix and meanwhile the out-of-phase vibration of the matrix. The square scatter and the connection bridge plates remain static. Because the vibration direction of the matrix is outward, the square scatter and the connection bridge plates remain localized. Then, we investigated the modes on the lower and upper edges of classic PC structure. From Figures 4(c) and 4(d), corresponding to the structure vibration of modes C and D are mainly extrudes and compresses the matrix; that is to say, there have been no localized conditions.Figure 4
Eigenmodes shapes and displacement vector fields of the modes marked in Figures2 and 3 correspond to points A, B, C, and D, respectively.
(a)
(b)
(c)
(d)Figure 5
The dispersion curves of the phononic crystal with three connection plates (a). The dispersion curves of the phononic crystal with two opposite connection plates (b). The dispersion curves of the phononic crystal with two adjacent connection plates (c). The dispersion curves of the phononic crystal with only one connection plate.
(a)
(b)
(c)
(d)From the analysis above, we can conclude that the occurrence of the large band gaps in low frequency is mainly attributed to the existence of the localized resonance modes resulting from the introduction of the four connection bridge plates formed microcavity structures. Since the connection bridge plates play the role of springs linking the square scatter and the matrix, their shape features may influence the band structure remarkably. In the following, we will study the effects of the geometrical parameters on the typical large low-frequency band gaps of the proposed PC structure.Furthermore, we also investigated the band gaps of microcavity structures with 1, 2, and 3 connection bridge plates to the matrix, respectively. The calculated dispersion curves are shown in Figure5. The variation of the number of connection bridge plates, however, does not significantly influence the band structures, and they show for all four scatter shapes very similar features. We can observe clearly existing incomplete band gaps X-M direction of the irreducible first Brillouin zone in Figures 5(a) and 5(b). Because of these two connection shapes, a connection plate plays a role of the spring extruding the resonator and another connection plate works as a spring to compress the resonator. So the vibration process can be seen as two mass-spring systems with out-of-phase vibration which leads to attenuation of frequency in X-M direction. In addition, because the direction of vibration on connection plates is not consistent in Figures 5(c) and 5(d), band gaps of these structures do not exist. This is slightly surprising and indicated that the connection form provides a more important design factor.


### 3.1. Effect of the Length of Connection Plates on the Low-Frequency Band Gap

The objective of the following is to investigate the effects of geometric dimension on the low-frequency band gap with microcavity structure. First, we investigate the influence of the length of the connection bridge plates on the band gap. Keeping the lattice constanta
=
3 mm, the thickness of the matrix is 0.3 mm and the width of connection bridge plate b
=
1 mm. The results are shown in Figure 6. We can observe that the lower edges of the first band gap move to a higher frequency range with the increase of the length of the connection bridge plates L from 0.2 to 1.8 mm, while the variable is minor. At the same time, the upper edges of the first band gap are larger and the increase of the frequency of the upper edge is faster than that of the lower edge. In fact, the length of the connection plates influences not only the frequency range of the band gap but also the value of the attenuation in the gap.Figure 6
The dispersion curves of the phononic crystals with different length of connection plates structures.


### 3.2. Effect of the Width of Connection Plates on the Low-Frequency Band Gap

Second, we investigate the influence of the width of connection plates on the band structure. Figure7 displays the evolution of the band gap as a function of the width b with a
=
3 mm and L
=
1.4 mm and the thickness of the matrix is 0.3 mm, the upper edge (black dashed line with square points) and lower edge (black solid line with circle points). It can be observed that, with the width of the connection plates increasing from 1.4 to 1.8 mm, the upper edge of the band gap rises rapidly from 305.4 to 318.4 kHz, while the lower edge of the band gap declines slightly. We can see that the width of the connection plates mainly affects the upper edge of the band gap and has limited effect on the lower edge. For the upper edge, with the increase of the width, the stiffness of the connection plates in the mass-spring system increases rapidly, quicker than that of the mass of the resonator, resulting in the resonance frequency rising accordingly. For the lower edge, whose relevant mode is concentrated on the matrix, the increase of the width has little influence on relevant mode.Figure 7
The evolution of the band gap as a function of the width of the connection plateb.


### 3.3. Effect of the Density of Connection Plates on the Low-Frequency Band Gap

Finally, we also research the effect of the density of the connection plates on the band gaps. Figure8 shows the evolution of the band gap as a function of the density with four connection plates; setting a
=
3 mm, L
=
10 mm, and b
=
2 mm, the red, black, and blue rectangular boxes denote the first, second, and third band gap, respectively. It can be seen that, with the density of the connection plates increasing from 800 to 11600 kg/m3, the widths of the first and the second band gap vary slightly, while the width of the third band gap rises rapidly. For the first and second band gap, with the increase of the density of the connection plate, the stiffness of the connection plates in the mass-spring system increases rapidly, while the mass of the resonator has no change, resulting in the resonance frequency minor declining accordingly.Figure 8
The evolution of the band gap as a function of the density of the connection plate.


## 3.1. Effect of the Length of Connection Plates on the Low-Frequency Band Gap

The objective of the following is to investigate the effects of geometric dimension on the low-frequency band gap with microcavity structure. First, we investigate the influence of the length of the connection bridge plates on the band gap. Keeping the lattice constanta
=
3 mm, the thickness of the matrix is 0.3 mm and the width of connection bridge plate b
=
1 mm. The results are shown in Figure 6. We can observe that the lower edges of the first band gap move to a higher frequency range with the increase of the length of the connection bridge plates L from 0.2 to 1.8 mm, while the variable is minor. At the same time, the upper edges of the first band gap are larger and the increase of the frequency of the upper edge is faster than that of the lower edge. In fact, the length of the connection plates influences not only the frequency range of the band gap but also the value of the attenuation in the gap.Figure 6
The dispersion curves of the phononic crystals with different length of connection plates structures.


## 3.2. Effect of the Width of Connection Plates on the Low-Frequency Band Gap

Second, we investigate the influence of the width of connection plates on the band structure. Figure7 displays the evolution of the band gap as a function of the width b with a
=
3 mm and L
=
1.4 mm and the thickness of the matrix is 0.3 mm, the upper edge (black dashed line with square points) and lower edge (black solid line with circle points). It can be observed that, with the width of the connection plates increasing from 1.4 to 1.8 mm, the upper edge of the band gap rises rapidly from 305.4 to 318.4 kHz, while the lower edge of the band gap declines slightly. We can see that the width of the connection plates mainly affects the upper edge of the band gap and has limited effect on the lower edge. For the upper edge, with the increase of the width, the stiffness of the connection plates in the mass-spring system increases rapidly, quicker than that of the mass of the resonator, resulting in the resonance frequency rising accordingly. For the lower edge, whose relevant mode is concentrated on the matrix, the increase of the width has little influence on relevant mode.Figure 7
The evolution of the band gap as a function of the width of the connection plateb.


## 3.3. Effect of the Density of Connection Plates on the Low-Frequency Band Gap

Finally, we also research the effect of the density of the connection plates on the band gaps. Figure8 shows the evolution of the band gap as a function of the density with four connection plates; setting a
=
3 mm, L
=
10 mm, and b
=
2 mm, the red, black, and blue rectangular boxes denote the first, second, and third band gap, respectively. It can be seen that, with the density of the connection plates increasing from 800 to 11600 kg/m3, the widths of the first and the second band gap vary slightly, while the width of the third band gap rises rapidly. For the first and second band gap, with the increase of the density of the connection plate, the stiffness of the connection plates in the mass-spring system increases rapidly, while the mass of the resonator has no change, resulting in the resonance frequency minor declining accordingly.Figure 8
The evolution of the band gap as a function of the density of the connection plate.


## 4. Conclusions

In this paper, the band structures of two-dimensional PCs with the microcavity structure are investigated using FEM. From these PC structures, a complete band gap can be obtained. The band gap is notably lowered, as expected. The microcavity structure formed with four connection plates bridge to the resonator possess many locally resonators lead to a low-frequency band gap. The displacement fields of eigenmodes are studied to reveal the physical mechanism for the existence of material’s periodic arrangement on the low-frequency band gaps. It was shown that the influence of the geometric dimension and connection plate adjacent shape is much higher than the influence of the materials properties. This provides a more important design factor in engineering application.



---

*Source: 239832-2015-10-18.xml*