Document Type : Research Paper
Authors
Department of Mechanical Engineering, Arak University, Arak, Iran
Abstract
Keywords

Mechanics of Advanced Composite Structures 3 (2016) 8998 

Semnan University 
Mechanics of Advanced Composite Structures journal homepage: http://MACS.journals.semnan.ac.ir 
Analytical Solution for Sound Radiation of Vibrating Circular Plates coupled with Piezoelectric Layers
K. Khorshidi ^{*}, M. Pagoli
Department of Mechanical Engineering, Arak University, Arak, Iran
Paper INFO 

ABSTRACT 
Paper history: Received: 20160325 Revised: 20160821 Accepted: 20160828 
In the present study, the classical plate theory (CPT) was used to study sound radiation of forced vibrating thin circular plates coupled with piezoelectric layers using simply supported and clamped boundary conditions. The novelty of the study consists of an exact closedform solution that was developed without any use of approximation. Piezoelectric, electrical potential loaded in the transverse direction satisfied the electric boundary conditions (open circuit) and Maxwell's electricity equation. It was assumed that no fluid loading occurred on the plate structure. The sound pressure and the sound power of the radiator were analytically obtained in a far field by using the Rayleigh integral. The proposed analytical method was validated using available data from the literature. Additionally, a few 2D plots of the directivity pattern were illustrated for thin circular plates coupled with piezoelectric layers. Finally, the effect of boundary conditions, piezoelectric thickness, and the piezoelectric layer on the acoustical parameters were examined and discussed in details. 



Keywords: Circular plates Sound power Classical plate theory Piezoelectric layer 


© 2016 Published by Semnan University Press. All rights reserved. 
Thin plates are important structural elements that are widely used in engineering applications. In particular, circular plates are employed in a wide range of engineering applications, such as aerospace structures, vehicles, ship and submarine bodies, airplane bodies and wings, and missiles.
According to classical plate theory (CPT), Wang et al. [1] analyzed the free vibration of piezoelectric coupled circular plates with clamped and simply supported boundary conditions. Their proposed model provided a good structure from which to obtain the distribution of the electric potential in piezoelectric layers. An analytical model for vibration of circular plates with piezoelectric layers was also proposed by Hagood and McFarland [2]. They assumed that the potential distribution of the circular plate in the radial direction was uniform. Heyliger and Ramirez [3] obtained free vibration characteristics of a piezoelectric circular plate, using a discrete layer employing the harmonic equations of motion. Due to the development of smart structures, elastic plates with piezoelectric elements have been studied by several researchers. In those studies, the piezoelectric layers were used as actuators and sensors to effectively control noise and vibration of these flexible structures. Therefore, the modeling and analysis of plates using piezoelectric layers are important. HosseiniHashemi et al. [4] and Khorshidi et al. [5] analyzed the vibration of piezoelectric, coupled, thick circular/annular functionally graded plates (FGPs) subjected to different combinations of soft simply supported, hard simply supported, and clamped boundary conditions, on the basis of Mindlin’s firstorder shear deformation theory (FSDT). Those studies comprehensively investigated the coupling effects between inplane and outofplane displacements and boundary conditions.
Unwanted vibrations in each of these structures can result in their destruction. In addition, vibration of the plates, especially thin plates, can often be a source of noise. Excessive noise can have a serious effect on the mechanical and electrical systems by producing vibration, stress, fatigue, and failure. This sound can be created by force or dynamic torque that is imposed on the plate. Therefore, studying the vibrational behaviors, especially the sound radiated by vibrating plates, to optimize structural design, remains a significant endeavor, and the formulation of sound power reduction appears to be necessary.
Early research efforts to calculate the radiated sound field for a flat plate in an infinite baffle dates back to the work of Rayleigh [6]. Lomas and Hayek [7] developed the Green function solution for the steadystate vibrations of an elastically supported rectangular plate coupled to a semiinfinite acoustic medium. Cremer and Heckl [8] analyzed the sound radiation of a planar source, using a Fourier transformation approach in Kspace (wave number space). Berry [9] studied the sound radiation from a rectangular plate with different boundary conditions. In his work, both free and forced vibrations were considered, and the acoustic emission in the far field and sound power were calculated using a onedimensional (1D) Fourier transform. Yoo [10] investigated the general characteristics of sound radiation, using a point force excited rectangular plate with different boundary edge conditions. In this study, a plate under guided boundary conditions and with the critical frequency only at high frequencies was related to radiation efficiency. This approach was similar to the results that Berry achieved for one force excitation. Rdzanek [11] used closed path integral techniques and stationary phase methods to find the active and passive acoustic powers. Exact acoustical analysis of sound radiation from the free vibration of rectangular Mindlin plates, was performed by Khorshidi [12] and HosseiniHashemi et al. [13]. Vibration and sound radiation of a plateended cylindrical shells were studied by Hongqiu [14]. Zhou et al. [15] calculated the sound pressure of a thin infinite plate in contact with a layered inhomogeneous fluid subject to single point excitation. Sound radiation by the vibrational modes of baffled, flat plates was analyzed by Frank and Paolo [16]. Zhang et al. [17] obtained sound radiation from a baffled rectangular plate, which was expressed using the plate displacement function as a 2D Fourier cosine series and the unknown Fourier expansion coefficients were determined by using the Rayleigh–Ritz procedure. Lee and Singh [18] studied the sound radiated by a thin annular plate in a far field. Lee and Singh [19] studied moderately thick, annular plates, using Mindlin plate theory in a far field. They also worked on sound radiation of thick plates and used the results to obtain a semianalytical method for calculating sound radiation of a broken disk rotor [20].
The analytical analysis of the sound radiation of a thin circular plate exists in the literature, but this study investigated an exact closedform formulation for the sound radiation of thin circular plates using piezoelectric layers. The novelty of the paper is that the exact closedform solution was developed without use of approximation. The exact closedform solution presented in this study can be used in the solution of benchmark problems for validation of future numerical methods.
First, it is assumed that no fluid loading occurs on the plate structure. To study the forced transverse vibration of thin circular plates coupled with piezoelectric layers, the equations of motion were derived based on classical plate theory. Structuralacoustic coupling was implemented for vibrating plate models. The radiation field of a vibrating plate with a specified distribution of velocity on the surface can be computed using the Rayleigh integral approach. The acoustic pressure distribution and the sound powers of the radiator were analytically obtained in its far field. The proposed analytical method was validated using available data from the literature. Additionally, a few 2D plots of the directivity pattern are illustrated for the thin circular plates coupled with piezoelectric layers. Finally, the effect of boundary conditions, piezoelectric thickness, and the piezoelectric layer on the acoustical parameters were examined and discussed in detail.
A flat, piezoelectric coupled circular plate, including one host layer in the middle and two identical piezoelectric layers bonded perfectly to the upper and lower surfaces of the host layer, with outer radius , host plate thickness , and piezoelectric layer thickness , are considered as shown in Figure 1.
Both the top and bottom surfaces of each piezoelectric layer are fully covered by electrodes, which are shortly connected. The thickness of the electrodes is assumed to be extremely small compared to the plate thickness. Thus, in the following formulation, the mechanical effects of the electrodes were neglected. Both piezoelectric layers were polarized perpendicular to the midplane in the positive direction of the plate’s thickness.
Figure 1. A circular plate coupled with piezoelectric layers, showing the coordinate and displacement systems. 
The plate geometry and dimensions were defined in an orthogonal cylindrical coordinate system , to extract mathematical formulations. The origin of the coordinate system was taken at the center of the plate in the middle plane. For convenience in the formulation, the suffixes ‘p’ and ‘h’ were used to denote each piezoelectric layer and the host structure, respectively.
According to the classical plate theory (CLP), the displacement field in the absence of inplane displacements are used as follows:
, 
(1ac) 
, 

, 
where , and are the displacements in the , radial, and tangential directions, respectively, and denotes transverse displacement. The linear strain–displacement relationships were introduced to describe the deformations of the plate. The strain components , and , at an arbitrary point of the plate, are given for small deformation as [5]
, 
(2ac) 
, 

. 
Based on Hooke's law, the stressdisplacement components in the host plate are expressed as follows:
, , , 
(3ac) 
where is Young's modulus, and is the shear modulus of the host plate. The constitutive relations in the piezoelectric layer can be written as
, , , 
(4ac) 
where and are the reduced material constraints of the piezoelectric medium for plane stress problems given by
(5ac) 
where , , , and are the moduli of elasticity under a constant electric field; and are the piezoelectric constants. The electric potential function is considered sinusoidal where it satisfied Maxwell's equations. Open circuit boundary conditions at the top and bottom of the piezoelectric layers are given by [4]
, 
(6) 
where is the electric potential across the piezoelectric middle layer. , , , and satisfy Maxwell's equations and open circuit boundary conditions (sensor conditions) that can be written as
, 
(7a, b) 
where is the electric displacement along the direction. The relationship between the electric field and electric displacement in the , and directions are given by
(8ac) 

, , . 
(9ac) 
In Eq. (9), and are reductions in the dielectric constants, and they can be written as
, . 
(10a, b) 
Eqs. (10a) and (10b) satisfy both Maxwell's equations and open circuit boundary conditions (sensor conditions), so that and are given by
. 
(11a, b) 
The expressions for the stress resultants and ( ) are [5]
(12ac) 
Based on the straindisplacement relations and stress distribution, the resultant bending moments, twisting moments, and shear forces in terms of and are obtained as follows [5]:
(13ae) 
where and .
The governing differential equations, based on the Kirchhoff plate theory in terms of the stress resultants using Eqs. (13ae), can be found as
(14) 
Note that all of the electrical variables must satisfy Maxwell’s equation, which requires that the divergence of the electric flux density vanishes at any point within the media. This condition can be satisfied approximately by enforcing the integration of the electric flux divergence across the thickness of the piezoelectric layers to be zero for any and , so the Maxwell’s equation after simplifying the result gives
(15) 
Considering Eqs. (14) and (15), the sixorder partial differential equation and the electric potential are obtained as follows
(16) 

(17) 

where
, , , , , , , 
(18ag) 

, , , . 
(19ad) 

After solving Eq. (17), the transverse displacement function of circular plates has taken the following form [4]
, 
(20) 
where is the transverse displacement amplitude in direction; is the natural frequency of the plate, and is the wavelength in direction. The transverse displacement amplitude is given by
(21) 
where
, , , 
(22ad) 
and
,
. 
(23ag) 
The clamped boundary conditions along the edges of the circular plate are as follows
. 
(24) 
Natural frequencies of piezoelectric coupled circular plates can be calculated by using the boundary conditions above as
(25) 
where
and is the sign of .
The equations of motion for a thin circular plate coupled with piezoelectric layers subjected to dynamic transverse loading (generated by the piezoelectric actuator) is given by
. 
(26) 
Using the modal superposition method
( ),
the forced vibration of the plate is expressed as
, 
(27) 
where is the modal shape function, and is the principal coordinate for the (m,n) modal of the plate. Substituting Eq. (17) into Eq. (27), the governing differential equations for forced vibration of the plate are obtained as
. 
(28) 
Multiplying Eq. (28) by , and integration of the summation of this equation over the plate area gives
, 
(29) 
where
, 
(30) 
. 
(31) 
Assuming zero initial conditions, Eq. (29) is solved as
. 
(32) 
Substituting Eq. (32) into ,
the transverse deflection of the plate due to forced vibration will be obtained as
. 
(33) 
Assuming free harmonic motion, the velocity distribution on the surface of the plate may be written as
, 
(34) 
where is the time, and . Here we assume that the circular plate radiator in flexural vibration is mounted on a flat rigid baffle of infinite extent. The coordinates shown in Figure 2 are positioned in the midplane surface of the plate. The acoustic pressure at the field point can be obtained by dividing the radiating surface of the flexural plate into infinitesimal elements ( ), where each element acts as a simple baffled source of strength having midpoint coordinates .
The distance between the midpoint of the infinitesimal elements ( ) and the observation point
( ) may be given by
, 
(35) 
where is the distance between the center of the spherical coordinates and the observation point.
The critical distance between the near and the far fields for a vibrating plate radiator as given by Khorshidi [12] may be approximated as
, 
(36) 
where is the wave length and is the frequency of the radiated sound that is also equal to the resonance frequency of the circular plate in the flexural vibration. Using the far field approximation, the distance in the acoustic pressure amplitude can be approximated as . The effect of distance on the phase of the acoustic pressure can be approximately expressed as
. 
(37) 
Based on the theory of a Rayleigh integral [13], the total far field acoustic pressure is
(38) 
Figure 2. A circular plate and its field point in the spherical coordinate system.
where , , and are the wave number, density of the medium, and speed of sound in the medium, respectively. Finally, the far field sound power of the radiator can be calculated using the following equation:
. 
(39) 
In this section, to ensure the correctness of computer programs, results were compared with the results obtained by Lee and Singh [18]. Therefore, the simple annular plate without a piezoelectric layer was considered. Material and geometry characteristics of the plate appear in Table 1.
It is assumed that the boundary conditions at r = b are clamped and at r = a are free. Therefore, the boundary conditions are expressed as

(40ad) 
Comparison between natural frequencies obtained by the method are presented in this paper, along with the exact solution, while the results that were obtained by Lee and Singh [18] are presented in Table 2. As shown in this table, there was very good agreement between the results obtained from the current method and the Lee and Singh results. After obtaining the natural frequencies and vibration modes, we can now obtain the sound radiation of the plate. Sound powers for two modes are calculated and compared with Lee and Singh [18] in Table 3. There are two reasons for errors in the results. One reason can be attributed to the methods used to solve the problem.
Table 1. Material and geometry characteristics of the plate. 

Outer radius (mm) 
151.5 
Inner radius b (mm) 
82.5 
Radii ratio 
0.54 
Thickness (mm) 
31.5 
Thickness ratio (h/a) 
0.21 
Mass density (kg/m^{3}) 
7905.9 
Young’s modulus (GPa) 
218 
Poisson ratio 
0.305 
Table 2. Frequencies of a thin annular plate with fixedfree boundary conditions. 

Diff. (%) 
Mode sequence (m,n) 

Lee and Singh [18] 
Present 

0.24 
5.465 
5.478 
(0,0) 
0.21 
5.580 
5.592 
(0,1) 
0.098 
6.091 
6.097 
(0,2) 
0 
7.360 
7.360 
(0,3) 
In reference vibration mode shape, shapes are approximated by a polynomial, whereas in the present study an exact solution was used to obtain the vibration mode shape. Another reason is related to plate thickness and the relationship between the frequency of the classical and Mindlin plate theories. When the thicknesstoradii ratio is more than 1/20, the classical thin plate theory does not provide accurate results.
After ensuring the accuracy of the formulation in this section, numerical results were calculated, according to the developed exact solution for the sound radiation of the vibrating circular plates coupled with piezoelectric layers. Calculations have been performed for the isotropic plate (host plate), using a Poisson ratio of 0.3 and a PZT4 piezoelectric layer by using commercial software, Mathematica (version 7). The results are presented in tabular and graphical forms for simplysupported and clamped boundary conditions, and plate and piezoelectric layer parameters. Material properties of the plate and piezoelectric layers are defined as listed in Table 4.
Table 3. Comparison of modal acoustic powers for two outofplane modes with fixedfree boundaries. 

Diff. (%) 
Method 
Mode sequence (m,n) 

13.1% 
84.2 
6.097 
Present 
(0,2) 
73.2 
4.85 
Lee and Singh [18] 

15.1% 
86.8 
7.360 
Present 
(0,3) 
73.7 
5.71 
Lee and Singh [18] 

Table 4. Material properties and dimensions of the plate and piezoelectric layer [5]. 

PZT4 
Host plate 
Material Property 

 
200 

132 
 

71 
 

115 
 

73 


 
0.3 
Poisson's ratio 

7500 
7800 
Mass density (kg/m^{3}) 

4.1 
 

14.1 
 

10.5 
 

7.124 
 

5.841 
 

In Figure 3, variation of sound pressure levels of the radiated plate coupled with a piezoelectric layer versus frequency are presented for the distance between the center of the spherical coordinates and the observation point R = 1.2 m, φ = 0.6 m, and ψ = 90. Acoustic frequency response function P (ω) for radius R = 1.2 m from the center of the plate and φ = 0, ψ = 90 are shown in Figure 3. The unit impulsive force was applied to the plate. As shown in Figure 3, peaks in this figure represent resonant frequencies (the undamped natural frequency of the plate is equal to the input excitation frequency). The natural frequencies of the plate were obtained using the graph acquired from the acoustic testing.
Twodimensional directivity patterns [12, 13] of the forced vibrations of a thin circular plate with a piezoelectric layer for the first four natural frequencies are shown in Figure 4. In this figure, the relationship between sound pressure distribution and the spherical coordinate's ψ is described. Results were obtained for a constant spherical coordinate's angle φ and ψ , which varies from π/2 to π/2. As shown in Figure 4, the sound pressure distribution for the mode (m = 0, n = 0) is uniform. The directivity pattern for the mode (m = 0, n = 1,2) and (m = 1, n = 0) dipped at ψ = 0, because ψ = 0 is the symmetry line of the sound pressure. In a mode shape that has a line of symmetry, this line affects the sound pressure and causes a dip.
6.1. Effect of Boundary Conditions on Sound Power
Acoustic sound power for the first four modes with clamped and simply supported boundary conditions are presented in Table 5. From the results displayed in Table 5, the sound powers of the vibrating plate under a clamped boundary condition were more than the sound powers of the vibrating plate under a simply supported boundary condition.
Figure 3. Variation of sound pressure with respect to frequanct (Acoustic frequency response function) for radius R = 1.2 m, φ = 0 , ψ = 90, given unit point force excitation (simply supported BCs).
(a)
(b)
(c)
(d)
Figure 4. Two dimensional directivity pattern: (a) m = 0, n = 0, 𝜔 = 902.5 rad/s; (b) m = 0, n = 1, 𝜔 = 1878.2 rad/s; (c) m = 0, n = 2, 𝜔 = 3081.13 rad/s; (d) m = 1, n = 0, 𝜔 = 3513.43 rad/s.
6.2. Effect of Piezoelectric Thickness on Sound Power
The influence of piezoelectric thickness, hp/2h, on sound powers of the radiated, vibrating simply supported circular plated coupled with piezoelectric layers is examined in Table 6. The results in Table 6 are presented for the thickness ratio, hp/2h, which varied from 1/12 to 1/5. As the piezoelectric thickness increased, the sound power decreased. This phenomenon occurred because the transverse displacement and the transverse velocity of the plate decreased with increased piezoelectric thickness.
6.3. Effect of the Piezoelectric Layer on the Sound Power
In this section, to show the effects of the piezoelectric layer on the sound radiation parameters, the sound power of a radiated vibrating clamped circular plate, with and without piezoelectric layers, for the first three vibrational modes were compared and presented in Table 7. According to the results with the piezoelectric layers, natural frequencies increased about 3.7%, and sound powers increased about 9%12% for the mode sequence (0,0), (0,1) and (0,2). Therefore, considering the inertia and stiffness of the piezoelectric layer, the vibration and sound transmission were required to obtain an accurate model for analysis.
Table 5. Modal acoustic powers for rgw ﬁrst four outofplane modes with simply supported and clamped boundary conditions. 

Mode sequence 
Boundary conditions 

78.88 
69.34 
(0,0) 
Simply supported 
87.68 
195.37 
(0,1) 

93.94 
360.08 
(0,2) 

73.76 
417.83 
(1,0) 

86.6 
143.63 
(0,0) 
Clamped 
90.54 
298.92 
(0,1) 

94.07 
490.37 
(0,2) 

75.58 
559.18 
(1,0) 
Table 6. Comparison between the sound powers of a simply supported circular plate coupled with piezoelectric layers with different piezoelectric thicknesses. 

Mode sequence 

(1,0) 
(0,2) 
(0,1) 
(0,0) 

73.85 
94.05 
87.75 
83.35 
1/12 
73.76 
93.94 
87.68 
78.88 
1/10 
73.68 
93.82 
87.61 
78.83 
1/8 
73.53 
93.71 
87.53 
75.78 
1/5 
Table 7. Effect of piezoelectric layers on the natural frequencies and sound powers 

Diff. (%) 
Plate with Piezo 
Plate without Piezo 
Parameter 
Mode sequence 
3.77 
902.494 
869.671 
(0,0) 

12.75 
105.021 
93.147 

3.72 
1878.2 
1810.79 
(0,1) 

10.29 
110.54 
100.23 

3.75 
3081.13 
2969.65 
(0,2) 

9.83 
113.96 
103.76 
In this study the classical plate theory (CPT) was used to investigate the transverse sound radiation of vibrating circular plates coupled with piezoelectric layers subjected to transverse external force. The exact closedform solution was obtained for the system with simply supported and clamped boundary conditions. The potential distribution function was earned by satisfying Maxwell's equation and boundary conditions. The frequencies, sound pressures, and sound powers of the system were presented in both tabular and graphical forms.
The significant advantages of the proposed closedform vibroacoustic equations consist of the following:
Based on comparison with previously published results, the accuracy of the present results were validated.
Acknowledgements
The authors gratefully acknowledge the funding by Arak University, under grant No 90/11542.
References
[1] Wang Q Quek ST, Sun CT, Liu X. Analysis of piezoelectric coupled circular plate. Smart Mater Struct 2001; 10: 229239.
[2] Hagood NW, McFarland AJ, Modeling of a piezoelectric rotary ultrasonic motor. IEEE Transactions on Ultrason, Ferroelectrics Freq Control 1995; 42: 210224.
[3] HeyligerPR, Ramirez G. Free vibration of laminated piezoelectric plates and discs. J Sound Vib 2000; 229: 935956.
[4] HosseiniHashemi Sh, Khorshidi K, Es’haghi M, Fadaee M, Karimi M. On the effects of coupling between inplane and outofplane vibrating modes of smart functionally graded circular/annular plates. Appl Math Modell 2012; 36: 1132–1147.
[5] Khorshidi K, Rezaei E, Ghadimi AA, Pagoli M. Active vibration control of circular plates coupled with piezoelectric layers excited by plane sound wave. Appl Math Modell 2015; 39(3–4): 1217–1228
[6] Rayleigh L. The Theory of Sound, Second edition, Reprinted by Dover, New York, 1945.
[7] Lomas NS, Hayek SI. Vibration and acoustic radiation of elastically supported rectangular plates. J Sound Vib 1977; 2: 1–25.
[8] Cremer L, Heckl M. StructureBorne Sound. Second edition, SpringerVerlag, New York, 1987.
[9] Berry A Guyader J, Nicolas, J. A general formulation for the sound radiation from rectangular, baffled plates with arbitrary boundary conditions. J Acoust Soc Am 1990; 88: 2792–2802.
[10] Yoo JW. Study on the general characteristics of the sound radiation of a rectangular plate with different boundary edge conditions. J Mech Sci Technol 2010; 24: 1111–1118.
[11] Rdzanek WP. The sound power of an individual mode of a clampedfree annular plate. J Sound Vib 2003; 261: 775–790.
[12] Khorshidi K. Vibroacoustic analysis of Mindlin rectangular plates resting on an elastic foundation. Scientia Iranica A 2011; 18: 45–52.
[13] HosseiniHashemi S, Khorshidi K, Rokni Damavandi Taher H. Exact acoustical analysis of vibrating rectangular plates with two opposite edges simply supported via Mindlin plate theory. J Sound Vib 2009; 322: 883–900.
[14] Hongqiu L, Guoping C, Linyan X. Structuralacoustic coupling and external sound pressure of a plateended cylindrical shell based on analytic method, Proc 2nd Int Conf Comput Eng Technol 2010; 5: 51–55.
[15] Zhou L, Zheng H, Hung KC. Sound radiation from a thin infinite plate in contact with a layered inhomogeneous fluid. Appl Acoust 2002; 63: 11771192.
[16] Frank F, Paolo G. Sound Radiation by Vibrating Structures. Second edition, Elsevier, Academic Press, 2007; 135241.
[17] Zhang X, Li WL. A unified approach for predicting sound radiation from baffled rectangular plates with arbitrary boundary conditions. J Sound Vib 2010; 329: 5307–5320.
[18] Lee M, Singh R. Analytical formulations for annular disk sound radiation using structural modes. J Acoust Soc Am 1994; 95: 3311–3323.
[19] Lee H, Singh R. Acoustic radiation from outofplane modes of an annular disk using thin and thick plate theories. J Sound Vib 2005; 282: 313–339.
[20] Lee H, Singh R. Determination of sound radiation from a simpliﬁed diskbrake rotor by a semianalytical method. Noise Control Eng J 2004; 52: 225–239.