Abstract

Radial vibration of the circular plate is presented using wave propagation approach and classical method containing Bessel solution and Hankel solution for calculating the natural frequency theoretically. In cylindrical coordinate system, in order to obtain natural frequency, propagation and reflection matrices are deduced at the boundaries of free-free, fixed-fixed, and fixed-free using wave propagation approach. Furthermore, radial phononic crystal is constructed by connecting two materials periodically for the analysis of band phenomenon. Also, Finite Element Simulation (FEM) is adopted to verify the theoretical results. Finally, the radial and piezoelectric effects on the band are also discussed.

1. Introduction

Wave propagation approach is a useful method on the aspects of calculating natural frequency for the analysis of elastic structures, such as rod, Euler beam, curved beam, and plate. It benefits us with a better understanding to analyze wave when propagating inside structures. In fact, as early as 1984, Mace analyzed the behaviors of wave by using wave propagation approach through dividing them into propagation attenuation matrices. And the reflection matrices under three boundary conditions were deduced, which built a theoretical foundation for wave method [1]. Meanwhile, the characteristics of propagation and reflection of Timoshenko beam were also discussed with the discontinuous behavior. And the frequency curves were obtained [2]. Harland et al. studied the motion of waves in tunable fluid-filled beams [3]. Based on Love theory, Kang et al. deduced the propagation matrices, reflection matrices, and coordinated matrices in curved beams [4]. Different from Kang et al., Lee made a detailed research on the behavior of waves of curved beam based on Flugge’s theory [5, 6]. Based on the high-order plate theory, Wang and Rose analyzed the properties of transient wave in an inhomogeneous beam. And the results were also verified by the experiment [7]. Tan and Kang considered the wave reflection and transmission in an axially strained, rotating Timoshenko shaft. Meanwhile, the influence of continuous condition and cross-section on natural frequency was also discussed [8]. Li and Gan analyzed the free vibration of cylindrical shell by using wave propagation approach based on Flugge’s thin shell theory [9, 10]. Also, Zhang calculated the frequencies of cross-ply laminated composite cylindrical shells and submerged cylindrical shells using wave approach based on Love’s theory [11, 12]. Using wave approach, the above scholars have done lots of analysis, which exhibit a unique advantage for the analysis of structural vibration.

In recent years, people have paid great attention to structure which consisted of two materials arraying periodically. At present, the research has been extended into radial phononic crystal structures. Shakeri et al. analyzed the behaviors of radial wave in functional graded radial periodical structures [13]. Torrent and Carbonell studied the propagation characteristic of acoustic in detail when propagating in radial periodical structure and found that it also had a low frequency band [14, 15]. After constructing and analyzing the two-dimensional radial phononic crystal, Xu et al. found that these structures had low frequency band that the starting frequency was zero [16]. Ma et al. studied the band characteristics of Lamb wave in two-layer radial phononic crystal plate, and it was found that crystal gliding in radial direction could open the lowest band gaps [17]. The above researches indicate that radial phononic crystal plate has the low frequency band which brings an attractive prospect of engineering application. For example, the natural frequency can be designed into the range of band gaps to isolate the vibration. Also, these structures can be regarded as a filter. Additionally, Shu et al. researched torsional vibration problem of generalized phononic crystal, finding that low frequency band was determined by radius. However, the high frequency band was caused by the periodicity of circular plate [18]. And Shu’s research shows that generalized phononic crystal can generate a huge band which almost covers the whole frequency range, providing a practical application for the rotating parts, such as gear. In the previous research, Shu et al. have studied the band gaps of flexural wave in the radial phononic crystal. Meanwhile, numerical results are also verified by FEM [19]. Herein, this paper is mainly focused on radial vibration of composite structures. In fact, other works on radial vibration of elastic solids under various loadings can also be found [2022].

The paper is organized into five sections. In Section 1, a brief introduction is given. In Section 2, classic Bessel solution and Hankel solution of radial vibration are given. Meanwhile, propagation and reflection matrices at three boundary conditions are derived. The results of natural frequency obtained by wave approach are also compared with the classical method. In Section 3, the band phenomenon is analyzed and the effect of parameters is discussed. Section 4 is the conclusion.

2. Theoretical Analysis for Natural Frequency

2.1. Solution of Radial Vibration Equation

Single circular plate is shown in Figure 1. Young modulus, density, and Poisson’s coefficient are . The thickness of piezoelectric ceramics is . The inner and external radius are and , respectively. With respect to the axisymmetric circular plates, radial equation of piezoelectric materials can be given by [23]

The expression of strain and radial displacement can be written as

Since piezoelectric material direction of polarization is along z-axle, piezoelectric equation can be expressed aswhere is radial displacement, and are radial stress and tangential stress, and are radial strain and tangential strain, and is radius.

Substituting (2)–(4) into (1) giveswhere is the square of radial wave velocity.

After simplifying (5), it can be obtained that

Thus, general solution of radial displacement can be written aswhere is the wave number. is the Bessel function of first kind; is the Bessel function of second kind. And the constants and can be determined by boundary conditions.

Therefore, radial stress can be obtained as

2.2. Classical Method

When both ends of circular plate are fixed, it has

Substituting (7) into (9) gives

When both ends of circular plate are free, it has

Substituting (8) into (11) gives

Similarly, when inner surface is fixed and outer surface is free, it has

Substituting (7) and (8) into (13) gives

After solving roots of (10), (12), and (14), the natural frequencies of circular plate structure can be obtained with the aforementioned boundary conditions.

2.3. Wave Propagation Method

One of the solutions of circular plate has been solved as shown in (7). And the other solution can be written aswhere and are the Hankel function of first and second kinds, which can be defined as

Considering (15), forward-traveling and negative-traveling waves can be written as

Then, (15) reduces to

When both ends of circular plate are fixed, it gives

Thus, the reflection matrices of both boundaries are

Substituting (2), (15), and (16) into (8), radial stress can be expressed as

When both ends of circular plate are free, it has

Thus, reflection matrices of both boundaries are

Forward-traveling and negative-traveling waves can be shown clearly in Figure 1. And the expression can be written as follows:

Considering (17), forward-traveling and negative-traveling waves can be obtained as

Similarly, the reflection matrices are (20a) and (24) for case of inner surface fixed and outer surface free.

2.4. Analysis and Discussion

After combining propagation and reflection matrices derived above, the natural characteristics of single circular plate can be analyzed effectively. Hence, it gives

And it yields

RESIN is selected as the material of circular plate. Material and structural parameters are tabulated in Table 1.

Substitute parameters listed in Table 1 into (10), (12), (14), and (28), and the natural frequency can be computed. Herein, red and black curves obtained by wave approach are real and imaginary wave solution, respectively. Also, it is important to note that the intersection points of the curves are natural frequency. After solving the inherent characteristic of circular plate structure, the characteristic curves calculated by classical method and wave approach are presented in Figure 2. It can be seen that these natural frequencies intersect at the same point in x-axis, which verifies the correctness of the numerical calculation.

3. Radial Vibration Analysis of Phononic Crystal

3.1. Model and Transfer Matrix

Figure 3 is the schematic diagram of radial phononic crystal circular plate. Herein, the white region represents RESIN materials and the gray region represents piezoelectric ceramic PZT4. The inner and outer radius are and , respectively. The radial span ratios are denoted with and . And represents the thickness.

Considering an eight-layer circular plate, it corresponds to the model . The interfaces between RESIN and PZT4 material are (see Figure 3). With regard to RESIN, Young’s modulus, density, and Poisson coefficient are . For the PZT4, Young’s modulus, density, and Poisson coefficient are . When radial wave propagates in RESIN layers 1, 3, 5, and 7, the velocity is and the wavenumber is . The velocity is and the wavenumber is when wave propagates in PZT4 layers 2, 4, 6, and 8.

With regard to radial vibration, there are two parameters which satisfy the continuous conditions at the interface of different materials. Namely, they are the radial displacement and radial stress. So transfer matrices of radial vibration can be derived as follows.

At the interface b (), it gives

Thus, it can be arranged in a matrices form

At the interface c (), it gives

Similarly, one has

Combining (30) and (32), one haswhere is the transfer matrices of radial vibration between the 3rd and the 1st ring.

As shown in Figure 3, replace and (in (30) and (32)) with and . Similarly with , we can obtain which is the transfer matrix between the fifth ring and the third ring.

Similarly, replace and with and , and the matrix can be obtained.

At last, submitting into (30), we can get . Therefore, with regard to the radial phononic crystal, one has

3.2. Numerical Result and Discussion

The radial phononic crystal circular plates are composed of two materials, namely, RESIN and piezoelectric ceramic PZT4. Material parameters and structural parameters are consistent with Table 2.

Additionally, piezoelectric parameters are set as  m2/N, ,  c/m, and  c/N.

Taking advantage of transfer matrices and boundary conditions, the transmission characteristics are calculated using MATLAB numerically. Figure 4 is the transmission curve of radial vibration with four periodic structures which is composited by single material RESIN and RESIN/PZT4.

Figure 4 is the comparison of numerical and FEM results. It indicates that there is no attenuation when radial wave propagates in a single material, while, for the phononic crystal, the transmission curve shows a big vibration attenuation band after introducing the periodicity. However, when radial wave propagates in the single material, there is no band gap which attenuates weakly and only owns −4 dB. Furthermore, ANSYS 14.5 is used to verify the radial propagation characteristics. Using harmonic analysis and keeping the boundary conditions of outer surface free, radial displacement is excited on the inner surface which is shown in Figure 4(b). Comparing numerical results with FEM results, it can be found that the bands are in good agreement with the aspects of position, width, and amplitude attenuation and other aspects, which indicates the correctness of theoretical calculations.

3.3. Effects of Structural Parameters

In this section, the influence of piezoelectric effect on radial wave has been discussed in detail. Furthermore, the influences of piezoelectric parameter and polar direction on the band have also been studied in case of open and short circuit for piezoelectric material.

3.3.1. Open and Short Circuit

With regard to piezoelectric material, the D-type equation of radial vibration can be expressed as

It corresponds to short circuit situation when field intensity is . Therefore, (35) can be simplified asThen, it haswhere denotes Poisson’s coefficient in the condition of short circuit. denotes Young modulus.

Piezoelectric material is in the open circuit when the charge is . Thus, (35) can be simplified as

Then, it giveswhere denotes Poisson’s coefficient in the condition of open circuit. denotes Young modulus.

By using transfer matrices derived in Section 2.1, vibration characteristics curves have been calculated for the case of short and open circuit. As depicted in Figure 5, the selected frequency band is 0–10 kHz. It indicates that the band changes obviously in the high frequency because of the piezoelectric effect. However, the short circuit curve moves toward low frequency slightly. From the above theoretical analysis by comparing (37) and (39), it can be observed that the elasticity modulus and Poisson’s coefficient of piezoelectric materials are different when they are in short or open circuit. For the case of open circuit, elasticity modulus is larger. Meanwhile, Poisson’s coefficient also changes. It is important to note that these two parameters can influence the value of and directly, which means that radial wavenumber changes. Thus, it is the natural mechanism of band moving.

3.3.2. Piezoelectric Parameters and Inner Radius

As to the open circuit case, the effects of electrical parameter and inner radius on radial wave band are depicted in Figure 6. Herein, Figure 6(a) describes the effect of piezoelectric parameter . However, Figure 6(b) displays the effect of inner radius .

From (39), it can be easily observed that determines the values of and which can influence the band position and attenuation. Figure 6(a) indicates that radial band has been expanded with the increase of (here, the negative sign only represents direction). However, the change is not obvious. Figure 6(b) displays that inner radius can influence the low frequency greatly. It is also clear from the figure that the band attenuates obviously at low frequency. However, the effect on high frequency is much smaller.

4. Conclusion

This paper presents the classical method and wave propagation approach to obtain the natural frequency. Based on radial vibration equation, propagation and reflection matrices are derived. It can be found that natural frequency of the single circular plate can be calculated conveniently by using wave propagation approach. Then, the results obtained by using wave approach are compared with the classical method at the boundary conditions of free-free, fixed-fixed, and fixed-free. It can be seen that the results calculated by these two methods coincide with each other.

Additionally, the band phenomenon of radial phonon crystal circular plate is also analyzed for a deeper understanding of radial vibration. FEM is used to verify the theoretical results. It can be found that the inner radius of circular plate has great influence on the low frequency. Piezoelectric coefficient, open circuit, and short circuit make the elastic modulus and Poisson’s ratio change. Then the radial wavenumber also changes, which is the natural mechanism of band removing for the case of open and short circuits.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work is financially supported by the National Natural Science Foundation of China (no. 51375105) and the Natural Science Foundation of Heilongjiang Province of China (E201418).