Study on Band Gap Characteristics of Boat-Shaped Phononic Crystal ()
1. Introduction
Phononic crystals (PCs) are periodic structure that can prohibit the propagation of acoustic or elastic waves with frequency in the range of phononic bandgap [1]. In recent years, theoretical and applied exploration of PCs has attracted extensive research interest due to the unique function of the band gap and has proposed a wealth of applications. Especially, it has a good application prospect in damping and noise reduction and elastic wave control, including but not limited to sound absorbers [2] [3], vibration isolators [4], acoustic cloaks [5] [6], filters [7], and waveguides [8]. It is generally believed that there are two mechanisms for generating band gaps in phononic crystals, namely Bragg scattering [9] and local resonance [10]. For the former, the interaction between unit cells or the periodicity of the structure leads to the appearance of the band gap; for the latter, the resonance properties of a single scatterer play a decisive role in the appearance of the band gap. Based on these two mechanisms, in order to obtain better bandgap properties, a host of academic research has been done on different periodic structures. Based on the Bragg scattering theory, Wang [11] studied the effect of a plate with cross cavities on the band gap. Wang [12] studied the effect of different inclusions in the same matrix on the band gap by the finite element method. Seongmin et al. [13] proposed a tapered phononic beam with a broad low-frequency band gap for flexural waves. Gao et al. [14] proposed and discussed the elastic wave attenuation of hollow metamaterial beam embedded acoustic black hole. M. Thota et al. [15] proposed and studied an origami phononic structure with tunable waveguiding. The boat-shaped phononic crystal we studied in this paper was proposed in a previous work of our group [16]. Compared with conventional Bragger scattering PCs and local resonance PCs, the advantage of Boat-shaped phononic crystal is the larger parameter spaces for design and modulation of bandgaps.
The structure of this paper is as follows: Firstly, we introduce the structural characteristics of Boat-shaped phononic crystal and investigate the band gap properties of this phononic crystal using the finite element method. Secondly, the parametric study was performed and the results of how parameters affect the bandgaps were discussed. Lastly, the conclusions are summarized.
2. Structure and Theoretical Basis
The structure of the Boat-shaped phononic crystal studied in this paper is shown in Figure 1(a). Due to the mirror symmetry of the structure along both the x-axis and the y-axis, the unit cell can be described by the following geometric parameters: The lengths of rods
and
are
= 15.4 mm and
= 15.7 mm respectively.
is the space angle between rod
and rod
,
is the angle between the plane where notes
,
and
are located and the plane where notes
,
and
are located. The values of
and
are
= 48.55˚ and
= 110.7˚ respectively; all rods have a radius R of 1 mm. The material parameters set during the numerical simulation are as follows:
= 1800 kg/m3, E = 0.86 GPa and
= 0.45 represent the density, Young’s modulus and Poisson ratio, respectively.
In a linear elastic, isotropic, homogeneous medium with infinite volume, ignoring the effect of damping, the vector form of the elastic wave equation can be written as:
(1)
In Equation (1), r is the spatial position vector, u is the mass displacement vector,
is the density of the medium,
and
are the Lame constants of the material,
denotes the characteristic circle frequency. Due to the spatial periodicity and symmetry of the phononic crystal, the eigenfrequency and eigenmode of the wave field have a certain periodicity. According to the Bloch theorem, the eigenwave field in the periodic structure can be expressed as:
(2)
In Equation (2), k denotes Bloch wave vector. The translation of the amplitude function
to the lattice vector R is also periodic, that is:
(3)
Due to the symmetry of the phononic crystal point group, the wave vector k can only take a value in the first irreducible Brillouin zone. Figure 1(b) shows the first Brillouin zone of Boat-shaped phononic crystal, and the shaded area denotes the irreducible Brillouin zone. Due to the function of setting periodic boundaries, all results in this paper are calculated by commercial code COMSOL. Figure 2 shows the band structure of Boat-shaped phononic crystal when all
Figure 1. (a) The unit cell of Boat-shaped phononic crystal.
,
,
and
are the labels at the node,
is the angle between rod
and rod
,
is the angle between the plane where notes
,
and
are located and the plane where notes
,
and
are located. (b) Brillouin zone of a rectangular lattice. Shading indicates irreducible Brillouin zone.
Figure 2. Band structure of the model under original size data. The yellow shaded area are the band gaps.
parameters are initialized. We find that the first band gap appears in the interval (1942 Hz, 2204 Hz) and the second band gap appears in the interval (11,216 Hz, 11,640 Hz). The band gap widths for low and high frequencies are 262 Hz and 424 Hz, respectively.
3. Numerical Simulations and Results
In this section, we discuss the effect of geometric parameters
,
,
and
on the low frequency first band gap and high frequency second band gap of boat-shaped phononic crystal. Noting that other geometric dimensions and material properties remain unchanged when the single parameter changes. As parameter
increase from 30˚ to 90˚, the lower bound of band gap has a significant decline while the band gap expands, as shown in Figure 3(a). Figure 3(b) shows the effect of
on band gap, the width of the band gap is basically unchanged when
is between 70˚ - 140˚, but when
increases from 140˚ to 150˚, the band gap shrinks and closes in 150˚. As geometric parameter
increase from 11 mm to 21 mm, the lower bound and upper bound of first band gap all show a downward trend and the bandwidth widen first and then narrow, and the result is shown in Figure 3(c). Figure 3(d) shows the effect of
on band gap, when
changes from 12 mm to 20 mm, the position of the band gap decreases, and the bandwidth shows the same pattern as in Figure 3(c). By analyzing Figure 3, we can see that the bandwidth is more sensitive to
, and
has the greatest effect on the bandgap position.
After analyzing the effect of parameters on the low frequency bandgap, now we study the relationship between these four parameters (
,
,
,
) and the high frequency band gap, the results are shown in Figure 4. We can find
Figure 3. Effect of geometric parameter
,
,
and
on the low frequency first band gap of Boat-shaped phononic crystal. (a), (b), (c) and (d) are change of first band gap as the respective change of
,
,
and
. The circle and square lines denote the upper and lower bounds of the low frequency bandgap, respectively.
Figure 4. Effect of geometric parameter
,
,
and
on the high frequency second band gap of Boat-shaped phononic crystal. (a), (b), (c) and (d) are change of second band gap as the respective change of
,
,
and
. The circle and square lines denote the upper and lower bounds of the high frequency bandgap, respectively.
that as
,
, and
increase within their respective ranges, the corresponding band gaps first open and then close, and bandwidths first widen and then narrow, as shown in Figure 4(a), Figure 4(b) and Figure 4(d). In contrast, as shown in Figure 4(c), as parameter
increase from 14 mm to 21 mm, the position of high frequency bandgap has a significant decrease while the bandwidth keeps increasing. The maximum of bandwidth of High frequency second bandgap occurs at
= 21 mm. Within the variation interval of
, the variation of the band gap position is also the largest. These show that both the bandgap position and the bandwidth are most sensitive to the parameter
.
4. Conclusion
In summary, we used the finite element method to study the band structure of the Boat-shaped phononic crystal and found that it has elastic wave band gaps at both high and low frequencies. In the parametric study of four key design parameters, the law of the influence of parameters on the band gap was found. For low frequency band gaps, the width of the band gap can be effectively changed by changing the size of
while
and
greatly affect the bandgap position. For high band gaps,
and
have a large effect on the band gap position. These laws will inspire us to design boat-shaped phononic crystal with more excellent or specific band gaps that can be applied in specific fields.