Abstract

In order to improve the performance and efficiency of the rotating cylindrical shell (RCS), one of the effective ways reduces the mass of the RCS. The scientific and effective method is to design the variable thickness of RCS (VTRCS) along the axial direction in response to this demand. Then, the vibration traveling wave characteristics of VTRCS are investigated by Chebyshev–Ritz method. The displacement field of the cylindrical shell is expanded in the form of the product of the Chebyshev polynomial and the boundary function. The kinetic and potential energies of the VTRCS are calculated based on the Sanders shell theory considering the effects of Coriolis forces and centrifugal forces. Also, the frequency equation of the VTRCS is obtained according to the Ritz method. The accuracy of the modeling method is verified by comparing the present results with literature that had done, and the convergence of the calculated results is studied. Finally, the free-vibration traveling wave characteristics of the VTRCS in different thickness variations are compared, and the parameters such as the rotational speed, thickness variation parameters, and aspect ratio on the free-vibration traveling wave characteristics of the VTRCS are discussed. The influence of thickness variation on the vibration characteristics is analyzed, which has significance for the lightweight design for VTRCS.

1. Introduction

The rotating cylindrical shell (RCS) is widely used in the industrial field, such as high-speed rotating centrifugal separator and rotor system of aeroengine. The RCS is increasingly demanding complex working condition in the engineering equipment; meanwhile, the wall thickness of the RCS shows the trend of more and more thin to meet the increasing of performance and efficiency, which lead to the structural vibration becoming more and more prominent. In addition, the RCS will have special traveling wave vibration phenomenon due to the influence of centrifugal force, Coriolis force, and the initial hoop tension compared with the static cylindrical shell; thus, it is of great significance to study the vibration characteristics of the RCS.

There are various solutions for the free vibration analysis of the RCS; for example, some scholars used Galerkin’s method to analyze the free vibration characteristics of the RCS [15], and the differential quadrature method was also used to study the traveling wave properties of the free vibration for the RCS [611]. In addition, Ritz method was relatively simple in calculation process and could maintain high accuracy, so it was widely used in the research of vibration characteristics of RCS. For instance, Yan et al. [12] solved an accurate solution for the vibration-damping characteristics of a constrained stand-off-layer damping cylindrical shell and obtained the loss factor and the frequency based on the Hamilton principle. Song et al. [13] presented an improved Rayleigh–Ritz method to investigate the vibration and damping behavior of thin short cylindrical shell. Nguyen et al. [14] established sandwich structures and derived equations for the free vibration of rotating shells by the first-order, shear-deformation shell theory. Abbaspour and Hosseini [15] employed the Ritz method with Chebyshev polynomials to perform the free vibration analysis of the graphene platelets reinforced piezoelectric cylindrical microshells with various boundary conditions. Qin et al. [1618] took Chebyshev polynomials as the admissible functions, employed the Rayleigh–Ritz method to derive the motion equations for the rotating shell-plate combination, investigated the vibration of an RCS coupled with an annular plate based on the Sanders shell theory and Mindlin plate theory, and calculated the strain energy of the shell and plate. Then, they studied the free vibration of rotating functionally graded carbon nanotube reinforced composite cylindrical shells with arbitrary boundary conditions. Some scholars utilized Ritz method and Gram-Schmidt polynomials to form approximate functions and analyzed the free vibration of RCS [19]. Furthermore, some scholars studied the nonlinear vibration of cylindrical shells of composite materials [2024].

In addition, the abovementioned studies did not consider the influence of the thickness variation of the RCS. However, sometimes, the RCS needs to be designed with variable thickness in practice; that is, the thickness varies along the axial direction, and this can further reduce the weight. Therefore, some scholars analyzed the traveling wave vibration characteristics of rotating cylindrical shells with variable thickness (VTRCS). Quoc et al. [25] used Galerkin method to study the vibration characteristics of VTRCS under thermal environment, but only one thickness variation form was considered in this study.

Based on the above discussion, three different thickness variations are considered in this study, and Chebyshev–Ritz method is used to investigate the free vibration traveling wave behavior of VTRCS. Furthermore, the influences of parameters such as rotational speed, thickness variation parameters, and aspect ratio of cylindrical shells on the traveling wave behavior of free vibration for the VTRCS are discussed, which has significance for the lightweight design of rotating cylindrical shell structures.

High speed rotary drum has important application in centrifuge, which is abstracted to the rotating cylindrical shell (RCS) to reduce its mass and improve its performance and efficiency. The scientific and effective method is to design the variable thickness of RCS along the axial direction. In addition, the research can be applied to the aircraft engine disk or other mechanical structures and industrial fields, such as submarine, ship, drilling oil and gas field development, and factory workshop. Therefore, the research in this work has important application value in the fields of national defense and military.

2. Theoretical Modeling

Firstly, a theoretical model is established to study the traveling wave vibration characteristics of VTRCS. The cylindrical shell structure rotates around its central axis with a rotational speed Ω, its length is L, and its average radius is R. The (x, θ, z) is an orthogonal coordinate system established on the surface of the cylindrical shell, and u, , and are the dispositions of any point on the cylindrical shell along the three directions of x, θ, and z, respectively, as shown in Figure 1. It is assumed that the thickness h(x) of the rotating cylindrical shell varies linearly along its axis, as shown in Figure 2, which can be divided into three variation forms, denoted as V1, V2, and V3, respectively.

The coordinates of the upper and lower surfaces of the cylindrical shell in the Z direction of the coordinate axis will change under the three different thickness changes. The specific expressions are as follows.

The thickness change form of V1 is expressed as follows:where h1(x) and h2(x), respectively, represent the coordinates of the upper and lower surfaces of the cylindrical shell in the Z direction of the coordinate axis; h0 represents the initial thickness of the cylindrical shell, that is, the thickness at x = 0; kh represents the thickness variation parameter.

The thickness change form of V2 is expressed as follows:

The thickness change form of V3 is expressed as follows:

2.1. Solution of Kinetic and Potential Energy

The Chebyshev–Ritz method is used to study the VTRCS in order to solve the natural frequency of the model. Firstly, the kinetic energy equation and potential energy equation of the cylindrical shell need to be given.

The velocity vector of any point on the VTRCS can be expressed as follows:where k is the unit vectors along z direction and r represents the displacement vector of the VTRCS at any point on the coordinate system (x, θ, z), which can be described as follows:where i, j, and k are, respectively, unit vectors along x, θ, and z directions.

The kinetic energy of the VTRCS can be calculated as follows:

Then, the kinetic energy equation of VTRCS is obtained by substituting equation (4) into equation (6), which is denoted as follows:

According to Sanders shell theory, the strain at any point on the rotating cylindrical shell is written as follows:where , , and are the strain components of the middle surface and , , and are the curvature components of the middle surface. They can be expressed by the following formula [26]:

The stress-strain relationship of the VTRCS can be expressed as follows:where .

The strain energy of the VTRCS caused by the deformation can be described as follows:

The strain energy of the VTRCS caused by centrifugal force [27] is expressed as follows:

Therefore, the total potential energy equation of the VTRCS with variable cross section can be expressed as follows:

The total potential energy solution of the VTRCS lays a foundation for studying its natural frequency and mode shape.

2.2. Solution of Natural Frequency and Mode Shape

The natural frequency and mode shape of the VTRCS are solved in this section, and the solution process is given on the basis of the kinetic energy equation and potential energy equation obtained by the Chebyshev–Ritz method in the previous section.

Firstly, the displacement field of the VTRCS can be expressed as follows:where is the natural frequency and n is the toroidal wave number of the traveling wave modal. U(x), V(x), and W(x) are modal functions, which are approximately expanded by the product of Chebyshev polynomials and their corresponding boundary functions, and the specific formula is expressed as follows:where ai, bj, and ck are unknown coefficients and , , and are the first kind Chebyshev expressions, which can be expressed in trigonometric terms as follows:

Chebyshev polynomials can achieve fast convergence speed and maintain high accuracy at a lower computational cost, but they are defined on the interval [−1, 1] and have orthogonality only on this interval, so coordinate transformation is needed, that is, . , , and represent boundary functions along the direction, which need to satisfy the corresponding geometric boundary conditions of the VTRCS, and the specific expressions are shown in Table 1.

Secondly, the displacement field equation (14) of the VTRCS is substituted into the kinetic energy equation (7), which can be further expressed as follows:

Similarly, equation (14) is substituted into equation (13), and the potential energy equation of the VTRCS can be expressed as follows:where A11, A12, and A66 represent the tensile stiffness and D11, D12, and D66 represent the bending stiffness. The specific expressions are denoted as follows:

Then, according to equations (17) and (18), the energy expression of the VTRCS can be obtained as follows:

The most likely approximate value is determined by minimizing the energy expression with respect to the unknown coefficient according to Rayleigh’s principle, so the derivative of the energy expression with respect to the unknown coefficient is performed as follows:

Then, equation (21) can be further transformed into an eigenvalue matrix, which is written as follows:where K is the stiffness matrix, M1 and M2 are the mass matrices, represents the natural frequency, and is the eigenvector composed of unknown coefficients, namely, the mode shape of the VTRCS, which is expressed as follows:

The natural frequencies and corresponding mode shapes can be obtained by solving equation (22).

3. Numerical Results and Discussion

In this section, the calculation results are compared with those in the existing literature, and their convergence is studied to verify the accuracy of the established model. Then, the free vibration traveling wave behavior of the VTRCS is analyzed through parameter research. Unless otherwise mentioned, the length of the VTRCS L = 1 m, the initial thickness h0 = 0.02 m, the thickness variation parameter kh = 0.5, the density ρ = 1072 kg/m3, the elastic modulus E = 172 GPa, the shear modulus G = 4.2 GPa, and the Poisson’s ratio  = 0.31, the average radius R = 0.2, the speed  = 25 r/s, the dimensionless frequency parameter  = , and the dimensionless speed  = .

3.1. Comparison and Convergence Research

Two cases are studied in this section to verify the accuracy and convergence of the proposed method.

3.1.1. Case One

The dimensionless frequency parameters (DFPs) of the backward wave and the forward wave for the RCS with uniform thickness under fixed (C-C) boundary conditions at both ends are given, respectively, and compared with the results in literature [27]. Also, the convergence research is conducted, and the calculation results under different number of truncation terms are listed, as shown in Tables 2 and 3. In this example, the aspect ratio of the RCS L/R = 10, the thick-aspect ratio h/R = 0.05, the Poisson’s ratio  = 0.3, the dimensionless frequency parameter  = , and the dimensionless speed  = .

It can be seen from Tables 2 and 3 that the results are basically consistent with the results in literature [27]. The study shows that the calculation results converge to a certain value with the increase of the number of truncated terms nmax, which have converged to a sufficiently accurate value when the number of truncated terms is 11, so nmax is 11 in the following calculation.

3.1.2. Case Two

The second case gives the natural frequency (Hz) of the RCS with uniform thickness under simply supported (S-S) boundary conditions at both ends as shown in Table 4, which is compared with the results in literature [28]. The cylindrical shell length L = 0.256 m, average radius R = 0.16 m, thickness h = 0.0025 m, elastic modulus E = 110 GPa, Poisson’s ratio  = 0.31, density ρ = 4480 kg/m3, axial half wave number m = 1, and rotational speed Ω = 20000 r/min.

It is observed in Table 4 that the results obtained by the proposed model are basically consistent with those in literature [28], and the maximum error is less than 1%.

In short, the above two examples verify the correctness and convergence of the proposed method.

3.2. Parameter Study

The effects of different thickness variations, rotational speed, thickness variation parameters, aspect ratio, and thickness-to-diameter ratio on the traveling wave behavior of free vibration are discussed in order to obtain the vibration characteristics of the VTRCS under different thickness variations, as shown in Figures 37.

The variation of the DFPs and with the circumferential wave number n of the VTRCS is shown in Figure 3. In the figure, V0: kh = 0, that is, the thickness of the RCS is uniform and remains unchanged in the x-axis direction.

Figure 3 shows that the DFPs and for the RCS under thickness change rise with the increase of the wave number n in the boundary conditions of the S-S or C-C. The DFPs of the VTRCS are the maximum under V0 thickness changes form, and it is the minimum under V1 thickness changes form. The DFPs of the VTRCS under V2 thickness variation form are higher than those under the other two thickness variations when V0 thickness variation form is not considered, which are close to those under V0 thickness variation form. It is also obvious from Figure 3 that the DFPs of the VTRCS under the four thickness changes firstly disperse from the same initial value and then converge to the same value with the increase of circumferential wave number n.

The variation of the DFPs and with the thickness parameter kh for the RCS under the three thickness forms when the rotational speed is, respectively, 0 r/s, 25 r/s, and 50 r/s, is shown in Figure 4.

It is seen from Figure 4 that the DFPs for the VTRCS decrease gradually with the increase of kh under different rotational speed and thickness, and it is minimum in the condition of the V2 thickness change. In addition, it can also be found that the rotational speed has almost no influence on the variation trend of the kh.

The variation of the DFP as the changing of Ω is shown in Figure 5 to further study the influence of rotational speed on the traveling wave for the VTRCS. The BW represents the backward wave, and FW represents the forward wave.

Figure 5 shows that the DFPs of backward wave rise with the increase of the rotational speed for the VTRCS, and the DFPs of forward wave reduce gradually with the increase of rotational speed. The DFP value of backward wave and forward wave under different rotational speed is the biggest under V0 thickness change form. The DFP value increases, and their difference becomes smaller under V2 thickness variation form when the boundary conditions change from S-S to C-C. In addition, it can also be found that the difference between the DFPs under different thicknesses becomes significantly larger when the circumferential wave number changes from 1 to 2.

Finally, the influence of geometric parameters on the DFPs of VTRCS is studied in Figures 6 and 7. Figure 6 gives the variation of the DFP for the VTRCS with the changing of the aspect ratio L/R, and Figure 7 shows the changing of the aspect ratio R/h0.

Figure 6 manifests that the DFPs of VTRCS decrease gradually with the increase of length-to-diameter ratio L/R. Particularly, the DFPs decrease rapidly when the length-to-diameter ratio is less than 3, and it decreases slowly when the length-to-diameter ratio is more than 3. It can be seen from Figure 7 that the backward wave and the forward wave both decrease first and then increase with the increase of the thickness-to-diameter ratio for the DFP for the VTRCS. It can also be found from Figure 7 that the difference between the DFPs under the four thickness variation forms becomes smaller as the boundary conditions change from S-S to C-C. Specifically, the influence of thickness variation forms on the DFPs is little when the thickness-to-diameter ratio is higher than 25.

4. Conclusion

(1)The vibration characteristics of the VTRCS are investigated based on the Sanders shell theory via the Chebyshev–Ritz method. The correctness and convergence of the modeling method are verified through comparison and convergence research, and it is proved that the proposed model can effectively predict the traveling wave behavior of the free vibration for the VTRCS.(2)The traveling wave DFP of the VTRCS rises with the increase of the wave number n under three thickness variation forms. The traveling wave DFP is the largest under V0, and it is higher than the other two kinds of thickness under V2 when V0 is not considered, which is close to the DFP in the V0. In addition, the DFP of the VTRCS first discretizes from the same initial value and then converges to the same value with the increase of circumferential wave number n under the four thickness variations.(3)The traveling wave DFP of the VTRCS gradually decreases with the increase of the thickness change parameter kh at different speeds and different thicknesses. Among them, the traveling wave DFP change is the smallest under V2 and is significantly smaller than that of the other two thickness changes. The rotational speed has little effect on the variation trend of traveling wave frequency with thickness.(4)The traveling wave DFP of VTRCS decreases gradually with the increase of aspect ratio L/R and decreases first and then increases with the increase of aspect ratio h0/R. The difference between the DFPs becomes smaller under the four thickness variations as the boundary conditions change from S-S to C-C. The variation of thickness has little effect on the traveling wave frequency when the thickness-to-diameter ratio is higher than a certain value.

Data Availability

The data used to support the findings of this study are currently under embargo while the research findings are commercialized. Requests for data, [6/12 months] after publication of this article, will be considered by the corresponding author.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

The authors gratefully acknowledge the financial supports for this research from National Natural Science Foundation of Changsha City in Hunan Province (Grant no. kq2208084), Natural Science Research of Education Department in Hunan Province (Grant no. 21C1235), National Natural Science Foundation of Hebei Province (Grant no. E2020202217), and Open Research Fund for Key Laboratory of National Defense Science and Technology for Particle Transport and Enrichment Technology (Grant no. STPTS202111).