# An Enhanced Plane Wave Expansion Method to Solve Piezoelectric Phononic Crystal with Resonant Shunting Circuits

Authors: Ziyang Lian; Shan Jiang; Hongping Hu; Longxiang Dai; Xuedong Chen; Wei Jiang
Journal: Shock and Vibration
 (2016)

---

## Abstract

An enhanced plane wave expansion (PWE) method is proposed to solve piezoelectric phononic crystal (PPC) connected with resonant shunting circuits (PPC-C), which is named as PWE-PPC-C. The resonant shunting circuits can not only bring about the locally resonant (LR) band gap for the PPC-C but also conveniently tune frequency and bandwidth of band gaps through adjusting circuit parameters. However, thus far, more than one-dimensional PPC-C has been studied just by Finite Element method. Compared with other methods, the PWE has great advantages in solving more than one-dimensional PC as well as various lattice types. Nevertheless, the conventional PWE cannot accurately solve coupling between the structure and resonant shunting circuits of the PPC-C since only taking one-way coupling from displacements to electrical parameters into consideration. A two-dimensional PPC-C model of orthorhombic lattice is established to demonstrate the whole solving process of PWE-PPC-C. The PWE-PPC-C method is validated by Transfer Matrix method as well as Finite Element method. The dependence of band gaps on circuit parameters has been investigated in detail by PWE-PPC-C. Its advantage in solving various lattice types is further illustrated by calculating the PPC-C of triangular and hexagonal lattices, respectively.


---

## Body

## 1. Introduction

The propagation of elastic or acoustic waves in phononic crystals (PCs) has attracted growing attention during the past two decades [1, 2], both because of its amazing physical properties and because of its potential applications. One of the most attractive properties of PCs is the band gaps, in which elastic or acoustic waves are attenuated significantly. According to the generation mechanisms of Bragg scattering and locally resonant (LR), the band gaps can be divided into Bragg band gap and LR band gap [3].The PC has the potential to be used as a vibration isolator because of the band gap, but one of main difficulties of isolating mechanical vibration lies in how to attenuate low-frequency vibration. Compared to the Bragg band gap, the LR band gap has a lower frequency and can be generated by a smaller PC, which enables the PC to isolate low-frequency vibration for the high-precision machine [3–5]. Yet a unit cell of conventional LR PCs commonly consists of a heavy core and a soft coat, as a mass-spring oscillator, to obtain a low resonance frequency. Therefore, a lower frequency LR band gap implies a heavier PC. By replacing the heavy mass-spring oscillator, a light inductor-capacitor oscillator can be constituted by the equivalent capacitance of piezoelectric phononic crystal (PPC) and the inductance of the resonant shunting circuit. The beam-like piezoelectric phononic crystal connected with resonant shunting circuits (PPC-C) was studied to obtain LR band gaps for controlling the propagation of vibration [6–10]. Furthermore, the PPC-C has another significant advantage that its band gaps can be tuned just by circuit parameters, but reconfiguring the structure is not needed [11–16]. This property will further enhance the application flexibility of the PPC-C.The PPCs have been researched by Transfer Matrix method [10], plane wave expansion (PWE) method [17, 18], Finite Difference Time Domain method [19], and Finite Element method [20, 21]. Nevertheless, among these methods, only Transfer Matrix method [6–9, 11–13] and Finite Element methods [15, 16, 22] have been used to study the PPC-C. Transfer Matrix method is just applicable to one-dimensional PPC-C, such as beam [7] and circle plate [23]. Up to now, two-dimensional plate-like PPC-C have been researched merely by Finite Element method [15, 16, 22]. Although the utilization of Finite Element method becomes popular nowadays, it is useful to have another method to handle this kind of problem.Two-dimensional and three-dimensional PPCs were investigated by the methods of Finite Difference Time Domain [19] and conventional PWE [17, 18]. However, without resonant shunting circuits are connected with piezoelectric patches of the PPC. Under this circumstance, we only need to consider two kinds of electric boundary conditions for the piezoelectric patches, that is, short and open circuits. When electrodes of one piezoelectric patch are shorted, electric field is equal to zero in a certain direction, but when these electrodes are open, normal component of electric displacement is equal to zero. Electric field or electric displacement can be expressed directly by displacements. Hence the band structure can be obtained by solving governing equation in terms of displacements. Yet for the PPC with circuits (PPC-C), electrical parameters, such as electric potential and electric displacement, are not contained in the equations of motion. Only one-way coupling, that is, from displacements to electrical parameters, can be taken into account. Firstly, displacements are always obtained from the equations of motion and then get electrical parameters according to circuit equation. Hence the influence of the circuit on structure movement is not taken into consideration. Thus the PPC-C cannot be solved accurately by conventional PWE.The PWE method is one of the most extensively used methods to calculate band structure because of its convenience [24]. Kushwaha et al. [25] firstly obtain the band structure of a PC by PWE. PWE method includes three steps as follows. Firstly, all parameters such as modulus, density, and Poisson’s ratio are expanded as the Fourier series in the reciprocal space. Secondly, equations of motion are transformed into the standard eigenvalue problem based on the Bloch theorem. Finally, band structure of the dispersion relations between the frequency and wave vector is obtained by solving the eigenvalue equation. It is convenient that wave vector just needs to traverse along the boundary of the irreducible region of the first Brillouin zone. Compared with Finite Element method, the PWE shows great advantages in dealing with scatterers of different geometries and different arrangements, because it does not require meshing and reconstructing the Finite Element matrix. But the conventional PWE has a disadvantage of the slow convergence rate, especially for systems of either very high or very low filling ratios or of large elastic mismatch. The fictitious band gaps always appear as redundant lines in the band structure. Nevertheless, an improved PWE has a good convergence and can provide much more accurate numerical results [26, 27].The main aim of the paper is to propose an enhanced PWE method to solve the PPC-C, which is abbreviated as PWE-PPC-C. A two-dimensional PPC-C model is established to demonstrate solving process of the PWE-PPC-C. The model is composed of two portions: elastic plate and elastic-piezoelectric composite plate. In general, the conventional PWE does not directly utilize the continuity conditions of forces and bending moment between two portions of the PPC. Moreover, the voltage appears in the expression of bending moment of elastic-piezoelectric composite plate but does not exist in that of elastic plate. To deal with this voltage mismatch, the voltage is included in equations of motion of these two portions of the PPC to ensure the continuity of bending moment at the interface between these two portions. And the voltage is introduced by multiplying an infinitely small quantity for the elastic plat. Then band structure is obtained by solving the equations of motion and circuit equation. The paper is organized as follows: Section2 presents governing equations of a PPC-C. The equations of motion are expanded by PWE-PPC-C and Bloch theorem in Section 3. Section 4 describes two approaches to solve the PPC with two kinds of common shunting circuits, respectively. The validity and correctness of the PWE-PPC-C are verified by comparisons with Transfer Matrix method and Finite Element method in Section 5. Finally, some conclusions are given in Section 6.


## 2. Governing Equations of a Two-Dimensional PPC-C

In order to demonstrate solving process of the PWE-PPC-C, we propose a two-dimensional PPC-C which consists of an isotropic elastic plate and rectangular piezoelectric patches. These piezoelectric patches are bonded on the upper and lower surfaces of the plate and are arrayed as orthorhombic system. Its lattice and arrangement are shown in Figures1 and 2(a). With The same polarization direction, each pair of parallel piezoelectric patches is connected by a resonant shunting circuit, which can be simplified by a complex impedance Z. Each unit cell of the PPC can be divided into two portions: the elastic plate as a matrix and the elastic-piezoelectric composite plate as a scatterer. Subscripts (1,2,3) of physical quantities correspond to coordinate system (x,y,z).Figure 1
Unit cell of the two-dimensional PPC-C of orthogonal lattice: (a) top view and (b) side view.
(a)
(b)Figure 2
Wigner-Seitz unit cell and its first Brillouin zone: (a) orthogonal lattice, (b) triangular lattice, and (c) hexagonal lattice.
(a)
(b)
(c)Based on the classical theory of plate [28], equation of motion of a thin plate with flexural vibration is given by (1)∂Q31∂x+∂Q32∂y=ρ0h0u¨3,where Q3ii=1,2, u3, ρ0, and h0 represent shear force, deflection, density, and thickness of the plate, respectively. A superimposed dot denotes the derivative with respect to time. The shear force-bending moment Miji,j=1,2 relations are(2)Q31=∂M11∂x+∂M21∂y,Q32=∂M12∂x+∂M22∂y.According to Kirchhoff assumption of thin plate, the strain-displacement relations with a compact matrix notation are(3)S1=-z∂2u3∂x2,S2=-z∂2u3∂y2,S6=2S12=-2z∂2u3∂x∂y.The governing equations of two portions and circuits of the PPC-C are presented, respectively, as follows.


### 2.1. Elastic Plate

For the isotropic elastic material, nontrivial stresses are(4)T1=λ+2GS1+λS2+λS3,T2=λS1+λ+2GS2+λS3,T3=λS1+λS2+λ+2GS3,T6=T12=GS6,where λ=Y1μ1/1+μ11-2μ1 and G=Y1/21+μ1. Y1 and μ1 are Young’s modulus and Poisson’s ratio of the isotropic elastic material, respectively. Stress relaxation of the thin shell, that is, T3=0, is performed on (4) and then yields(5)T1=Y-1S1+μ1S2,T2=Y-1μ1S1+S2,T6=T12=GS6,where Y-1=Y1/1-μ12. The bending moments can be written as(6)M11=∫-h/2h/2T1zdz=-De∂2u3∂x2-υ1∂2u3∂y2,M22=∫-h/2h/2T2zdz=-υ1∂2u3∂x2-De∂2u3∂y2,M12=M21=∫-h/2h/2T12zdz=-γ1∂2u3∂x∂y,where flexural stiffness of the elastic plate De=Y-1h3/12, υ1=Deμ1, and γ1=De1-μ1. Substituting (6) into (1)–(3), equation of motion (1) can be expressed in terms of displacement as(7)-η1∂2u3∂t2=∂2∂x2De∂2u3∂x2+υ1∂2u3∂y2+∂2∂y2De∂2u3∂y2+υ1∂2u3∂x2+2∂2∂x∂yγ1∂2u3∂x∂y,where η1=ρ1h is the area density of the elastic plate with volume density ρ1 and thickness h.


### 2.2. Elastic-Piezoelectric Composite Plate

The constitutive equations of the piezoelectric patches are [29](8)S1=s11T1+s12T2+d31E3,S2=s12T1+s11T2+d31E3,S6=s66T6,D3=d31T1+d31T2+ε33E3,where D3 is electric displacement in z direction; s11, s12, and s66=2(s11-s12) are compliance constants; and d31 and ε33 are piezoelectric and dielectric constants, respectively.With components in other directions being ignored, the electric field inz direction is approximately equal to [30](9)E3=-Vhp,h2≤z≤h2+hp,Vhp,-h2+hp≤z≤-h2,where V is the voltage on the two outer electrodes of the elastic-piezoelectric composite plate and hp is the thickness of the piezoelectric patches. Using (8) yields (10)T1=Y-2S1+μ2S2-d311+μ2E3,T2=Y-2μ2S1+S2-d311+μ2E3,T6=Y221+μ2S6,D3=d31Y21-μ2S1+S2+ε-33E3,where Y2=1/s11 and μ2=-s12/s11 are Young’s modulus and Poisson’s ratio of the piezoelectric material under condition of constant electric field. Y-2=Y2/1-μ22 and ε-33=ε331-k312 with k312=2d312Y2/1-μ2ε33 are equivalent Young’s modulus and equivalent two-dimensional dielectric permittivity. From (3), (9), and (10), the bending moments of the composite plate are given by(11)M11=∫-h/2+hph/2+hpT1zdz=-Ds∂2u3∂x2-υ2∂2u3∂y2+τ2V,M22=∫-h/2+hph/2+hpT2zdz=-υ2∂2u3∂x2-Ds∂2u3∂y2+τ2V,M12=M21=∫-h/2+hph/2+hpT12zdz=-γ2∂2u3∂x∂y,where flexural stiffness of the composite plate Ds=De+Dp, υ2=μ1De+μ2Dp, τ2=Y-2d31(1+μ2)(hp+h), γ2=(1-μ1)De+(1-μ2)Dp, and Dp=2Y-2(h/2+hp)3-(h/2)3/3. Substituting (11) into (2), equation of motion (1) can be rewritten as(12)-η2∂2u3∂t2=∂2∂x2Ds∂2u3∂x2+υ2∂2u3∂y2-τ2V+∂2∂y2Ds∂2u3∂y2+υ2∂2u3∂x2-τ2V+2∂2∂x∂yγ2∂2u3∂x∂y,where η2=ρ1h+2ρ2hp is the area density of the composite plate. ρ2 denotes volume density of piezoelectric patch. Although the voltage is a constant on electrodes of a pair of piezoelectric patches, the voltage is retained in the differential equation of motion to ensure that the bending moments meet the continuity condition between two portions of the unit cell.


### 2.3. Resonant Shunting Circuits

The electric currentI flowing out from one of two outer electrodes of the elastic-piezoelectric composite plate is(13)I=-Q˙e=ApD3h2+hps,where Laplace operator s=iω with imaginary unit i and angular frequency ω. Ap=l1l2 is the area of the electrode. l1 and l2 are length and width of a piezoelectric patch, respectively. The total current flowing to the impendence Z is twice the current output from one electrode since two piezoelectric patches in one unit cell are connected in parallel. Therefore the voltage and current satisfy Ohm’s law(14)2I=VZ.Using (13) and (14) yields the voltage(15)V=2ZApD3h2+hps.


## 2.1. Elastic Plate

For the isotropic elastic material, nontrivial stresses are(4)T1=λ+2GS1+λS2+λS3,T2=λS1+λ+2GS2+λS3,T3=λS1+λS2+λ+2GS3,T6=T12=GS6,where λ=Y1μ1/1+μ11-2μ1 and G=Y1/21+μ1. Y1 and μ1 are Young’s modulus and Poisson’s ratio of the isotropic elastic material, respectively. Stress relaxation of the thin shell, that is, T3=0, is performed on (4) and then yields(5)T1=Y-1S1+μ1S2,T2=Y-1μ1S1+S2,T6=T12=GS6,where Y-1=Y1/1-μ12. The bending moments can be written as(6)M11=∫-h/2h/2T1zdz=-De∂2u3∂x2-υ1∂2u3∂y2,M22=∫-h/2h/2T2zdz=-υ1∂2u3∂x2-De∂2u3∂y2,M12=M21=∫-h/2h/2T12zdz=-γ1∂2u3∂x∂y,where flexural stiffness of the elastic plate De=Y-1h3/12, υ1=Deμ1, and γ1=De1-μ1. Substituting (6) into (1)–(3), equation of motion (1) can be expressed in terms of displacement as(7)-η1∂2u3∂t2=∂2∂x2De∂2u3∂x2+υ1∂2u3∂y2+∂2∂y2De∂2u3∂y2+υ1∂2u3∂x2+2∂2∂x∂yγ1∂2u3∂x∂y,where η1=ρ1h is the area density of the elastic plate with volume density ρ1 and thickness h.


## 2.2. Elastic-Piezoelectric Composite Plate

The constitutive equations of the piezoelectric patches are [29](8)S1=s11T1+s12T2+d31E3,S2=s12T1+s11T2+d31E3,S6=s66T6,D3=d31T1+d31T2+ε33E3,where D3 is electric displacement in z direction; s11, s12, and s66=2(s11-s12) are compliance constants; and d31 and ε33 are piezoelectric and dielectric constants, respectively.With components in other directions being ignored, the electric field inz direction is approximately equal to [30](9)E3=-Vhp,h2≤z≤h2+hp,Vhp,-h2+hp≤z≤-h2,where V is the voltage on the two outer electrodes of the elastic-piezoelectric composite plate and hp is the thickness of the piezoelectric patches. Using (8) yields (10)T1=Y-2S1+μ2S2-d311+μ2E3,T2=Y-2μ2S1+S2-d311+μ2E3,T6=Y221+μ2S6,D3=d31Y21-μ2S1+S2+ε-33E3,where Y2=1/s11 and μ2=-s12/s11 are Young’s modulus and Poisson’s ratio of the piezoelectric material under condition of constant electric field. Y-2=Y2/1-μ22 and ε-33=ε331-k312 with k312=2d312Y2/1-μ2ε33 are equivalent Young’s modulus and equivalent two-dimensional dielectric permittivity. From (3), (9), and (10), the bending moments of the composite plate are given by(11)M11=∫-h/2+hph/2+hpT1zdz=-Ds∂2u3∂x2-υ2∂2u3∂y2+τ2V,M22=∫-h/2+hph/2+hpT2zdz=-υ2∂2u3∂x2-Ds∂2u3∂y2+τ2V,M12=M21=∫-h/2+hph/2+hpT12zdz=-γ2∂2u3∂x∂y,where flexural stiffness of the composite plate Ds=De+Dp, υ2=μ1De+μ2Dp, τ2=Y-2d31(1+μ2)(hp+h), γ2=(1-μ1)De+(1-μ2)Dp, and Dp=2Y-2(h/2+hp)3-(h/2)3/3. Substituting (11) into (2), equation of motion (1) can be rewritten as(12)-η2∂2u3∂t2=∂2∂x2Ds∂2u3∂x2+υ2∂2u3∂y2-τ2V+∂2∂y2Ds∂2u3∂y2+υ2∂2u3∂x2-τ2V+2∂2∂x∂yγ2∂2u3∂x∂y,where η2=ρ1h+2ρ2hp is the area density of the composite plate. ρ2 denotes volume density of piezoelectric patch. Although the voltage is a constant on electrodes of a pair of piezoelectric patches, the voltage is retained in the differential equation of motion to ensure that the bending moments meet the continuity condition between two portions of the unit cell.


## 2.3. Resonant Shunting Circuits

The electric currentI flowing out from one of two outer electrodes of the elastic-piezoelectric composite plate is(13)I=-Q˙e=ApD3h2+hps,where Laplace operator s=iω with imaginary unit i and angular frequency ω. Ap=l1l2 is the area of the electrode. l1 and l2 are length and width of a piezoelectric patch, respectively. The total current flowing to the impendence Z is twice the current output from one electrode since two piezoelectric patches in one unit cell are connected in parallel. Therefore the voltage and current satisfy Ohm’s law(14)2I=VZ.Using (13) and (14) yields the voltage(15)V=2ZApD3h2+hps.


## 3. Expansion by PWE-PPC-C Method and Bloch Theorem

The voltage multiplying an infinitesimal quantityτ1 is introduced into equation of motion of elastic plate in order that equations of motion have a unified form for two portions of the PPC. From (7) and (12), the unified form of equations of motion is (16)-η∂2u3∂t2=∂2∂x2D∂2u3∂x2+υ∂2u3∂y2-τV+∂2∂y2D∂2u3∂y2+υ∂2u3∂x2-τV+2∂2∂x∂yγ∂2u3∂x∂y,where η, D, υ, γ, and τ represent η1, De, υ1, γ1, and τ1 for elastic plate and η2, Ds, υ2, γ2, and τ2 for composite plate, respectively.Due to the periodicity,η, D, υ, γ, and τ can be expanded as Fourier series:(17)ηr=∑GηGe-iG·r,Dr=∑GDGe-iG·r,υr=∑GυGe-iG·r,γr=∑GγGe-iG·r,τr=∑GτGe-iG·r,where r=n1a1+n2a2 denotes two-dimensional lattice vector in primitive space, in which a1=a1i and a2=a2j are basis vectors. a1 and a2 are lattice constants of orthorhombic lattice along x-axis and y-axis, respectively. n1 and n2 are both integers. G=n1b1+n2b2 are the two-dimensional reciprocal lattice vectors. b1 and b2 are basis vectors of reciprocal space along x-axial and y-axial directions, respectively, which can be written as b1=2πi/a1 and b2=2πj/a2.Letg(r)=∑Gg(G)e-iG·r, where g can be any of the parameters η, D, υ, γ, and τ. As a Fourier coefficient, g(G) is defined by(18)gG=1A∬AgreiG·rd2r,where A=a1a2 denotes the area of a unit cell. The integration gives(19)gG=g2f+g11-f,G=0,g2-g1PG=ΔgPG,G≠0,where g1 and g2 represent any one kind of parameters, among η, D, υ, γ, and τ, of elastic and composite portions, respectively. f=Ap/A is defined as filling ratio of the PPC. P(G)=(1/A)∬ApeiG·rd2r is called structure factor, which depends only on the geometry of the piezoelectric patches but has nothing to do with its arrangement. For the elastic-piezoelectric composite plate, the structure factor is(20)PG=fpGx,l1pGy,l2,where subscripts x and y denote components of the vector on the x-axis and y-axis, respectively. Function p(a,b) is defined as(21)pa,b=1,a=0,sinab/2ab/2,a≠0.To satisfy the Bloch theorem, solutions of (16), the displacement and the voltage, are expanded as follows:(22)ur,t=eiωt∑GuKGe-iK+G·r,Vr,t=eiωt∑Gu4Ge-iK+G·r,where K is Bloch wave vector. Substituting (17) and (22) into (16) gives(23)ω2∑G′ηG-G′uKG′=∑G′DG-G′K+G′x2K+Gx2uKG′+∑G′υG-G′K+G′y2K+Gx2uKG′+∑G′τG-G′K+Gx2u4G′+∑G′DG-G′K+G′y2K+Gy2uKG′+∑G′υG-G′K+G′x2K+Gy2uKG′+∑G′τG-G′K+Gy2u4G′+2∑G′γG-G′K+G′xK+G′yK+GxK+GyuKG′.


## 4. Two Kinds of Resonant Shunting Circuits

### 4.1. Pure Inductive Circuit

The complex impedanceZ(ω)=iωL if each resonant shunting circuit only consists of a single inductor L. Then (10) and (15) yield(24)χ1+ω2χ2V=-ω2χ3∂2u3∂x2+∂2u3∂y2,where χ1=1, χ2=-LCp, and χ3=-LApd31Y2(h+2hp)/(1-μ2) with equivalent static capacitance Cp=2ε-33Ap/hp of the two parallel piezoelectric patches. Hence estimated oscillation frequency fr=1/2πLCp for the pure inductive resonant shunting circuit. Expanding χjj=1,2,3 as (17) and substituting (22) into (24) yield(25)∑G′χ1G-G′+ω2∑G′χ2G-G′u4G′=ω2∑G′χ3G-G′K+G′x2+K+G′y2uKG′.Then (23) and (25) can be written as matrix form(26)HuKG′+Eu4G′=ω2FuKG′,P+ω2Qu4G′=ω2WuKG′,where H, E, F, P, Q, and W are matrices of (n1+n2+1)×(n1+n2+1) and (27)H=∑G′DG-G′K+G′x2K+Gx2+K+G′y2K→+Gy2+∑G′υG-G′K+G′y2K+Gx2+K+G′x2K+Gy2+2∑G′γG-G′K+G′xK+G′yK+GxK+Gy,E=∑G′-τG-G′K+Gx2+K+Gy2,F=∑G′εG-G′,P=∑G′χ1G-G′,Q=∑G′χ2G-G′,W=∑G′χ3G-G′K+G′x2+K+G′y2.In order to facilitate solving eigenvalues, (26) is further simplified into standard form (28)F0W-Q-1HE0PuKG′u4G′=ω2uKG′u4G′,where 0 represents zero matrix with same order as H.


### 4.2. Resistive-Inductive Circuit

Complex impedanceZ(ω)=RL+iωL if each resonant shunting circuit is composed of a resistor and an inductor in series. For the resistive-inductive circuit, when RL<2L/C, the oscillation frequency can be estimated as (29)fr=12π1LCp-RL24L2.Equations (10) and (15) give(30)χ1+ωχ4+ω2χ2V=-ωχ5+ω2χ3∂2u3∂x2+∂2u3∂y2,where χ4=iRLCp and χ5=iRLApd31Y2(h+2hp)/(1-μ2). Expanding χjj=1,2,3,4,5 as (17) and substituting (22) into (30) yield(31)∑G′χ1G-G′+ω∑G′χ4G-G′+ω2∑G′χ2G-G′u4G′=ω∑G′χ5G-G′+ω2∑G′χ3G-G′K+G′x2+K+G′y2uKG′.Equations (23) and (31) may be expressed as matrix form (32)HuKG′+Eu4G′=ω2FuKG′,P+ωS+ω2Qu4G′=ωΤ+ω2WuKG′,where S=∑G′χ4(G-G′) and T=∑G′χ5(G-G′). Equation (32) can also be written as standard form [31](33)00F000W-QI0000I00-1HE000P-TS00I0000IZ1Z2Z3Z4=ωZ1Z2Z3Z4,where Z1=uK(G′), Z2=u4(G′), Z3=ωZ1, and Z4=ωZ2.In order to improve the convergence rate and the accuracy of the numerical calculation, technique of the improved PWE [26] is adopted by PWE-PPC-C. Namely, a parameter matrix g is replaced by 1/g-1, where the index −1 means inverse matrix. All parameters η, D, υ, γ, and τ are converted in the same way. It can be noted that the parameter τ1 cannot be equal to zero. The band structure is obtained when K has been traversed along the boundary of the irreducible region of the first Brillouin zone as shown in Figure 2(a).


## 4.1. Pure Inductive Circuit

The complex impedanceZ(ω)=iωL if each resonant shunting circuit only consists of a single inductor L. Then (10) and (15) yield(24)χ1+ω2χ2V=-ω2χ3∂2u3∂x2+∂2u3∂y2,where χ1=1, χ2=-LCp, and χ3=-LApd31Y2(h+2hp)/(1-μ2) with equivalent static capacitance Cp=2ε-33Ap/hp of the two parallel piezoelectric patches. Hence estimated oscillation frequency fr=1/2πLCp for the pure inductive resonant shunting circuit. Expanding χjj=1,2,3 as (17) and substituting (22) into (24) yield(25)∑G′χ1G-G′+ω2∑G′χ2G-G′u4G′=ω2∑G′χ3G-G′K+G′x2+K+G′y2uKG′.Then (23) and (25) can be written as matrix form(26)HuKG′+Eu4G′=ω2FuKG′,P+ω2Qu4G′=ω2WuKG′,where H, E, F, P, Q, and W are matrices of (n1+n2+1)×(n1+n2+1) and (27)H=∑G′DG-G′K+G′x2K+Gx2+K+G′y2K→+Gy2+∑G′υG-G′K+G′y2K+Gx2+K+G′x2K+Gy2+2∑G′γG-G′K+G′xK+G′yK+GxK+Gy,E=∑G′-τG-G′K+Gx2+K+Gy2,F=∑G′εG-G′,P=∑G′χ1G-G′,Q=∑G′χ2G-G′,W=∑G′χ3G-G′K+G′x2+K+G′y2.In order to facilitate solving eigenvalues, (26) is further simplified into standard form (28)F0W-Q-1HE0PuKG′u4G′=ω2uKG′u4G′,where 0 represents zero matrix with same order as H.


## 4.2. Resistive-Inductive Circuit

Complex impedanceZ(ω)=RL+iωL if each resonant shunting circuit is composed of a resistor and an inductor in series. For the resistive-inductive circuit, when RL<2L/C, the oscillation frequency can be estimated as (29)fr=12π1LCp-RL24L2.Equations (10) and (15) give(30)χ1+ωχ4+ω2χ2V=-ωχ5+ω2χ3∂2u3∂x2+∂2u3∂y2,where χ4=iRLCp and χ5=iRLApd31Y2(h+2hp)/(1-μ2). Expanding χjj=1,2,3,4,5 as (17) and substituting (22) into (30) yield(31)∑G′χ1G-G′+ω∑G′χ4G-G′+ω2∑G′χ2G-G′u4G′=ω∑G′χ5G-G′+ω2∑G′χ3G-G′K+G′x2+K+G′y2uKG′.Equations (23) and (31) may be expressed as matrix form (32)HuKG′+Eu4G′=ω2FuKG′,P+ωS+ω2Qu4G′=ωΤ+ω2WuKG′,where S=∑G′χ4(G-G′) and T=∑G′χ5(G-G′). Equation (32) can also be written as standard form [31](33)00F000W-QI0000I00-1HE000P-TS00I0000IZ1Z2Z3Z4=ωZ1Z2Z3Z4,where Z1=uK(G′), Z2=u4(G′), Z3=ωZ1, and Z4=ωZ2.In order to improve the convergence rate and the accuracy of the numerical calculation, technique of the improved PWE [26] is adopted by PWE-PPC-C. Namely, a parameter matrix g is replaced by 1/g-1, where the index −1 means inverse matrix. All parameters η, D, υ, γ, and τ are converted in the same way. It can be noted that the parameter τ1 cannot be equal to zero. The band structure is obtained when K has been traversed along the boundary of the irreducible region of the first Brillouin zone as shown in Figure 2(a).


## 5. Numerical Results and Discussion

For the two-dimensional PPC, geometrical parameters of unit cell are listed in Table1. Epoxy is chosen as the material of the elastic plate with density ρ1=1180 kg/m3, Poisson’s ratio μ1=0.37, and Young’s modulus Y1=4.35 GPa. PZT-5H is considered as the material of the piezoelectric patches. Its density ρ2=7500 kg/m3 and other material parameters are given as(34)s11=16.5,s12=-4.78×10-12m2/N,d31=-274×10-12C/N,ε33=3400ε0,where ε0 is permittivity of free space. The infinitesimal τ1=1.0×10-10 or less can ensure the convergence.Table 1
Geometrical parameters of unit cell of the PPC.
a
1 (mm)
a
2 (mm)
l
1 (mm)
l
2 (mm)
h (mm)
h
p (mm)
80
60
40
30
5
0.2The effectiveness of the PWE-PPC-C method is verified by other computational methods. Firstly, we compare PWE-PPC-C with Transfer Matrix method in solving a one-dimensional PPC-C. As we know, Transfer Matrix method can accurately solve the one-dimensional PPC-C and obtain its dispersion relation. With same structure and circuit parameters as [7], the PPC-C is solved by Transfer Matrix method and PWE-PPC-C, respectively. In PWE-PPC-C, 101 reciprocal vectors (n1=50) are employed to ensure the convergence sufficiently. The dispersion relation curves calculated by these two methods are shown in Figures 3(a) and 3(b), respectively. For results of Transfer Matrix method, the first LR band gap locates at 422.3–424.5 Hz and the Bragg band gap lies in 708–987 Hz; for those of PWE-PPC-C, the first LR band gap is 422.3–424.2 Hz and the Bragg band gap is 708–986 Hz. The comparisons are listed in Table 2. They show that results of PWE-PPC-C are consistent with those of Transfer Matrix method. Yet the voltage mismatch between two portions of the PPC-C leads to calculation error which has not been completely eliminated since τ1 cannot be zero. So the fictitious band gaps still appear as redundant lines in the band structure of PWE-PPC-C.Table 2
Comparisons between PWE-PPC-C and other methods.
Models
Methods
LR band gaps
Bragg band gaps
One-dimensional PPC-C
Transfer Matrix method [7]
422.3–424.5 Hz
708–987 Hz
PWE-PPC-C
422.3–424.2 Hz
708–986 Hz
Two-dimensional PPC-C
Finite Element method [22]
556–563 Hz
752–948 Hz
PWE-PPC-C
559–566 Hz
754–949 HzFigure 3
Dispersion relation of the one-dimensional PPC-C calculated by (a) Transfer Matrix method and (b) PWE-PPC-C, wherea=80 mm, RL=0, and L=0.419 H.
(a)
(b)The PWE-PPC-C inherits the advantage of the PWE in effectively solving more than one-dimensional PC. As the next step, PWE-PPC-C is further verified by Finite Element method in solving two-dimensional PPC-C. Figure4 illustrates band structure of a two-dimensional PPC-C obtained by PWE-PPC-C, where structural, circuit, and other relevant parameters are all the same as those of [22]. Piezoelectric patches are arrayed as tetragonal lattice with lattice constants a1=a2=80 mm. Band structure with all propagation directions is obtained conveniently by PWE-PPC-C as shown in Figure 4. Band structures of two particular orientations, θ=0 and π/4, are computed by Finite Element method [22]. It can be noted that location and width of band gaps of these two orientations obtained by PWE-PPC-C and Finite Element method are in good agreement. Comparisons on band gaps between PWE-PPC-C and Finite Element method are further listed in Table 2. Therefore, PWE-PPC-C can correctly solve two-dimensional PPC-C and obtain its band structure with all propagation directions.Figure 4
Band structure of the two-dimensional PPC-C of tetragonal lattice, wherea=80 mm and L=0.2 H.The proposed PWE-PPC-C is further used to study a two-dimensional PPC-C with orthogonal lattice. Both unit cell and piezoelectric patches are rectangular as shown in Figure1. Piezoelectric patches are arrayed as orthogonal lattice as shown in Figure 2(a). Geometrical parameters of unit cell of the PPC are listed in Table 1. Circuit parameters RL=50 Ω and L=0.25 H. Here n1=n2=15; thus the number of reciprocal vectors is 961. Figure 5 illustrates the band structure of the PPC-C that a complete LR gap exists in all directions and spans from 573 Hz to 584 Hz; two incomplete Bragg gaps appear between 756 and 960 Hz in Γ→X direction and between 1336.2 and 1656.8 Hz in Y→Γ direction, respectively. Directions Γ→X and Y→Γ of reciprocal space correspond to horizontal and vertical directions of primitive space, respectively. It can be found from Figure 5 that the Bragg gap frequency of vertical direction is higher than that of horizontal direction. This is because lattice constant of vertical direction, a2, is smaller than that of horizontal direction, a1. It can be inferred that the Bragg gap frequency decreases with increasing of lattice constant in a certain direction. Therefore, the PPC-C with orthorhombic as well as tetragonal lattice can be calculated by the proposed PWE-PPC-C.Figure 5
Band structure of the two-dimensional PPC-C of orthogonal lattice, whereRL=50 Ω and L=0.25 H.We further investigate the effect of circuit parameters on band gap by the proposed PWE-PPC-C. The LR band gap is often used to isolate vibration because of the low frequency. Therefore, the following discussions are focused on the first LR band gap. Firstly, the inductance is fixed atL=0.25 H. Band structures of the two-dimensional PPC-C of orthogonal lattice with different resistances are shown in Figures 5, 6(a), and 6(b), where resistance RL=50, 200, and 400 Ω, respectively. The center frequency and bandwidth of the LR band gap both decrease with increasing of the resistance. Then the resistance is fixed at RL=50 Ω. Band structures of the two-dimensional PPC-C of orthogonal lattice with different inductances are shown in Figures 6(c), 5, and 6(d), where the inductance L=0.2, 0.25, and 0.30 H, respectively. These figures illustrate that the larger the inductance, the lower the center frequency and the narrower the bandwidth of the LR band gap.Figure 6
LR band gaps of the two-dimensional PPC-C with different circuit parameters: (a)RL=200 Ω, L=0.25 H; (b) RL=400 Ω, L=0.25 H; (c) RL=50 Ω, L=0.20 H; (d) RL=50 Ω, L=0.30 H.
(a)
(b)
(c)
(d)PWE-PPC-C is further applied to calculate dependence of the LR band gap distribution on circuit parameters of the two-dimensional PPC. The LR band gap distribution and estimated oscillation frequency versus the inductance for different resistance are shown in Figure7. The center frequency and bandwidth of the LR band gap both increase with decreasing of the inductance. The tendency is same as that reported in related references. Figure 8 shows that distribution of the LR band gap and estimated oscillation frequency versus the resistance for different inductance. The center frequency and bandwidth of the LR band gap both decrease with increasing of the resistance. Transmission factor is often used to measure the wave propagation properties of the PPC with finite period. Nevertheless, curves of transmission factor versus the frequency show that the larger the resistance, the wider the band gap. The reason is that damping effect is enhanced by increasing of the resistance; thus band gap of near field wave is broadened. One also can find from Figures 7 and 8 that the frequency of LR band gaps is lower than the estimated oscillation frequency of the resonant shunting circuits. This is because the oscillation frequency depends on the dynamic capacitance but not on the static capacitance [32].Figure 7
Distribution of the LR band gap of the two-dimensional PPC-C versus the inductance for different resistance.Figure 8
Distribution of the LR band gap of the two-dimensional PPC-C versus the resistance for different inductance.Finally, the advantage of PWE-PPC-C in effectively calculating the PPC-C of various lattices is demonstrated. Figure9 displays the band structure of the PPC with piezoelectric patches arranged as triangular lattice (see Figure 2(b)), where the lattice constant a=80 mm and the radius of the piezoelectric patches R=30 mm. It can be found that a completely tiny LR band gap exists between 375 and 384 Hz and two Bragg band gaps locate between 1685.9 and 1686.6 Hz in Γ→K direction and between 1131.3 and 1541.5 Hz in Γ→M direction, respectively. Moreover, Figure 10 displays band structure of a PPC-C of hexagonal lattice. The hexagonal lattice (a=80 mm) consists of elastic matrix and piezoelectric patches (R=30 mm) as shown in Figure 2(c). A Bragg band gap appears between 1019 and 1142 Hz along P→Q direction. Therefore, from band structures of four kinds of lattices as shown in Figures 4, 5, 9, and 10, it can be found that LR band gap is a complete band gap in all directions if it exists. But Bragg band gaps only exist in some directions. Moreover, Bragg band gap appears in different direction for different lattice type of the PPC-C.Figure 9
Band structure of the two-dimensional PPC-C of triangular lattice.Figure 10
Band structure of the two-dimensional PPC-C of hexagonal lattice.


## 6. Conclusions

An enhanced plane wave expansion method is proposed to solve the PPC connected with resonant shunting circuits (PWE-PPC-C). A two-dimensional PPC-C model has been presented to demonstrate in detail solving process of PWE-PPC-C. Piezoelectric patches of the two-dimensional PPC are arrayed as orthorhombic lattice with two different lattice constants in horizontal and vertical directions. The PWE-PPC-C method has been validated by two ways. On one hand, PWE-PPC-C has consistent result with Transfer Matrix method in solving one-dimensional PPC-C. On the other hand, PWE-PPC-C agrees well with Finite Element method in solving a two-dimensional PPC-C of tetragonal lattice.The PPC-C of orthorhombic lattice is studied by PWE-PPC-C. Both LR band gaps and Bragg band gaps can be calculated by PWE-PPC-C effectively. Dependence of the first LR band gap on the resistance and the inductance is further investigated by PWE-PPC-C in detail.The band structures of the PPC-C of triangular and hexagonal lattices are obtained by PWE-PPC-C, which demonstrate the advantage of PWE-PPC-C in solving PPC-C of various arrangements. Some redundant lines appear near the LR band gap in the figure of band structure, whose generation mechanism will be further clarified in the future work.



---

*Source: 4015363-2016-03-14.xml*