Shock and Vibration
Volume 2018 (2018), Article ID 6181204, 14 pages
https://doi.org/10.1155/2018/6181204

Research Article

Free Vibration Analysis of Rings via Wave Approach

Wang Zhipeng, Liu Wei, Yuan Yunbo, Shuai Zhijun, Guo Yibin, and Wang Donghua

College of Power and Energy Engineering, Harbin Engineering University, Harbin 150001, China

Correspondence should be addressed to Shuai Zhijun; shuaizhijun@hrbeu.edu.cn

Received 1 February 2018; Revised 28 March 2018; Accepted 2 April 2018; Published 14 May 2018

Academic Editor: Tai Thai

Copyright © 2018 Wang Zhipeng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Free vibration of rings is presented via wave approach theoretically. Firstly, based on the solutions of out-of-plane vibration, propagation, reflection, and coordination matrices are derived for the case of a fixed boundary at inner surface and a free boundary at outer surface. Then, assembling these matrices, characteristic equation of natural frequency is obtained. Wave approach is employed to study the free vibration of these ring structures. Natural frequencies calculated by wave approach are compared with those obtained by classical method and Finite Element Method (FEM). Afterwards natural frequencies of four type boundaries are calculated. Transverse vibration transmissibility of rings propagating from outer to inner and from inner to outer is investigated. Finally, the effects of structural and material parameters on free vibration are discussed in detail.

1. Introduction

Noise and vibration problem of rings has been a hotspot that makes it attract lots of attention because of the wide applications, such as gear transmission system and aircraft structures. Since most of components within these structures can be regarded as a simple model of ring structure, dynamic properties of these structures in recent years have been the subject of present studies. Based on the Stodola–Vianello method, Gutierrez et al. [1] considered the free vibration of annular membranes with continuously variable density by using Rayleigh method, differential quadrature method, and finite elements simulation. By adopting the dynamic stiffness method, Jabareen and Eisenberger [2] investigated the natural characteristics of the nonhomogeneous annular structure. Wang [3] discussed the natural frequency with the boundaries of fixed, free, and simply supported. And the researches indicated clearly that complex modes would switch correspondingly with the radius increasing gradually. Based on the first-order shear deformation theory, Roshan and Rashmi [4] analyzed the free vibration of axisymmetric sandwich circular plate with relatively stiff thickness. Oveisi and Shakeri [5] constructed a sandwich composite circular plate containing piezoelectric material, finding that the feedback control gain had an effective control for suppressing the transverse vibration. Hosseini-Hashemi et al. [6] presented the free vibration of functionally graded circular plate with stepped thickness via Mindlin plate theory. By employing harmonic differential quadrature method for obtaining numerical solution of circular plate, Civalek and Uelker [7] made a further analysis for the behavior of the boundary conditions of fixed and simply supported. Moreover, the results obtained by finite difference method were compared with those calculated by harmonic differential quadrature method, whose feasibility is verified by Bakhshi Khaniki and Hosseini-Hashemi [8]. Liu et al. [9] construct the composite thin annular plate by wave approach. However, composite structures with multilayer are not studied. Different boundaries and parameter effects are also not analyzed.

Additionally, in terms of the vibration of structures, there is an alternative method called “wave approach” which is very efficient and widely used for calculating the natural frequency of structures, such as beams, plates, rings, and periodic structure by describing waves in the matrix form. For example, as early as 1984, Mace [10] applied wave approach to analyze the wave behaviors in Euler beam. By dividing the waves into propagation and attenuation matrices, he derived the reflection matrices under three boundary conditions, which established a theoretical foundation for wave approach. Adopting wave approach, Mei [11] analyzed the flexural vibration with added mass for Timoshenko beam, and the influence of lumped mass on natural frequency is also discussed in detail. Kang et al. [12] divided the real and imaginary parts of wave solutions of curved beam into four cases, and they calculated the natural frequencies by combining propagation, transmission and reflection matrices. Lee et al. [13, 14] considered the power flow when wave propagated in curved beam. Furthermore, they applied the Flugge theory to analyze the free vibration of a single curved beam, and their result was compared with Kang et al., which verified the correctness of the numerical results. Huang et al. [15] investigated the free vibration of planar rotating rings. The effect of cross section on natural frequency was also discussed. Bahrami and Teimourian [16] studied the free vibration of composite plates consisting of two layers, and they also made a comparison between classical results and wave propagation results. Tan and Kang [17] concentrated on the free vibration of rotating Timoshenko shaft with axial force and discussed the effect of continuous condition and cross section on natural frequencies. From the wave point of view, Bahrami and Teimourian [18] analyzed the free vibration of nanobeams for the first time. Ilkhani et al. [19] studied the free vibration of thin rectangular plate. It should be noted that the above scholars have done lots of studies for free vibration of structures by using wave approach, while few reports for the analysis of natural frequency for transverse vibration of rings can be found. In fact, it is well known that the nature of vibration is the propagation of waves. Analyzing free vibration in terms of wave propagation and attenuating can have a better understanding for us. Moreover, one advantage of using wave approach to analyze the free vibration is its conciseness of matrices that makes the natural frequencies be calculated easily. Wave approach is a strong tool for studying the behavior of wave transmission and reflection in waveguides, providing a practical engineering application such as filters.

The emphasis of this paper is focused on the free vibration of rings. This paper is organized into five parts. Section 1 is introduction. In Section 2, propagation, coordination, and reflection matrices are deduced in forms of matrix. In addition, the characteristic equation of natural frequency is obtained using classical method and wave approach. In Section 3, natural frequencies of rings are calculated by combining these matrices. Meanwhile, vibration transmissibility of rings propagating from outer to inner and from inner to outer is obtained. In Section 4, the influence of structural and material parameters on natural frequencies is discussed. Section 5 is the conclusion.

2. Theoretical Analysis

2.1. Classical Method for Free Vibration

2.1.1. Solution of Transverse Vibration

Consider sandwich rings consisting of two different materials depicted in Figure 1. Adhesive can be employed for connecting the rings. Same material is selected for the first and third layers. The other material is selected for the second layer. Radius of the first and third layers is and . Radius of the intermediate layers is and . Radial span of the first and third layers is . Radial span of the second layer is . is bending deflection. is thickness. At the boundaries of and , positive–going and negative–going wave vectors are , , , , , , , . Also, considering another two boundaries at and , positive–going and negative–going wave vectors are , , , . In cylindrical coordinates, the radius is assumed to be large enough compared to thickness which means that it satisfies the small deformation theory. Transverse solution is given by Wang [3]:

where , , , are constants which are determined by boundaries. and are Bessel functions of first and second kinds, respectively. and are modified Bessel functions of first and second kinds. is wave number, and is stiffness.

Figure 1: Composite rings.

2.1.2. Solution of Classical Bessel Method

With regard to rings subjected to bending excitation, expression of transverse displacement, rotational angle, shear force, and bending moment within the first and third layers can be written as

Applying fixed boundary condition at obtains

where .

Free boundary condition is selected at ; then

where .

In order to obtain the natural frequencies, substituting (A.9) into (5) and combining (3)-(4), it reduces to

where the specific theoretical derivation of in (6) is presented in the Appendix. And each element is defined as

Therefore, (6) can be written as a determinant:

where (8) is the characteristic equation of natural frequency. By searching the root, natural frequency of rings can be calculated with a fixed boundary at inner surface and a free boundary at outer surface.

2.1.3. Solution of Classical Hankel Method

The solution is obtained in (1). However, it also can be expressed in a Hankel form:

where and are the Hankel functions of second and first kinds, respectively. They can be defined as

Similarly, expression of parameters within the first and third layers can be written as

Natural frequencies of transverse vibration can be calculated using classical Hankel method. Characteristic equation of natural frequency can be deuced like the process of (3)–(8). In order to avoid repeating, herein, it is ignored.

2.2. Wave Approach for Free Vibration

In this section, the solution is presented in terms of cylindrical waves for this ring. Meanwhile, positive–going propagation, negative–going propagation, coordination, and reflection matrices are also deduced. By combining these matrices, natural frequencies are calculated using wave approach.

2.2.1. Propagation Matrices

Wave propagates along the positive–going and negative–going directions when propagating within structures, as is shown in Figure 1. Waves will not propagate at the boundaries but only can be reflected. Moreover, parameters are continuous for the connection. In recent years, many researchers describe the waves in the matrix forms [8–17].

By considering (11), positive–going waves can be described as

These wave vectors are related by

Substituting matrices (15a), (15b), (16a), (16b), (17a), (17b) into (18a)–(18c), positive–going propagation matrices are obtained as

Similarly, negative–going waves can be rewritten as

These wave vectors are related by

Substituting matrices (20a), (20b), (21a), (21b), (22a), (22b) into (23a)–(23c), negative–going propagation matrices are obtained:

2.2.2. Reflection Matrices

Keeping the boundary condition of fixed, thus, displacements and rotational angle are taken as

The relationship of incident wave and reflected wave is related by

Substituting (25) into (26), the reflection matrices can be obtained as follows:

Keeping the boundary condition of free gives

The relationship of incident wave and reflected wave is

Substituting (28) into (29), reflection matrices are calculated as

2.2.3. Coordination Matrices

By imposing the geometric continuity at yields

Equations (31) can be rewritten as

According to the continuity at , shear force and bending moment are required that

Equations (33) can be written as

2.2.4. Characteristic Equation of Natural Frequency

Combining propagation matrices, reflection matrices, and coordination matrices derived in Section 2.2, natural frequencies of composite rings can be calculated smoothly. Figure 1 presents a clear description of incident and reflected waves. Thus, the wave matrices described by (18a)–(18c), (23a)–(23c), (26), (29), (32), and (34) are assembled as

In order to obtain the natural frequency, (35) can be rewritten in a matrix form

is a matrix of . If (36) has solution, it requires that

By solving the roots of characteristic equation (37), one can calculate the real and imaginary parts. It is important here to note that the natural frequencies can be found by searching the intersections in -axis.

3. Numerical Results and Discussion

In this section, free vibration of rings is calculated by using wave approach, and the results are also compared with those obtained by classical method. Material RESIN is selected for the first and third layers. Material STEEL is selected for the middle layers. Material and structural parameters are given in Table 1.

Table 1: Material and structural parameters.

Material parameters
Density (kg/m³)
Young modulus (Pa)
Poisson’s ratio

I (RESIN)
1180

0.3679

II (STEEL)
7780

0.3

Structural parameters

h

(mm)
45
125
1

Based on Bessel and Hankel solutions calculated by classical method theoretically, natural frequency curves are presented by solving characteristic equation (8) depicted in Figure 2. Furthermore, (37) is calculated using wave approach. It can be seen that the real and imaginary parts intersect at multiple points simultaneously in -axis. It is important, here, to note that the roots of the characteristic curves are natural frequencies when the values of longitudinal coordinates are zero.

Figure 2: Natural frequency obtained by classical method and wave approach.

In Figure 2, two different natural frequencies can be clearly presented in the range of 450–1500 Hz, that is, 1244.22 Hz and 1443.31 Hz. However, the values are very small in the range of 0–450 Hz. In order to find whether the values in this range also intersect at one point, three zoomed figures are drawn for the purpose of better illustration about the natural frequencies of characteristic curves which are described in Figure 3.

(a) 0–40 Hz

(b) 40–180 Hz

(c) 180–450 Hz

Figure 3: Characteristic curves in the range of 0–450 Hz.

Natural frequencies calculated by these two methods are compared. Modal analysis is carried out by FEM. The natural frequencies are presented in Table 2 from which it can be observed that the first five-order modes calculated by these three methods are in good agreement. Obviously, it also can be found that natural frequencies obtained by ANSYS software are larger than the results calculated by classic method and wave approach, which is mainly caused by the mesh and simplified solid model in FEM. However, these errors are within an acceptable range, which verifies the correctness of theoretical calculations. To assess the deformation of rings, Figure 4 is employed to describe the mode shape. It can be found that the maximum deformations of the first three mode shapes occur in the outermost surface. The fourth and fifth mode shapes appear in the innermost surface.

Table 2: Results calculated by classical method, wave approach, and FEM.

Method
1st mode
2nd mode
3rd mode
4th mode
5th mode

Classical Bessel
37.65 Hz
167.54 Hz
414.27 Hz
1244.22 Hz
1443.31 Hz

Classical Hankel
37.65 Hz
167.54 Hz
414.27 Hz
1244.22 Hz
1443.31 Hz

Wave approach
37.65 Hz
167.54 Hz
414.27 Hz
1244.22 Hz
1443.31 Hz

FEM
37.76 Hz
168.30 Hz
415.19 Hz
1247.90 Hz
1448.11 Hz

(a)

(b)

(c)

(d)

(e)

Figure 4: Mode shapes of natural frequencies. (a) First mode. (b) Second mode. (c) Third mode. (d) Fourth mode. (e) Fifth mode.

Adopting FEM method, the first five natural frequencies are calculated for four type boundaries, as is shown in Table 3. It shows that the first natural frequency is 37.76 Hz (Min) at the case of inner boundary fixed and outer boundary free. The first natural frequency is 143.40 Hz (Max) at the case of inner and outer boundaries both free.

Table 3: Comparison of free vibration by FEM for four type boundaries.

Different boundaries
1st mode
2nd mode
3rd mode
4th mode
5th mode

Inner free, outer free
143.40 Hz
334.10 Hz
569.55 Hz
1336.50 Hz
1845.16 Hz

Inner fixed, outer free
37.76 Hz
168.30 Hz
415.19 Hz
1247.90 Hz
1448.10 Hz

Inner free, outer fixed
70.44 Hz
321.96 Hz
572.94 Hz
1388.86 Hz
1848.67 Hz

Inner fixed, outer fixed
101.07 Hz
413.28 Hz
1272.37 Hz
1479.74 Hz
2927.38 Hz

Harmonic Response Analysis of rings is carried out by using ANSYS 14.5 software. RESIN is chosen for the first and third layer. The second layer is selected as STEEL. Element can be selected as Solid 45, which is shown in red and blue in Figure 5(a). Through loading transverse displacement onto the innermost layer and picking the transverse displacement onto the outermost layer, vibration transmissibility of rings propagating from inner to outer is obtained by using formula . Similarly, through loading transverse displacement onto the outermost layer and picking the transverse displacement onto the innermost layer, vibration transmissibility propagating from outer to inner is obtained by using formula .

(a) The meshing mode

(b) Vibration response

Figure 5

Figure 5(b) indicates that there is no vibration attenuation in the range of 0–1500 Hz when transverse vibration propagates from outer to inner. Also, four resonance frequencies appear, namely, 70.44 Hz, 321.96 Hz, 572.94 Hz, 1388.86 Hz, which coincide with the first four-order natural frequencies in Table 3 at the case of innermost layer free and outermost layer fixed. Compared with the case of vibration propagation from outer to inner, there is vibration attenuation when vibration propagates from inner to outer. In addition, five resonance frequencies also appear, namely, 37.76 Hz, 168.30 Hz, 415.19 Hz, 1247.9 Hz, and 1448.1 Hz, which coincide with the results obtained by wave approach, classical Hankel, and classical Bessel methods shown in Table 2.

4. Effects of Structural and Material Parameters

4.1. Structural Parameters

The effects of structural parameters such as thickness, inner radius, and radial span are investigated in Figure 6. Adopting single variable principle, herein, only change one parameter. Figure 6(a) shows clearly that, with thickness increasing, the first modes change from 37.76 Hz to 188.15 Hz, and the remaining three modes increase obviously, which indicates that thickness has great effect on the first four natural frequencies. In fact, characteristic equation of natural frequency is determined by thickness, density, and elastic modulus, which is shown by the expression of wave number and stiffness . Therefore, thickness is used to adjust the natural frequency directly through varying wave number in (36).

(a) Thickness

(b) Inner radius

(c) Radial span

Figure 6: Effect of structural parameters.

From the wave number , it can be found that inner radius is not related to the natural frequency. Thus, inner radius almost has no effect on the natural frequency shown in Figure 6(b).

In Figure 6(c), there are five different types analyzed for the radial span ratios of RESIN and STEEL, that is,  :  : ,  :  : ,  :  : ,  :  : ,  :  : , respectively. When radial span is equal to 1, this means that the size of RESIN and STEEL is 1 : 1, namely, . For this case, the total size of composite ring is max, so the mode is min. Additionally, symmetrical five types cause the approximate symmetry of Figure 6(c). It also can be found that as radial span increases, natural frequencies appear as a similar trend, namely, decrease afterwards increase.

4.2. Material Parameters

Adopting single variable principle, density of middle material STEEL is replaced by the density of PMMA, Al, Pb, Ti. Similar with the study on effects of structural parameters, the effects of density and elastic modulus are studied for the case of keeping the material and structural parameters unchanged. Also, material parameters of PMMA, Al, Pb, and Ti are presented in Table 4.

Table 4: Material parameters.

Method
Density (kg/m³)
Young modulus (Pa)
Poisson’s ratio

PMMA
1062
0.32 × 10¹⁰
0.3333

Al
2799
7.21 × 10¹⁰
0.3451

Pb
11600
4.08 × 10¹⁰
0.3691

Ti
4540
11.7 × 10¹⁰
0.32

Figure 7(a) indicates that as density increases, the first mode decreases but not very obviously. However, the second, third, and fourth modes reduce significantly. Figure 7(b) shows that when elastic modulus increases gradually, the first mode increases but not significantly. The second, third, and fourth modes increase rapidly.

(a) Density

(b) Elastic modulus

Figure 7: Effect of material parameters.

5. Conclusion

This paper focuses on calculating natural frequency for rings via classical method and wave approach. Based on the solutions of transverse vibration, expression of rotational angle, shear force, and bending moment are obtained. Wave propagation matrices within structure, coordination matrices between the two materials, and reflection matrices at the boundary conditions are also deduced. Additionally, characteristic equation of natural frequencies is obtained by assembling these wave matrices. The real and imaginary parts calculated by wave approach intersect at the same point with the results obtained by classical method, which verifies the correctness of theoretical calculations.

A further analysis for the influence of different boundaries on natural frequencies is discussed. It can be found that the first natural frequency is Min 37.76 Hz at the case of inner boundary fixed and outer boundary free. In addition, it also shows that there exists vibration attenuation when vibration propagates from inner to outer. However, there is no vibration attenuation when vibration propagates from outer to inner. Structural and material parameters have strong sensitivity for the free vibration.

Finally, the behavior of wave propagation is studied in detail which is of great significance to the design of natural frequency for the vibration analysis of rotating rings and shaft systems.

Appendix

Derivation of the Transfer Matrix

Due to the continuity at , the following is obtained:

Equation (A.1) can be organized as

where , and each element is defined as

Hence, (A.2) can be written as

Similarly, by imposing the geometric continuity at , the following is obtained:

Arranging (A.5) yields

and each element is defined as

Equation (A.6) can be simplified as

Combining (A.4) and (A.8) gives

where is the transfer matrix of flexural wave from inner to outer.

Conflicts of Interest

There are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The research was funded by Heilongjiang Province Funds for Distinguished Young Scientists (Grant no. JC 201405), China Postdoctoral Science Foundation (Grant no. 2015M581433), and Postdoctoral Science Foundation of Heilongjiang Province (Grant no. LBH-Z15038).

References

R. H. Gutierrez, P. A. A. Laura, D. V. Bambill, V. A. Jederlinic, and D. H. Hodges, “Axisymmetric vibrations of solid circular and annular membranes with continuously varying density,” Journal of Sound and Vibration, vol. 212, no. 4, pp. 611–622, 1998.

M. Jabareen and M. Eisenberger, “Free vibrations of non-homogeneous circular and annular membranes,” Journal of Sound and Vibration, vol. 240, no. 3, pp. 409–429, 2001.

C. Y. Wang, “The vibration modes of concentrically supported free circular plates,” Journal of Sound and Vibration, vol. 333, no. 3, pp. 835–847, 2014.

L. Roshan and R. Rashmi, “On radially symmetric vibrations of circular sandwich plates of non-uniform thickness,” International Journal of Mechanical Sciences, vol. 99, article no. 2981, pp. 29–39, 2015.

A. Oveisi and R. Shakeri, “Robust reliable control in vibration suppression of sandwich circular plates,” Engineering Structures, vol. 116, pp. 1–11, 2016.

S. Hosseini-Hashemi, M. Derakhshani, and M. Fadaee, “An accurate mathematical study on the free vibration of stepped thickness circular/annular Mindlin functionally graded plates,” Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 37, no. 6, pp. 4147–4164, 2013.

Ö. Civalek and M. Uelker, “Harmonic differential quadrature (HDQ) for axisymmetric bending analysis of thin isotropic circular plates,” Structural Engineering and Mechanics, vol. 17, no. 1, pp. 1–14, 2004.

H. Bakhshi Khaniki and S. Hosseini-Hashemi, “Dynamic transverse vibration characteristics of nonuniform nonlocal strain gradient beams using the generalized differential quadrature method,” The European Physical Journal Plus, vol. 132, no. 11, article no. 500, 2017.

W. Liu, D. Wang, and T. Li, “Transverse vibration analysis of composite thin annular plate by wave approach,” Journal of Vibration and Control, p. 107754631773220, 2017.

B. R. Mace, “Wave reflection and transmission in beams,” Journal of Sound and Vibration, vol. 97, no. 2, pp. 237–246, 1984.

C. Mei, “Studying the effects of lumped end mass on vibrations of a Timoshenko beam using a wave-based approach,” Journal of Vibration and Control, vol. 18, no. 5, pp. 733–742, 2012.

B. Kang, C. H. Riedel, and C. A. Tan, “Free vibration analysis of planar curved beams by wave propagation,” Journal of Sound and Vibration, vol. 260, no. 1, pp. 19–44, 2003.

S.-K. Lee, B. R. Mace, and M. J. Brennan, “Wave propagation, reflection and transmission in curved beams,” Journal of Sound and Vibration, vol. 306, no. 3-5, pp. 636–656, 2007.

S. K. Lee, Wave Reflection, Transmission and Propagation in Structural Waveguides [Ph.D. thesis], Southampton University, 2006.

D. Huang, L. Tang, and R. Cao, “Free vibration analysis of planar rotating rings by wave propagation,” Journal of Sound and Vibration, vol. 332, no. 20, pp. 4979–4997, 2013.

A. Bahrami and A. Teimourian, “Free vibration analysis of composite, circular annular membranes using wave propagation approach,” Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 39, no. 16, pp. 4781–4796, 2015.

C. A. Tan and B. Kang, “Free vibration of axially loaded, rotating Timoshenko shaft systems by the wave-train closure principle,” International Journal of Solids and Structures, vol. 36, no. 26, pp. 4031–4049, 1999.

A. Bahrami and A. Teimourian, “Nonlocal scale effects on buckling, vibration and wave reflection in nanobeams via wave propagation approach,” Composite Structures, vol. 134, pp. 1061–1075, 2015.

M. R. Ilkhani, A. Bahrami, and S. H. Hosseini-Hashemi, “Free vibrations of thin rectangular nano-plates using wave propagation approach,” Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 40, no. 2, pp. 1287–1299, 2016.