Abstract

This paper gives an analytical approach for investigating free and forced transverse vibrations of a clamped-free rotating annular disk and obtaining the acoustic pressure distribution around the spinning disk. In the beginning, a modal analysis based on new analytical methods is carried out to find natural frequencies and mode shapes of the disk. Forced vibration of the disk is then investigated using Galerkin’s method. An analytical approach based on Laplace transformation is used to obtain time-dependent coefficients of the transverse response. A passive control strategy is examined for reducing the amplitudes of the transverse vibrations. The properties of the absorbers are examined in order to obtain the best performance. Rayleigh integral method and Durbin’s numerical Laplace transform inversion technique are adopted to compute the acoustic pressure around the rotating disk. Finally, a parametric study is performed and the effects of the design parameters as well as rotational conditions on the vibrational responses and the sound pressure of the spinning disk are examined.

1. Introduction

Many researchers have studied the dynamic behavior of the rotating disks and the acoustic pressure radiated by them because of their wide application in the engineering [1]. The importance of vibrations and acoustic radiation of rotating disks can be found in a different array of practical applications such as circular saws, hard disks, gas turbines, automobile parts, and aerospace structures [2]. In the past decades, some numerical, analytical, and experimental studies have been performed to examine the dynamic analysis of and the acoustic radiation from rotating disks. For example, Qiu et al. [3] examined an active control strategy to control the transverse vibration of a circular disk. The vibration and the noise reduction of an optical disk derived using a vibration absorber was studied by Heo et al. [4], and the required fundamental natural frequency of the absorber was obtained using a finite-element model. Ciğeroğlu and Özgüven [5] proposed a new model for the vibration analysis of turbine blades with dry friction dampers including both macroslip and microslip models representing dry friction dampers. Koo [6] analyzed the vibration and the critical speeds of polar orthotropic rotating annular disks employing the Rayleigh–Ritz method. In-plane free vibration of circular annular disks was studied by Bashmal et al. [7]. They presented a generalized formulation for the in-plane modal characteristics of circular annular disks under combinations of all possible classical boundary conditions. Hashemi et al. [8] performed vibration analysis of rotating thick plates based on Mindlin plate theory combined with second-order strain-displacement assumptions employing finite-element formulations. Damped vibrations of the double-sided fixed beams placed on a rotational disk were analyzed by Żółkiewski [9]. Younesian et al. [10] analytically studied vibration of a hollow circular plate subjected to a rotating peripheral force adopting Galerkin’s approach. The influence of shaft’s bending on the coupling vibration of a flexible blade-rotor system was analyzed by Li et al. [11]. Elastic stress analysis of a rotating annular disk made of functionally graded material (FGM) with variable thickness was studied by Jalali and Shahriari [12] employing the finite difference method. Bagheri and Jahangiri [13] studied the in-plane free vibration of the functionally graded rotating disks with variable thickness. An accurate solution for the in-plane vibration analysis of rotating circular panels with general edge restraints was proposed by Lyu et al. [14]. The natural frequencies and mode shapes of rotating turbo-machinery components from both rotating and stationary reference frames were experimentally analyzed by Presas et al. [15]. Yang et al. [16] analyzed thermos-elastic coupling vibration and stability of rotating circular plate in friction clutch. Nonlinear vibration analysis of turbine bladed disks with mid-span dampers was studied by Ferhatoglu et al. [17] utilizing Harmonic Balance Method, Alternating Frequency/Time approach, and Newton–Raphson method. Large amplitude vibrations of thin-walled rotating laminated composite cylindrical shell with arbitrary boundary conditions were examined by Li et al. [18] considering the large amplitude vibration of rotating shells with geometric nonlinearity. The coupling vibration characteristics of a rotating disk-beam system with the dovetail interfaces were analyzed by She et al. [19] based on a continuum model. On the other hand, the acoustic radiation from rotating disks or other engineering structures has attracted the attention of many researchers [20]. For instance, the acoustic radiation from out-of-plane modes of an annular disk using thin and thick plate theories was investigated by Lee and Singh [21]. They proposed a semianalytical procedure in which the disk surface velocity is numerically defined by a finite-element model. Reduction of flow-induced vibration and noise of an optical disk drive was performed by Cheng et al. [22]. Maeder et al. [23] studied numerical analysis of sound radiation from rotating disks using a simplified form of the Rayleigh integral known as the lumped parameter model. Recently, the average radiation efficiency of rotating annular plates in the rotating frame was studied analytically by Wang et al. [24] employing Galerkin's method and Rayleigh integral technique. Norouzi and Younesian [25] analyzed the transient and the steady-state sound radiation of nonlinear plates using analytical approaches.

Surveying the literature shows that studies have been addressed the vibrations of the rotating disk do not deal with proposing an analytical model to analyze the flexural vibrations, to mitigate the transverse response, and to reduce the sound radiated from the spinning disk. So, the main contribution of this paper is to represent a closed-form analytical solution in order to cover this gap in the literature. The results of this study can be significant in analyzing the disk brake squeal phenomena or in the study of the noise generated from the railway wheels.

2. Problem Statement

2.1. Analytical Model

An isotropic, homogeneous, thin annular disk of inner radius a, outer radius b, thickness h, Young’s modulus E, mass density ρ, and Poisson’s ratio ν subjected to a constant external load F is shown in Figure 1. The load is applied on the surface of the plate at the coordinate of (rp, θp). The disk has a constant angular speed Ω about the z-axis and is assumed to be rigidly clamped at its inner radius and free at outer radius.

Based on the classical plate theory, the transverse motion of the disk has no significant effect on the in-plane membrane forces of rotation. Moreover, the effects of the gravity and in-plane vibrations are also negligible [26]. So, the governing equation of transverse motion of the disk in cylindrical coordinates (r, θ, z) can be expressed as [2]where is the transverse deflection, c is the viscous damping, and is the rigidity of the disk. Moreover, and are the in-plane membrane stresses due to the centrifugal effect. also represents Dirac’s delta function. The boundary conditions for the clamped-free disk arein which Q and M are the shear force and the bending moment, respectively. One can employ Galerkin’s expansion method and assume (r, θ, t) as a periodic function of period 2π and write the response of equation (1) aswhere Amn(t) and Bmn(t) are time-dependent coefficients and Wm(r) is mode shapes of the disk obtained from free vibration analysis.

2.2. Free Vibration Analysis

One can let F = 0 in equation (1) in order to study the free vibration of the annular disk shown in Figure 1. So, the governing equation can be shown as

Using separation of variables technique (SV) and considering as a harmonic function, one may suppose the response of equation (4) aswhere and are the natural frequencies and the corresponding mode shapes, respectively, and m represents the circular mode number and n describes the diametric mode number. Additionally, using periodic property of the response, can also be expressed as

By substituting equation (6) into (5) and the result in equation (3), an ordinary differential equation describing W(r) is obtained as

The corresponding boundary conditions can be derived as follows:

For solving equation (7) with boundary conditions (8), one may multiply both sides of equation (7) by and perform integral from a to b for r and 0 to 2π for θ and then use the orthogonality property of the mode shapes of the disk [1, 20] and obtainHere, and are arbitrary parameters that show the upper bound for the mode numbers and they are both specified. Equation (9) forms an eigenvalue-eigenfunction problem. The solution of equation (9) can be written as follows [27]:where cj and RP (j= 1, 2,…,P) are unknown linear independent coefficients and the rest, respectively, and P is a certain positive integer which is chosen large enough such that the rest has a negligible error. Therefore, neglecting RP in equation (9), the application of conditions (8) leads to yield two linear equations of P + 1 unknown coefficients of c0, c1, c2, …, cP as follows:with

On the other hand, using conditions (8), one haswith

Equations (11) through (14) give four expressions about the unknowns of ci. But, according to equation (10), without considering RP,m, one needs obtaining P + 1 coefficients. So, another P − 3 independent equation should be produced. Inserting polynomial equation (10) into (9) and multiplying both sides of the result by ri (i = 4, 5, 6,…,P) and then integrating with respect to r between a and b lead towithin which the operator is

Therefore, equations (11), (13), and (15) form a system of P + 1 linear algebraic equations for P + 1 unknown coefficients ci (i = 0, 1,…,P), which can be rewritten in a compact form aswhere

To obtain a nontrivial solution of the system of equation (19), the determinant of the coefficient matrix has to vanish. Then, one gets a characteristic equation in eigenvalues λ as

Solving equation (20), one can obtain the natural frequencies of the disk. After finding the natural frequencies, the coefficient matrix c is obtained from equation (18) and, finally, equation (10) gives the mode shapes.

2.3. Forced Vibration Analysis

Recalling equations (1) and (3), one can replace equation (3) in (1) and multiply both sides of the result first by and then by and finally perform integration from 0 to 2π in θ and from a to b in r to give a system of coupled ordinary differential equations as follows:

Equation (21) can be solved employing Laplace transform technique. By taking Laplace transform from system of equation (21), one hasin which

Solving the algebraic equation (22) leads to obtain and as

The Laplace transform inversion scheme is now employed to obtain Amn(t) and Bmn(t) aswhere

And, is the unit imaginary number. Replacing equation (25) in conjunction with Wm(r) acquired in the previous section into equation (3), the transverse vibrational response of the disk is determined.

2.4. Vibration Absorption

In this section, the reduction in the vibrational deflection of the annular rotating disk is proposed by adding TMDs to the surface of the disk properly. Each TMD has the mass of mq, the stiffness of kq, and the viscous damping of cq (q= 1, 2, 3,…,Q), respectively, and is placed on the surface of the rotating disk at the coordinate of (rq, θq) as shown in Figure 2.

Based on equation (1), the equation of transverse motion of the disk-TMDs system can be expressed aswhere is the force acting from qth TMD on the surface of the disk. The equation of the motion of the qth TMD can also be shown as follows:in which uq=uq(t) is the displacement, is the natural frequency, and is the damping ratio of the dynamic absorber, respectively. Moreover, is the deflection of the disk at the point where qth TMD is connected. Putting Q = 1 and considering an absorber of mass ma, the natural frequency of , and the damping ratio of connected to the point of (ra, θa) on the surface of the disk, one can writewhere is the contact force created between the dynamic absorber and the disk. Equation (29) should be solved considering equation (28). Utilizing Galerkin’s approach, the solution can be expressed aswhere Emn(t) and Fmn(t) are unknown coefficients. Substituting equation (30) into (29) and multiplying both sides of the result first by and performing integration from 0 to 2π in θ and from a to b in r and then by and performing the same integrations constructs a system of coupled ordinary differential equations as follows:in which . Considering equations (29) and (28) and letting , one can obtain

Equation (32) is a set of ordinary differential equations which can be solved using numerical methods. The Runge–Kutta method is employed here in order to solve equation (32) and to determine the parameters of Eij(t), Fij(t), and ua(t). So, obtaining the parameters of Eij(t), Fij(t), and ua(t), one can obtain the transverse vibrational response of the system of disk-absorber using equation (30).

2.5. Acoustic Radiation

Figure 3 shows the disk of Figure 1 radiating the sound pressure. The sound pressure distribution P(r,θ, z, t) at an arbitrary point G can be obtained based on Rayleigh integral as follows [25, 28]:where denotes the acceleration of the disk, and are the mass density and the sound speed of the acoustic medium, respectively, and is the distance between the observation point G and an element on the surface of the disk at the location of (r, θ). Laplace transform technique can be employed again to obtain the acoustic pressure of equation (33). By taking Laplace transform from both sides of equation (33), one hasin which . Employing the Laplace transform inversion scheme, one can acquire the acoustic pressure. Durbin’s approach is used to obtain Laplace transform inversion of equation (34) as [29]where is an arbitrary real number greater than all the real parts of the singularities of and should be defined in the interval [0, 2τ]. For sufficient accuracy, the suggested value of “” is given with the appropriate sign by [30]where Nτ is the number of points in the time signal.

3. Results and Discussion

In this section, a parametric study is performed in order to investigate the effects of different parameters on the vibrational response and the acoustic pressure radiated from the spinning annular disk. In Table 1, the natural frequencies of the rotating disk obtained from the present study have been compared with the results from other references. The agreement between the results shows the accuracy of the proposed analytical approach. Table 2 also contains the geometrical properties of the disk and the physical properties of the acoustic medium.

Figures 4(a) through 4(f) show the mode shapes and the natural frequencies of the rotating annular disk for the first six modes. Time-dependent coefficients of Amn(t) and Bmn(t) are plotted in Figure 5 against the time. The figures have been plotted using two methods: the first graph in each figure by using equation (25) and the second one employing the numerical Runge–Kutta approach. There are excellent agreements between the results.

Figures 6(a)6(d) compare vibrational response of the disk observed at the point with coordinate of for and for different loading amplitudes. It can be found from these figures that increasing the amplitude of the external load causes the dynamical response to be more regular. For fully forced vibration of the disk (Figure 6(d)), one can see that the period of the response (0.4 sec) is the same as the rotational period.

The effect of the dynamic absorber attaching on the surface of the disk has been illustrated in Figures 7(a)7(d). The figures have been plotted with and without the absorber existence and for different frequencies of the TMD. It can be seen in Figure 7 that the efficiency of the TMD is variable in different natural frequencies. For finding the best natural frequency of the dynamic absorber, the root mean square (RMS) of the response has been shown against the natural frequency of the TMD in Figure 8. It is clear from Figure 8 that the efficiency of the absorber has the maximum value in a certain natural frequency (about 2100 Hz in this study). Moreover, the frequencies in which the deflection is not restrained well (about 1700 Hz and 2800 Hz), actually, are the natural frequencies of the system of disk-absorber. Figure 8 also shows that the performance of the TMD is impaired by increasing its natural frequency after a certain frequency (about 4 KHz in our study).

Figure 9 displays RMS of the transverse response versus the natural frequency of the TMD which is located in the different places on the surface of the disk. It is found from Figure 9 that the absorber location plays an important role in the performance of the vibration absorber. Figure 9 also shows that for the natural frequencies of the TMD after a certain frequency (about 3000 Hz in this study), the location of the TMD has no significant effect on the disk response. Moreover, at any location where the absorber placed on, the TMD has a maximum and a minimum effectiveness.

As it is shown in Figure 8, one can choose the best TMD from the natural frequency point of view. On the other hand, Figure 9 gives the best location on the surface of the disk in which the attached TMD has the best efficiency (407.5 mm from the center of the disk in this study). In fact, considering Figures 8 and 9, one can choose the best characters of TMD with the optimum performance. Figure 10 shows the time history of the transverse response of the disk before and after using the best dynamic absorber. From Figure 10, one can see a coincidence between two results at the beginning. This low quality of the TMD at the initial times is because the respite TMD needs to be adapted in the system.

The acoustic pressure radiated from the disk was plotted against the time in Figure 11. This figure compares the analytical model of present study with numerical models. A good agreement between the results is observed specially after the transient response.

Figure 12 focuses on the effect of annulus radiuses on the acoustic pressure. It is clear from the figure that increasing the ratio of outer to inner radius of the disk (i.e., b/a) causes dominant increasing in sound pressure. Figures 13(a) and 13(b) describe the effect of the rotational speed on the acoustic response. In Figure 13(a), this effect has been illustrated for different external forces while Figure 13(b) shows these evolutions when the thickness of the disk changes. Form Figures 13(a) and 13(b), one can find a certain frequency (about 800 Hz in our study) in which the sound radiated from the spinning disk generates its minimum pressure. This frequency is independent of any changes in the force amplitude and the disk thickness.

Figures 14(a) through 14(c) demonstrate time snapshots of the acoustic pressure distribution at the plane parallel to the surface of the rotating disk. The figures were generated when the disk has been rotated a half cycle, one, and one-and-a-half cycles, respectively. It can be found from the figure that the acoustic pressure has the minimum and the maximum values in certain points around the disk, but the field is completely symmetric. Comparing Figures 14(a) through 14(c), one can see the changes in the sound pressure intensity and the distribution pattern. It is also clear from Figures 14(a) and 14(c) that the sound pressure will be amplified in the higher cycles of rotation, but the distribution pattern will be periodic.

The time snapshots of the acoustic pressure field at the plane parallel to the surface of the system of dynamic absorber-rotating disk showing the location of the TMD have been illustrated in Figures 15(a) through 15(c). Comparing these figures with Figures 14(a) through 14(c), one can obtain the acoustic pressure decreasing after using TMD on the surface of the rotating disk. Furthermore, the acoustic field at the plane parallel to the surface of the disk is asymmetric when the dynamic absorber is attached. One can observe the reduction in the acoustic pressure radiated from the vibrating disk after using TMD.

4. Conclusions

The analytical modeling of free and forced vibration and the acoustic radiation mitigation for an annular rotating disk was studied in this paper. Galerkin’s expansion method was employed to obtain the vibrational response of the disk. The natural frequencies and the mode shapes of the disk were found by analytical methods. Time-dependent coefficients were obtained analytically, and the results were compared with the numerical simulation results. A tuned-mass-damped (TMD) mechanism was added to the main system in order to reduce the amplitude of the vibrations. The influence of the different parameters on the efficiency of the TMD was examined. Acoustic pressure radiated from the rotating annulus was computed employing Rayleigh integral method, Laplace transform technique, and Durbin’s Laplace transform inversion scheme. The most important results can be mentioned as follows:(i)The acoustic pressure radiated from the vibrating rotating annular disk is increased by increasing the ratio of outer-inner radius.(ii)The rotating disk generates different sound levels in different spinning speeds. Moreover, there is a certain speed in which the sound pressure radiated from the disk has a minimum value. This merit rotational speed is independent of forcing conditions.(iii)Rotating disks with higher thicknesses are noisier than the thin disks.(iv)The acoustic pressure field generated by the disk is symmetric, but adding TMD to the system makes the sound field asymmetric.

Data Availability

The data that support the findings of this study are available from the corresponding author, upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest in preparing this article.