# Wave Propagation in L-Shape Beams with Piezoelectric Shunting Arrays

Authors: Shengbing Chen; Yubao Song; Hao Zhang
Journal: Shock and Vibration
 (2019)

---

## Abstract

Piezoelectric shunting arrays are employed to control the elastic wave propagation in L-shape beams. Unlike straight beams where longitudinal and flexural waves usually propagate independently, these waves are coupled in an L-shape beam. Based on transfer matrix method and Bloch theorem, dispersion curves and vibration transmissibility are evaluated and analyzed. A locally resonant gap is produced on the flexural and longitudinal waves, respectively, whose locations are nonoverlapped if the shunt damping is void. However, the longitudinal wave band gap can be completely overlaid by the flexural one when a proper shunting resistance is involved. With the decreasing of shunting inductance, the locations of longitudinal and flexural wave gaps both go up to higher frequencies which agree with the variation of resonant frequencies, but they are less affected by shunting resistance. As the resistance increases, the width of the band gaps grows, whereas the attainable maximum attenuation within the band gaps shows a significant decrease. Also, finite element simulations are performed to validate the numerical predictions, which demonstrate that the resulting transmissibility of displacements agree well with the band gaps.


---

## Body

## 1. Introduction

Flexible structures are extensively applied in large aerospace products, which represents a trend as well as great challenge work for industrial and academic communities. Due to low structural damping and lack of air drag, these vibrations or elastic waves propagating in the structures are lightly decayed that leads to fatigue, instability, or other problems [1, 2]. Beam element as one of the principal components of space structures, such as truss and framework, has received considerable attention of engineers and scientists in the noise and vibration control field [3–5]. Therefore, much effort has been exerted on vibration mitigation of flexible beams. Nevertheless, this vibration or elastic wave propagation in the beam-type structures also attracts significant attention, especially when the vibrationless environment is desired, such as payload of high resolution satellites and platform for high accuracy experiments.In the past two decades, elastic wave propagation in periodic composite materials or structures, named phononic crystals (PCs), has received great attention [6–9]. The primary property of PCs is band-gap behavior, i.e., the propagation of elastic waves whose frequencies lie in the band gap will be forbidden. The seminal work of Liu et al. sparked passionate research interest in locally resonant (LR) PCs [8]. They examined three-dimensional PCs consisting of cubic arrays of coated lead spheres immersed in an epoxy matrix and proposed a new type of band-gap formation mechanism, i.e., locally resonant band gap. More recently, the development of LR PCs has been extended to newly emerging field: acoustic metamaterials [10–12]. Acoustic metamaterials are generally regarded as materials with artificial microstructures that possess unusual physical properties such as band gaps, negative refraction, acoustic cloaking, etc.Structures periodically shunted by electrical resonant circuits can also produce locally resonant band gaps. As a result, piezoelectric shunting arrays can be utilized to suppress the vibration or elastic wave propagation in flexible beams. Particularly, the band gaps can be tuned to the desired frequencies conveniently via adjustment of the circuit parameters. The original periodic shunting concept was proposed by Thorp et al. who utilized an array of shunted piezoelectric patches that mounted on a rod to create band gaps centered at the resonant frequency of the shunting circuits [13]. Wang et al. experimentally investigated the low-frequency locally resonant gaps induced by resonant shunting arrays [14]. Airoldi et al. proposed the idea of producing multiple band gaps in a beam with multiresonant shunts [15]. Also, they designed tunable acoustic metamaterials with periodic arrays of resonant shunts [16]. We also performed some theoretical and experimental investigations on flexible structures with arrays of shunting circuits, and some useful conclusions are arrived [17, 18]. These studies demonstrate that the equivalent mechanical properties of host structures can be modified by piezoelectric shunts and suggest the application for developing tunable band gaps. On the whole, longitudinal and flexural waves in beams with shunting arrays are generally examined independently. Thus, longitudinal and flexural band gaps are usually analyzed separately. In practical structures, however, different beams commonly join together by welding or thread connection with certain intersection angle to form a framework or truss. Actually, the flexural and longitudinal waves are usually coupled and propagate simultaneously in such beam structures. Hence, the goal of this work is to investigate the elastic wave propagation in an L-shape beam with piezoelectric shunting arrays and to achieve a broadband suppression of flexural and longitudinal vibrations.In this paper, piezoelectric shunting arrays are periodically bonded to the surfaces of an aluminum L-shape beam to produce a locally resonant band gap over the desired frequency range. Based on the transfer matrix (TM) method and Bloch theorem, numerical investigations are performed in detail. Specifically, the structure configuration and band-gap calculation based on the TM method are described in Section2. Firstly, layout of the assembled beam is illustrated and effective stress equation of shunted piezoelectric patches is derived. Secondly, transfer matrix is obtained from interface continuity conditions. Then, propagation constant is calculated by applying Bloch periodic conditions to the transfer matrix. Also, vibration transmissibility is derived to quantify the capacity of vibration suppression. Numerical results and discussions are demonstrated in Section 3. The dispersion curves are illustrated, where locally resonant band gap is characterized. Parametric analyses of shunting circuits are carried out, and some rules are summarized. Transmissibility of displacements obtained from the TM method and finite element method (FEM) is compared to verify the theoretical calculations. The variation of vibration transmissibility with intersection angle and shunting parameters is also demonstrated. Finally, conclusions are summarized in Section 4.


## 2. Structure Configuration and Mathematical Model

### 2.1. Configuration

The assembled composite beam consists of a substrate L-shape beam and two arrays of shunted piezoelectric patches (Figure1), which are oppositely mounted on the upper and lower surfaces of the beam. Each piezoelectric patch is connected to a resistive-inductive circuit. As a result, the composite beam behaves as a one-dimensional waveguide propagating different waves. For simplicity, we only consider the longitudinal wave and flexural wave in this paper. Apparently, the L-shape beam can be considered as two straight beams joining together with an intersection angle of θ. In addition, the parameters of piezoelectric patches and shunting circuits are identical. Hence, the waveguide will be regarded as superlattice structures, whose wave propagation properties can be conveniently characterized through dynamics analyses of a unit cell (Figure 2) combined with the application of Bloch theorem. The segments without and with piezoelectric patches are denoted as A and B, respectively.Figure 1
Sketch of the composite beam: (a) L-shape beam with piezoelectric patches and (b) connection pattern of the shunting circuit.Figure 2
Sketch of the unit cell.


### 2.2. Modeling of Shunted Piezoelectric Patches

Piezoelectric patches are polarized along with thez-axis, whose surfaces are assumed to be free of constraint except those perpendicular to the x-axis. For convenience, the directions along with x, y, and z axes are denoted by 1, 2, and 3 in the piezoelectric equations, respectively. As a result, the reduced piezoelectric equations are [13](1)S1=s11ET1+d31E3,D3=d31T1+ε33TE3,where S1 and T1 are mechanical strain and stress, respectively; s11E is compliance coefficient at constant electric field; d31 is piezoelectric strain constant coupling electrical and mechanical properties; D3 is electric displacement on the electrodes; ε33T is the permittivity at constant stress; and E3 is the electric field intensity.In the equation set (1), eliminating T1 gives(2)D3=d31s11ES1+ε33SE3,where ε33S is the permittivity at constant strain, which can be expressed as(3)ε33S=ε33T−d312s11E.The electric charge on the electrodes of each piezoelectric patch can be given by(4)Qs=∬AsD3dA,where As is the area of the electrode and As=bl.Substituting equation (2) into equation (4) yields(5)E3=Qsε33SAs−d31bs11Eε33SAs∫0lS1dx.On the other hand, the piezoelectric patch is attached to a resistive-inductive circuit, so applying Ohm’s law yields(6)Qs=−E3hpsR+Ls,where s is the complex variable of Laplace transformation.The beam propagates longitudinal and flexural waves, so mechanical strain in the piezoelectric patch can be written as(7)S1=∂ux,t∂x−h+hp2⋅∂2wx,t∂x2,where ux,t and wx,t are longitudinal and flexural displacements, respectively.Substituting equations (6) and (7) into equation (5) yields(8)E3=−sd31bR+Lss11EhpsR+LsCpS+1ul,t−u0,t−h+hp2∂wl,t∂x−∂w0,t∂x,where CpS is the capacitance of the piezoelectric patch at constant strain and(9)CpS=ε33SAshp.Substituting equation (8) into equation (1) yields(10)T1=S1s11E+Ta,where(11)Ta=sd312bR+Lss11E2hpsR+LsCpS+1ul,t−u0,t−h+hp2∂wl,t∂x−∂w0,t∂x.


### 2.3. Governing Equations of Wave Motion

On the assumption of Euler–Bernoulli beam, the governing equation of longitudinal and flexural motions of the beam element can be given by(12)Ei∂2ux,t∂x2−ρi∂2ux,t∂t2=0,EiIi∂4wx,t∂x4+ρiAi∂2wx,t∂t2=0,where Ei is Young’s modulus, ρi is the density, Ii is the inertia moment of cross section, and Ai is the area of cross section. Here, i = 1, 2, which represent segments A and B, respectively.If harmonic motions of the beam are assumed, the longitudinal and transverse displacements can be written as(13)ux,t=Uxeiωt,wx,t=Wxeiωt,where ω is the angular frequency and i=−1.Substituting equation (13) into equation (12) and solving the differential equations yield(14)Ux=c1cosλix+c2sinλix,Wx=c3cosγix+c4sinγix+c5coshγix+c6sinhγix,where ci,i=1,2,…,6 are undetermined coefficients. λi and γi are wavenumbers of longitudinal and flexural waves, which can be given by(15)λi=ρiEiω,γi=ρiAiω2EiIi4.Therefore, the vector of generalized displacements and forces can be rewritten as(16)Zx=Hxc,where Zx=UxNxWxϑxMxQxT. ϑx is the slope of the flexural displacement. Nx, Mx and Qx are tension, moment, and shear force, respectively. c=c1c2c3c4c5c6T and Hx is the coefficient matrix (refer to Appendix for details).The continuity conditions between segments A and B in thenth cell can be expressed as(17)Zn,A0=Zn,B0,where the first subscript denotes cell number and the second is segment notation.Similarly, the continuity conditions betweennth and (n + 1)th cells also can be written as(18)Zn+1,Al−a=Zn,Bl.Hence, combining equations (17) and (18) gives(19)cn+1=Tcn,where T is the transfer matrix of the unit cell and(20)T=HA−1l−aHBlHB−10HA0.On the other hand, applying the Bloch theorem to the periodic beam gives(21)cn+1=eμcn,where μ is propagation constant, whose real part α and imaginary part β are called attenuation constant and phase constant, respectively. If α = 0, the corresponding wave can propagate without decay; otherwise, the corresponding wave will be attenuated.Combining equations (19) and (21) yields(22)T−eμcn=0.Therefore, the propagation constant can be obtained by solving eigenvalue of the transfer matrix.


### 2.4. Derivation of Vibration Transmissibility

When the L-shape beam is composed of finite periods, it can be divided into two components (Figure3), which consists of M and N cells, respectively. The constraint of the beam is free-free and an excitation force F is applied on the left end with incident angle δ. As a result, the boundary conditions of the L-shape beam are summarized as follows:(1)
On the left end(23)N1L=Fcosδ,M1L=0,Q1L=Fsinδ.(2)
On the right end(24)N2R=0,M2R=0,Q2R=0.Figure 3
Components of the L-shape beam.Here, the first subscript is component number, and the second represents location, i.e., on the left end (L) or right end (R).Substituting equation (14) into equations (23) and (24) and assembling them into matrix form yield(25)Ψc1,1=Fcosθ0Fsinθ,(26)Φc2,N=000,where the first subscript is component number and the second is cell number.On the other hand, the right end motion of component 1 can be expressed as(27)Z1R=HBlc1,M.Substituting equation (19) into equation (27) gives(28)Z1R=HBl⋅TM−1c1,1.Similarly(29)Z2R=HBlc2,N=HBl⋅TN−1c2,1.At the joint of components 1 and 2, it has the relation(30)Z1R=ΘZ2L=ΘHAl−ac2,1,where(31)Θ=−cosθ0−sinθ0000−cosθ000−sinθsinθ0−cosθ0000001000000100sinθ000−cosθ.Hence, combining equations (28), (29), and (30) yields(32)c2,N=Πc1,1,where(33)Π=TN−1ΘHAl−a−1HBl⋅TM−1.Finally, equations (25), (26), and (32) can be cast into assembled form(34)MX=Y,where(35)M=Ψ00ΦΠ−I,X=c1,1c2,N,Y=dcosθdsinθ0000000000T,and I is unit matrix.The solution of equation (34) is(36)X=M−1Y.Consequently, the longitudinal and flexural displacements at the left and right ends can be solved as(37)U1LW1LU2RW2R=XTΛ,where(38)Λ=cosλ1l−a0cosλ2l0cosλ1l−a0sinλ2l00cosγ1l−a0cosγ2l0sinγ1l−a0sinγ2l0coshγ1l−a0coshγ2l0sinhγ1l−a0sinhγ2l.As the L-shape beam is finite, there will be no energy transmission, if damping is void. Therefore, the vibration transmissibility is defined as the transmissibility of displacements, namely,(39)TLL=20lgU2RU1L,(40)TFF=20lgW2RW1L,(41)TFL=20lgU2RW1L,(42)TLF=20lgW2RU1L.Equations (39) and (40) present the displacement transmissibility of longitudinal and flexural waves, while equations (41) and (42) demonstrate the conversions between different wave modes. When the transmissibility is greater than zero, the displacement at the detector end will be larger than that of the source end. Hence, the vibration will be amplified. On the contrary, when the transmissibility is smaller than zero, the displacement at the detector end will be smaller than that of the source end. Hence, the vibration will be decayed.


## 2.1. Configuration

The assembled composite beam consists of a substrate L-shape beam and two arrays of shunted piezoelectric patches (Figure1), which are oppositely mounted on the upper and lower surfaces of the beam. Each piezoelectric patch is connected to a resistive-inductive circuit. As a result, the composite beam behaves as a one-dimensional waveguide propagating different waves. For simplicity, we only consider the longitudinal wave and flexural wave in this paper. Apparently, the L-shape beam can be considered as two straight beams joining together with an intersection angle of θ. In addition, the parameters of piezoelectric patches and shunting circuits are identical. Hence, the waveguide will be regarded as superlattice structures, whose wave propagation properties can be conveniently characterized through dynamics analyses of a unit cell (Figure 2) combined with the application of Bloch theorem. The segments without and with piezoelectric patches are denoted as A and B, respectively.Figure 1
Sketch of the composite beam: (a) L-shape beam with piezoelectric patches and (b) connection pattern of the shunting circuit.Figure 2
Sketch of the unit cell.


## 2.2. Modeling of Shunted Piezoelectric Patches

Piezoelectric patches are polarized along with thez-axis, whose surfaces are assumed to be free of constraint except those perpendicular to the x-axis. For convenience, the directions along with x, y, and z axes are denoted by 1, 2, and 3 in the piezoelectric equations, respectively. As a result, the reduced piezoelectric equations are [13](1)S1=s11ET1+d31E3,D3=d31T1+ε33TE3,where S1 and T1 are mechanical strain and stress, respectively; s11E is compliance coefficient at constant electric field; d31 is piezoelectric strain constant coupling electrical and mechanical properties; D3 is electric displacement on the electrodes; ε33T is the permittivity at constant stress; and E3 is the electric field intensity.In the equation set (1), eliminating T1 gives(2)D3=d31s11ES1+ε33SE3,where ε33S is the permittivity at constant strain, which can be expressed as(3)ε33S=ε33T−d312s11E.The electric charge on the electrodes of each piezoelectric patch can be given by(4)Qs=∬AsD3dA,where As is the area of the electrode and As=bl.Substituting equation (2) into equation (4) yields(5)E3=Qsε33SAs−d31bs11Eε33SAs∫0lS1dx.On the other hand, the piezoelectric patch is attached to a resistive-inductive circuit, so applying Ohm’s law yields(6)Qs=−E3hpsR+Ls,where s is the complex variable of Laplace transformation.The beam propagates longitudinal and flexural waves, so mechanical strain in the piezoelectric patch can be written as(7)S1=∂ux,t∂x−h+hp2⋅∂2wx,t∂x2,where ux,t and wx,t are longitudinal and flexural displacements, respectively.Substituting equations (6) and (7) into equation (5) yields(8)E3=−sd31bR+Lss11EhpsR+LsCpS+1ul,t−u0,t−h+hp2∂wl,t∂x−∂w0,t∂x,where CpS is the capacitance of the piezoelectric patch at constant strain and(9)CpS=ε33SAshp.Substituting equation (8) into equation (1) yields(10)T1=S1s11E+Ta,where(11)Ta=sd312bR+Lss11E2hpsR+LsCpS+1ul,t−u0,t−h+hp2∂wl,t∂x−∂w0,t∂x.


## 2.3. Governing Equations of Wave Motion

On the assumption of Euler–Bernoulli beam, the governing equation of longitudinal and flexural motions of the beam element can be given by(12)Ei∂2ux,t∂x2−ρi∂2ux,t∂t2=0,EiIi∂4wx,t∂x4+ρiAi∂2wx,t∂t2=0,where Ei is Young’s modulus, ρi is the density, Ii is the inertia moment of cross section, and Ai is the area of cross section. Here, i = 1, 2, which represent segments A and B, respectively.If harmonic motions of the beam are assumed, the longitudinal and transverse displacements can be written as(13)ux,t=Uxeiωt,wx,t=Wxeiωt,where ω is the angular frequency and i=−1.Substituting equation (13) into equation (12) and solving the differential equations yield(14)Ux=c1cosλix+c2sinλix,Wx=c3cosγix+c4sinγix+c5coshγix+c6sinhγix,where ci,i=1,2,…,6 are undetermined coefficients. λi and γi are wavenumbers of longitudinal and flexural waves, which can be given by(15)λi=ρiEiω,γi=ρiAiω2EiIi4.Therefore, the vector of generalized displacements and forces can be rewritten as(16)Zx=Hxc,where Zx=UxNxWxϑxMxQxT. ϑx is the slope of the flexural displacement. Nx, Mx and Qx are tension, moment, and shear force, respectively. c=c1c2c3c4c5c6T and Hx is the coefficient matrix (refer to Appendix for details).The continuity conditions between segments A and B in thenth cell can be expressed as(17)Zn,A0=Zn,B0,where the first subscript denotes cell number and the second is segment notation.Similarly, the continuity conditions betweennth and (n + 1)th cells also can be written as(18)Zn+1,Al−a=Zn,Bl.Hence, combining equations (17) and (18) gives(19)cn+1=Tcn,where T is the transfer matrix of the unit cell and(20)T=HA−1l−aHBlHB−10HA0.On the other hand, applying the Bloch theorem to the periodic beam gives(21)cn+1=eμcn,where μ is propagation constant, whose real part α and imaginary part β are called attenuation constant and phase constant, respectively. If α = 0, the corresponding wave can propagate without decay; otherwise, the corresponding wave will be attenuated.Combining equations (19) and (21) yields(22)T−eμcn=0.Therefore, the propagation constant can be obtained by solving eigenvalue of the transfer matrix.


## 2.4. Derivation of Vibration Transmissibility

When the L-shape beam is composed of finite periods, it can be divided into two components (Figure3), which consists of M and N cells, respectively. The constraint of the beam is free-free and an excitation force F is applied on the left end with incident angle δ. As a result, the boundary conditions of the L-shape beam are summarized as follows:(1)
On the left end(23)N1L=Fcosδ,M1L=0,Q1L=Fsinδ.(2)
On the right end(24)N2R=0,M2R=0,Q2R=0.Figure 3
Components of the L-shape beam.Here, the first subscript is component number, and the second represents location, i.e., on the left end (L) or right end (R).Substituting equation (14) into equations (23) and (24) and assembling them into matrix form yield(25)Ψc1,1=Fcosθ0Fsinθ,(26)Φc2,N=000,where the first subscript is component number and the second is cell number.On the other hand, the right end motion of component 1 can be expressed as(27)Z1R=HBlc1,M.Substituting equation (19) into equation (27) gives(28)Z1R=HBl⋅TM−1c1,1.Similarly(29)Z2R=HBlc2,N=HBl⋅TN−1c2,1.At the joint of components 1 and 2, it has the relation(30)Z1R=ΘZ2L=ΘHAl−ac2,1,where(31)Θ=−cosθ0−sinθ0000−cosθ000−sinθsinθ0−cosθ0000001000000100sinθ000−cosθ.Hence, combining equations (28), (29), and (30) yields(32)c2,N=Πc1,1,where(33)Π=TN−1ΘHAl−a−1HBl⋅TM−1.Finally, equations (25), (26), and (32) can be cast into assembled form(34)MX=Y,where(35)M=Ψ00ΦΠ−I,X=c1,1c2,N,Y=dcosθdsinθ0000000000T,and I is unit matrix.The solution of equation (34) is(36)X=M−1Y.Consequently, the longitudinal and flexural displacements at the left and right ends can be solved as(37)U1LW1LU2RW2R=XTΛ,where(38)Λ=cosλ1l−a0cosλ2l0cosλ1l−a0sinλ2l00cosγ1l−a0cosγ2l0sinγ1l−a0sinγ2l0coshγ1l−a0coshγ2l0sinhγ1l−a0sinhγ2l.As the L-shape beam is finite, there will be no energy transmission, if damping is void. Therefore, the vibration transmissibility is defined as the transmissibility of displacements, namely,(39)TLL=20lgU2RU1L,(40)TFF=20lgW2RW1L,(41)TFL=20lgU2RW1L,(42)TLF=20lgW2RU1L.Equations (39) and (40) present the displacement transmissibility of longitudinal and flexural waves, while equations (41) and (42) demonstrate the conversions between different wave modes. When the transmissibility is greater than zero, the displacement at the detector end will be larger than that of the source end. Hence, the vibration will be amplified. On the contrary, when the transmissibility is smaller than zero, the displacement at the detector end will be smaller than that of the source end. Hence, the vibration will be decayed.


## 3. Numerical Results and Discussions

As an example, aluminum and PZT-5H are selected as the materials of substrate beam and piezoelectric patches, respectively. The density, Young’s modulus, and Poisson’s ratio of aluminum are 2700 kg/m3, 70 GPa, and 0.3, respectively. The material parameters of PZT-5H are listed in Table 1. The geometrical parameters are summarized as follows: l = 30 mm, b = 20 mm, a = 40 mm, h = 2 mm, and hp = 0.6 mm.Table 1
Material parameters of PZT-5H.
Parameters
Value
Density
7500 kg/m−3
s
11
E
16.5 × 10−12 m2/N
d
31
−2.74 × 10−10 C/N
ε
33
T
3.01 × 10−8 F/m


### 3.1. Dispersion Relations

The band gaps generated in the composite beam are not only produced by local resonances of the shunts but also affected by the shunt damping when electrical resistance is involved. In order to demonstrate the band gaps induced purely by local resonances of the shunting circuits, the electrical resistanceR is set to be zero at first. Hence, the propagation constant (when L = 0.4 H) of the periodic composite beam is illustrated in Figure 4. One locally resonant gap (shaded area in inset) is induced in longitudinal and flexural waves, respectively. The band gap of longitudinal wave is a pure attenuation zone (PAZ) ranging from 1528 Hz to 1558 Hz, in which the propagation constant is a pure real number. The maximum attenuation occurs at the lower bounding frequency of the band gap, which is a typical feature observed in conventional locally resonant crystals. Also, the phase constant curve near the lower bounding frequency is rather flat, which means that the group velocity of longitudinal wave is close to zero. However, the band gap of flexural wave shows great difference. It is the combination of a pure attenuation zone (PAZ) and a complex attenuation zone (CAZ), whose propagation constant is a complex number. Unlike longitudinal wave which owns one single wave mode, the flexural wave is more complicated which possesses two wave modes, i.e., propagation wave and evanescent wave. Generally, the evanescent wave is neglected when the waveguide is large enough to decay it quickly. However, there is an exception when a CAZ is produced by the coupling effect of the propagation wave and evanescent wave. In the CAZ, the propagation wave and evanescent wave are coupled together and transfer energy to each other. This is the reason why the propagation constant is complex in the CAZ, though no damping is involved. Actually, there is no vibration energy transmitting between cells in the CAZ [19]. The resulting band gap of flexural wave is located in a lower frequency scope, ranging from 1465 Hz to 1521 Hz. Obviously, the maximum attenuation does not appear at the lower bounding frequency but coincides with the interface of PAZ and CAZ instead. In addition, the phase constant curve near the lower bounding frequency is not as flat as that of the longitudinal wave. That is to say, in the flexural wave band gap, some of the vibration energy is easier to be confined in the matrix, rather than almost completely restricted to the vibrators. Though both the flexural wave band gap and longitudinal wave band gap are induced by the same resonant circuits, the locations of the two band gaps are not overlapped. The reason is that the dynamic capacitance of a piezoelectric patch is variable with constraint conditions. In other words, the resonance frequencies of shunting circuits are different between longitudinal and flexural cases. The actual resonance frequency is affected by the motion of the substrate. Generally, the real effective capacitance lies between the capacitance at constant strain and that at constant stress, i.e.,(43)Ceff∈CpS,CpT.Figure 4
Propagation constant of the composite beam.
(a)
(b)The variations of longitudinal wave gap (LG) and flexural wave gap (FG) with shunting inductance are illustrated in Figure5. With the decreasing of inductance L, the locations of LG and FG both go up to higher frequencies. It is aforementioned that locally resonant gap is produced by local resonance of scatters and thus coincides with the resonance frequency of local vibrators. As for electrical resonators, the resonance frequency is(44)fr=12πLCeff.Figure 5
Variation of band gaps with shunting inductance.Evidently, equation (44) also demonstrates that the frequency of band gaps will increase with the decreasing of L. Moreover, it is observed that the two band gaps (LG and FG) keep distinct and nonoverlapping.The effect of shunt damping (R = 100 Ω) on the band gaps is demonstrated in Figure 6. The impacts on longitudinal and flexural wave band gaps are quite similar. Firstly, the attainable maximum attenuation within the band gaps shows a significant decrease when shunting resistance is involved. Evidently, the reason is that local resonance of oscillators will be suppressed by resistance damping. In addition, the phase fluctuations within the band gaps become much smoother. Secondly, the damping effect at the bounding frequencies of band gaps is substantial, resulting in much broader attenuation zones than those merely induced by local resonances. Moreover, the damping capacities at the upper and lower bounding frequencies are similar in magnitude. It is worthy of note that flexural and longitudinal wave band gaps become overlapped. Nevertheless, the frequencies corresponding to the maximum attenuation constant are inconsistent, whose difference Δf is about 39 Hz.Figure 6
Propagation constant with shunt damping.
(a)
(b)The variation of band gaps with shunting resistance (L = 0.4 H) is shown in Figure 7. The location of band gaps keeps stationary, while the band width increases with the increase of resistance R. This is also in accordance with equation (42), which demonstrates that the resonance frequency is unaffected by shunting resistance. In addition, as the shunting resistance increases, longitudinal and flexural wave gaps are gradually overlapped. For example, the two band gaps will be partially overlapped when R = 20 Ω (see inset of Figure 7).Figure 7
Variation of band gaps with shunting resistance.


### 3.2. Vibration Transmissibility

Finite element method (FEM) is employed to validate the theoretical calculation of TM method. An L-shape beam whose intersection angleθ = 90° is selected as the benchmark model. Both two components of the beam consist of three cells, i.e, M = N = 3, which are built in commercial FEM software COMSOL as shown in Figure 8. The substrate beam and piezoelectric patches are meshed with quadratic hexahedral elements, while the shunting circuits are defined in SPICE circuit editor. An excitation force F of unit amplitude is applied on the left end of the beam, whose incident angle δ is 45°. The analysis type is frequency response, and direct SPOOLES is selected as the linear system solver.Figure 8
Finite element model.The comparison between finite element simulation and theoretical calculation is illustrated in Figure9. The vibration transmissibility without shunting effect, i.e., the circuits are open, is presented in Figures 9(a) and 9(b). On the transmission curves of flexural wave, five peaks appear in the frequency range of 0∼2000 Hz. In contrast, two transmission dips and three peaks occur on the transmission curves of longitudinal wave. A peak denotes that the displacement of the detector end is much larger than that of the source end. So, the beam strongly amplifies the vibration corresponding to the peak. On the contrary, a dip denotes that the displacement of the detector end is much smaller than that of the source end. That is to say, the beam strongly decays the vibration corresponding to the dip. When the shunting circuits are functioning, the two flexural wave transmission peaks (Figure 9(c)) around the band-gap location are suppressed and become much smoother. That is to say, the shunting arrays can effectively decay the flexural wave that falls in the band gap. Similarly, the shunting arrays can also smooth the transmission dips (Figure 9(d)) at the location of the band gap of longitudinal wave. On the whole, the theoretical result by TM method and numerical result by FEM match well, whose small difference may be caused by different modeling postulates. The model of FEM is based on two-dimensional plane-stress assumption, whereas the model of TM method is based on Euler–Bernoulli beam assumption.Figure 9
Comparison between TM method and FEM. (a, b) Vibration transmissibility of flexural and longitudinal waves with open circuits. (c, d) Vibration transmissibility of flexural and longitudinal waves with closed circuits, whereL = 0.4 H and R = 800 Ω.
(a)
(b)
(c)
(d)The variation of vibration transmissibility with intersection angle is demonstrated in Figure10, which includes three cases, i.e., acute angle (θ = 60°), blunt angle (θ = 120°), and flat angle (θ = 180°). For comparison, the attenuation constants with closed shunting circuits are illustrated in Figures 10(a) and 10(b). As the shunting resistance is rather large (R = 800 Ω), the attenuation zones of flexural and longitudinal waves are very wide with the compromise of a smaller decay magnitude. On the whole, the variation of vibration transmissibility with intersection angle is dramatic. Not only the location of transmission peaks and dips but also the shape of curves changes a lot. However, the invariable is that the attenuation zones on the transmission curves keep stationary, which always coincide with the location of band gaps. Note that the flexural and longitudinal waves are coupled when θ = 60° and θ = 120°. Actually, when the intersection angle is not flat, the longitudinal wave and flexural wave are coupled, and longitudinal waves can be converted to flexural waves whose attenuation effect is stronger. Consequently, the attenuation on vibration transmissibility of longitudinal waves is still substantial, even though the attenuation constant of longitudinal wave is small. When θ = 180°, the beam becomes straight, so the flexural and longitudinal waves will be decoupled. The attenuation of flexural wave induced by shunting arrays remains distinct as shown in Figure 10(g). In contrast, the influence of shunting circuits on vibration transmissibility of longitudinal wave is less notable (Figure 10(h)). The reason for this is that the attenuation constant of longitudinal wave is much smaller and there is no transmission peak located within the frequency scope of band gap. In a word, the attenuation effect of shunting arrays is always remarkable with different intersection angles, particularly when the angle is not flat.Figure 10
Variation of vibration transmissibility with intersection angleθ. For comparison, both attenuation constants and vibration transmissibility of flexural and longitudinal waves with open circuits (OC) and closed circuits (CC, L = 0.4 H and R = 800 Ω) are presented.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)The variation of vibration transmissibility with circuit parameters is illustrated in Figure11. The attenuation effect of shunting circuit is weak when the frequency is far away from the band gap, so the graphs are merely plotted over the range of 1000∼2000 Hz. Figures 11(a) and 11(b) demonstrate the variation of transmission curves with shunting resistance. The effect of shunting circuits is apparent. Compared with the open circuits (OC), the two transmission peaks of flexural wave are squeezed away from each other when the circuits are closed. The reason is that resonant shunting circuit will reduce and increase the effective modulus of the piezoelectric patches before and after the resonant frequency, respectively [20]. Note that the shunting circuits generate a new peak on the longitudinal transmission curve. The peak is induced by resonance of shunting circuits, which should coincide with the bounding frequency of band gaps. With the increasing of resistance R, the transmission peaks become lower, while the transmission dips get higher. The reason is that the shunt damping will increase with resistance and suppress the resonant oscillation of shunting circuits. However, the influence on locations is small because the resonant frequency of shunting circuit is less affected by resistance. Figures 11(c) and 11(d) demonstrate the variation of transmission curves with shunting inductance. It is observed that the location of transmission trough is significantly affected by the inductance L. As the inductance increases, the trough goes to a lower frequency. The reason is that the resonant frequency (refer to equation (44)) of shunting circuits will decrease with the increase of inductance L. This rule agrees well with the result in Figure 5.Figure 11
Variation of vibration transmissibility with circuit parameters.L = 0.4 H in (a, b). R = 800 Ω in (c, d).
(a)
(b)
(c)
(d)


### 3.3. Mode Conversion

The longitudinal and flexural waves are coupled (whenθ < 180°), so the two wave modes can convert to each other. In order to characterize the mode conversion from longitudinal wave to flexural wave, the incident angle is set to be zero (δ = 0). Hence, there is only longitudinal wave excitation. The vibration transmissibility is plotted in Figure 12. Specifically, the transmissibility of longitudinal wave is demonstrated in Figures 12(a) and 12(c) with different intersection angles (θ = 60° and θ = 120°), respectively. The characterization of longitudinal wave converted to flexural wave is presented in Figures 12(b) and 12(d). It is observed that flexural wave has been produced, though the excitation force is fully longitudinal. Both the longitudinal wave and the converted flexural wave are decayed by the shunting circuits (CC). With different intersection angles (acute and blunt), the attenuation of shunt circuits is rather similar.Figure 12
Mode conversion from longitudinal to flexural waves with open circuits (OC) and closed circuits (CC,L = 0.4 H and R = 800 Ω).
(a)
(b)
(c)
(d)Similarly, to characterize the conversion from flexural wave to longitudinal wave, the incident angle is set to be right angle (δ = 90°). Therefore, there is merely flexural wave excitation. The vibration transmissibility is presented in Figure 13. The transmissibility of flexural wave is demonstrated in Figures 13(a) and 13(c) with different intersection angles (θ = 60° and θ = 120°), respectively. The characterization of flexural wave converted to longitudinal wave is presented in Figures 13(b) and 13(d). One can find that longitudinal wave has also been induced, though the excitation force is completely flexural. Also, both the flexural wave and the converted longitudinal wave are simultaneously decayed by the shunting circuits (CC). However, some of the transmission peaks of longitudinal wave are not fully developed.Figure 13
Mode conversion from flexural to longitudinal waves with open circuits (OC) and closed circuits (CC,L = 0.4 H and R = 800 Ω).
(a)
(b)
(c)
(d)In a word, the flexural wave and longitudinal wave propagating in the L-shape beam are coupled and coexisting, even though the excitation at the source end is single mode.


## 3.1. Dispersion Relations

The band gaps generated in the composite beam are not only produced by local resonances of the shunts but also affected by the shunt damping when electrical resistance is involved. In order to demonstrate the band gaps induced purely by local resonances of the shunting circuits, the electrical resistanceR is set to be zero at first. Hence, the propagation constant (when L = 0.4 H) of the periodic composite beam is illustrated in Figure 4. One locally resonant gap (shaded area in inset) is induced in longitudinal and flexural waves, respectively. The band gap of longitudinal wave is a pure attenuation zone (PAZ) ranging from 1528 Hz to 1558 Hz, in which the propagation constant is a pure real number. The maximum attenuation occurs at the lower bounding frequency of the band gap, which is a typical feature observed in conventional locally resonant crystals. Also, the phase constant curve near the lower bounding frequency is rather flat, which means that the group velocity of longitudinal wave is close to zero. However, the band gap of flexural wave shows great difference. It is the combination of a pure attenuation zone (PAZ) and a complex attenuation zone (CAZ), whose propagation constant is a complex number. Unlike longitudinal wave which owns one single wave mode, the flexural wave is more complicated which possesses two wave modes, i.e., propagation wave and evanescent wave. Generally, the evanescent wave is neglected when the waveguide is large enough to decay it quickly. However, there is an exception when a CAZ is produced by the coupling effect of the propagation wave and evanescent wave. In the CAZ, the propagation wave and evanescent wave are coupled together and transfer energy to each other. This is the reason why the propagation constant is complex in the CAZ, though no damping is involved. Actually, there is no vibration energy transmitting between cells in the CAZ [19]. The resulting band gap of flexural wave is located in a lower frequency scope, ranging from 1465 Hz to 1521 Hz. Obviously, the maximum attenuation does not appear at the lower bounding frequency but coincides with the interface of PAZ and CAZ instead. In addition, the phase constant curve near the lower bounding frequency is not as flat as that of the longitudinal wave. That is to say, in the flexural wave band gap, some of the vibration energy is easier to be confined in the matrix, rather than almost completely restricted to the vibrators. Though both the flexural wave band gap and longitudinal wave band gap are induced by the same resonant circuits, the locations of the two band gaps are not overlapped. The reason is that the dynamic capacitance of a piezoelectric patch is variable with constraint conditions. In other words, the resonance frequencies of shunting circuits are different between longitudinal and flexural cases. The actual resonance frequency is affected by the motion of the substrate. Generally, the real effective capacitance lies between the capacitance at constant strain and that at constant stress, i.e.,(43)Ceff∈CpS,CpT.Figure 4
Propagation constant of the composite beam.
(a)
(b)The variations of longitudinal wave gap (LG) and flexural wave gap (FG) with shunting inductance are illustrated in Figure5. With the decreasing of inductance L, the locations of LG and FG both go up to higher frequencies. It is aforementioned that locally resonant gap is produced by local resonance of scatters and thus coincides with the resonance frequency of local vibrators. As for electrical resonators, the resonance frequency is(44)fr=12πLCeff.Figure 5
Variation of band gaps with shunting inductance.Evidently, equation (44) also demonstrates that the frequency of band gaps will increase with the decreasing of L. Moreover, it is observed that the two band gaps (LG and FG) keep distinct and nonoverlapping.The effect of shunt damping (R = 100 Ω) on the band gaps is demonstrated in Figure 6. The impacts on longitudinal and flexural wave band gaps are quite similar. Firstly, the attainable maximum attenuation within the band gaps shows a significant decrease when shunting resistance is involved. Evidently, the reason is that local resonance of oscillators will be suppressed by resistance damping. In addition, the phase fluctuations within the band gaps become much smoother. Secondly, the damping effect at the bounding frequencies of band gaps is substantial, resulting in much broader attenuation zones than those merely induced by local resonances. Moreover, the damping capacities at the upper and lower bounding frequencies are similar in magnitude. It is worthy of note that flexural and longitudinal wave band gaps become overlapped. Nevertheless, the frequencies corresponding to the maximum attenuation constant are inconsistent, whose difference Δf is about 39 Hz.Figure 6
Propagation constant with shunt damping.
(a)
(b)The variation of band gaps with shunting resistance (L = 0.4 H) is shown in Figure 7. The location of band gaps keeps stationary, while the band width increases with the increase of resistance R. This is also in accordance with equation (42), which demonstrates that the resonance frequency is unaffected by shunting resistance. In addition, as the shunting resistance increases, longitudinal and flexural wave gaps are gradually overlapped. For example, the two band gaps will be partially overlapped when R = 20 Ω (see inset of Figure 7).Figure 7
Variation of band gaps with shunting resistance.


## 3.2. Vibration Transmissibility

Finite element method (FEM) is employed to validate the theoretical calculation of TM method. An L-shape beam whose intersection angleθ = 90° is selected as the benchmark model. Both two components of the beam consist of three cells, i.e, M = N = 3, which are built in commercial FEM software COMSOL as shown in Figure 8. The substrate beam and piezoelectric patches are meshed with quadratic hexahedral elements, while the shunting circuits are defined in SPICE circuit editor. An excitation force F of unit amplitude is applied on the left end of the beam, whose incident angle δ is 45°. The analysis type is frequency response, and direct SPOOLES is selected as the linear system solver.Figure 8
Finite element model.The comparison between finite element simulation and theoretical calculation is illustrated in Figure9. The vibration transmissibility without shunting effect, i.e., the circuits are open, is presented in Figures 9(a) and 9(b). On the transmission curves of flexural wave, five peaks appear in the frequency range of 0∼2000 Hz. In contrast, two transmission dips and three peaks occur on the transmission curves of longitudinal wave. A peak denotes that the displacement of the detector end is much larger than that of the source end. So, the beam strongly amplifies the vibration corresponding to the peak. On the contrary, a dip denotes that the displacement of the detector end is much smaller than that of the source end. That is to say, the beam strongly decays the vibration corresponding to the dip. When the shunting circuits are functioning, the two flexural wave transmission peaks (Figure 9(c)) around the band-gap location are suppressed and become much smoother. That is to say, the shunting arrays can effectively decay the flexural wave that falls in the band gap. Similarly, the shunting arrays can also smooth the transmission dips (Figure 9(d)) at the location of the band gap of longitudinal wave. On the whole, the theoretical result by TM method and numerical result by FEM match well, whose small difference may be caused by different modeling postulates. The model of FEM is based on two-dimensional plane-stress assumption, whereas the model of TM method is based on Euler–Bernoulli beam assumption.Figure 9
Comparison between TM method and FEM. (a, b) Vibration transmissibility of flexural and longitudinal waves with open circuits. (c, d) Vibration transmissibility of flexural and longitudinal waves with closed circuits, whereL = 0.4 H and R = 800 Ω.
(a)
(b)
(c)
(d)The variation of vibration transmissibility with intersection angle is demonstrated in Figure10, which includes three cases, i.e., acute angle (θ = 60°), blunt angle (θ = 120°), and flat angle (θ = 180°). For comparison, the attenuation constants with closed shunting circuits are illustrated in Figures 10(a) and 10(b). As the shunting resistance is rather large (R = 800 Ω), the attenuation zones of flexural and longitudinal waves are very wide with the compromise of a smaller decay magnitude. On the whole, the variation of vibration transmissibility with intersection angle is dramatic. Not only the location of transmission peaks and dips but also the shape of curves changes a lot. However, the invariable is that the attenuation zones on the transmission curves keep stationary, which always coincide with the location of band gaps. Note that the flexural and longitudinal waves are coupled when θ = 60° and θ = 120°. Actually, when the intersection angle is not flat, the longitudinal wave and flexural wave are coupled, and longitudinal waves can be converted to flexural waves whose attenuation effect is stronger. Consequently, the attenuation on vibration transmissibility of longitudinal waves is still substantial, even though the attenuation constant of longitudinal wave is small. When θ = 180°, the beam becomes straight, so the flexural and longitudinal waves will be decoupled. The attenuation of flexural wave induced by shunting arrays remains distinct as shown in Figure 10(g). In contrast, the influence of shunting circuits on vibration transmissibility of longitudinal wave is less notable (Figure 10(h)). The reason for this is that the attenuation constant of longitudinal wave is much smaller and there is no transmission peak located within the frequency scope of band gap. In a word, the attenuation effect of shunting arrays is always remarkable with different intersection angles, particularly when the angle is not flat.Figure 10
Variation of vibration transmissibility with intersection angleθ. For comparison, both attenuation constants and vibration transmissibility of flexural and longitudinal waves with open circuits (OC) and closed circuits (CC, L = 0.4 H and R = 800 Ω) are presented.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)The variation of vibration transmissibility with circuit parameters is illustrated in Figure11. The attenuation effect of shunting circuit is weak when the frequency is far away from the band gap, so the graphs are merely plotted over the range of 1000∼2000 Hz. Figures 11(a) and 11(b) demonstrate the variation of transmission curves with shunting resistance. The effect of shunting circuits is apparent. Compared with the open circuits (OC), the two transmission peaks of flexural wave are squeezed away from each other when the circuits are closed. The reason is that resonant shunting circuit will reduce and increase the effective modulus of the piezoelectric patches before and after the resonant frequency, respectively [20]. Note that the shunting circuits generate a new peak on the longitudinal transmission curve. The peak is induced by resonance of shunting circuits, which should coincide with the bounding frequency of band gaps. With the increasing of resistance R, the transmission peaks become lower, while the transmission dips get higher. The reason is that the shunt damping will increase with resistance and suppress the resonant oscillation of shunting circuits. However, the influence on locations is small because the resonant frequency of shunting circuit is less affected by resistance. Figures 11(c) and 11(d) demonstrate the variation of transmission curves with shunting inductance. It is observed that the location of transmission trough is significantly affected by the inductance L. As the inductance increases, the trough goes to a lower frequency. The reason is that the resonant frequency (refer to equation (44)) of shunting circuits will decrease with the increase of inductance L. This rule agrees well with the result in Figure 5.Figure 11
Variation of vibration transmissibility with circuit parameters.L = 0.4 H in (a, b). R = 800 Ω in (c, d).
(a)
(b)
(c)
(d)


## 3.3. Mode Conversion

The longitudinal and flexural waves are coupled (whenθ < 180°), so the two wave modes can convert to each other. In order to characterize the mode conversion from longitudinal wave to flexural wave, the incident angle is set to be zero (δ = 0). Hence, there is only longitudinal wave excitation. The vibration transmissibility is plotted in Figure 12. Specifically, the transmissibility of longitudinal wave is demonstrated in Figures 12(a) and 12(c) with different intersection angles (θ = 60° and θ = 120°), respectively. The characterization of longitudinal wave converted to flexural wave is presented in Figures 12(b) and 12(d). It is observed that flexural wave has been produced, though the excitation force is fully longitudinal. Both the longitudinal wave and the converted flexural wave are decayed by the shunting circuits (CC). With different intersection angles (acute and blunt), the attenuation of shunt circuits is rather similar.Figure 12
Mode conversion from longitudinal to flexural waves with open circuits (OC) and closed circuits (CC,L = 0.4 H and R = 800 Ω).
(a)
(b)
(c)
(d)Similarly, to characterize the conversion from flexural wave to longitudinal wave, the incident angle is set to be right angle (δ = 90°). Therefore, there is merely flexural wave excitation. The vibration transmissibility is presented in Figure 13. The transmissibility of flexural wave is demonstrated in Figures 13(a) and 13(c) with different intersection angles (θ = 60° and θ = 120°), respectively. The characterization of flexural wave converted to longitudinal wave is presented in Figures 13(b) and 13(d). One can find that longitudinal wave has also been induced, though the excitation force is completely flexural. Also, both the flexural wave and the converted longitudinal wave are simultaneously decayed by the shunting circuits (CC). However, some of the transmission peaks of longitudinal wave are not fully developed.Figure 13
Mode conversion from flexural to longitudinal waves with open circuits (OC) and closed circuits (CC,L = 0.4 H and R = 800 Ω).
(a)
(b)
(c)
(d)In a word, the flexural wave and longitudinal wave propagating in the L-shape beam are coupled and coexisting, even though the excitation at the source end is single mode.


## 4. Conclusions

Band gaps and vibration transmissibility of the proposed composite beam which consists of an aluminum L-shape beam and periodic arrays of shunted piezoelectric patches are analyzed theoretically. It is found that a locally resonant band gap is generated on the flexural and longitudinal waves, respectively. The two band gaps are not overlapped if the shunt damping is void. Though they are induced by the same resonant circuits, the dynamic capacitance of a piezoelectric patch is variable with the boundary constraints. With the decreasing of shunting inductance, the locations of longitudinal and flexural wave gaps both go up to higher frequencies, which agree with the variation of resonant frequencies. When a proper shunting resistance is involved, the two band gaps can be completely overlapped. However, the attainable maximum attenuation within the band gaps shows a significant decrease. Finite element simulations are performed to validate the numerical predictions, which demonstrate that the resulting vibration transmissibility agrees well with the band gaps. Although the variation of vibration transmissibility with intersection angle is dramatic, the attenuation zones on the transmission curves keep stationary, which always coincide with the location of band gaps. In comparison with the straight beam, the attenuation of longitudinal vibration in the L-shape beam is more efficient because part of the longitudinal vibration will be transformed into flexural vibration, whose attenuation constant is much larger. In summary, piezoelectric shunting arrays can generate a broadband locally resonant gap and effectively decay the flexural and longitudinal waves propagating in the L-shape beam.



---

*Source: 6264251-2019-04-11.xml*