An analytical investigation on the dynamic stability of a rotor filled with liquid

Li, Bobo; Wang, Guangding; Yuan, Huiqun

doi:10.21595/jve.2018.19886

Journal of Vibroengineering

Browse Journal

Submit article

Published: 30 September 2018

Check for updates

An analytical investigation on the dynamic stability of a rotor filled with liquid

Bobo Li¹

Guangding Wang²

Huiqun Yuan³

^{1, 2}School of Mechanical Engineering and Automation, Northeastern University, Shenyang, China

³School of Science, Northeastern University, Shenyang, China

Corresponding Author:

Huiqun Yuan

Cite the article Download PDF

Downloads 1606

WoS Core Citations 2

CrossRef Citations 3

Abstract

This paper deals with the dynamic stability of a rigid rotor arbitrarily filled with liquid. On the basis of the established coupled-field equations of the rotor system, the general whirling eigenequation, which is a quartic complex coefficients equation, is derived. In order to obtain the solutions of the general whirling eigenequation, a mathematical method is proposed. To illustrate the precision of calculating results, a comparison is carried out between the present analysis and the numerical results. The results show that two calculation results are in good agreement. Then the stability of the rotor system is analyzed. It is shown that the dynamic instability occurs at a particular bound of the spinning speed. Moreover, the effects of system parameters, such as fluid-fill ratio and mass ratio, on the unstable regions are discussed.

An analytical investigation on the dynamic stability of a rotor filled with liquid

Highlights

A mathematical method for obtaining the solutions of the general whirling eigen-equation was proposed.
The whirling frequency equation of the rigid liquid-filled rotor was presented to predict the system stability.
The dynamic instability occurs at a particular bound of the spinning speed.
The effects of some main parameters on system stability were discussed.

1. Introduction

Rotor filled with liquid is applied widely in engineering. The inside of rotor is hollow, and the chamber contains liquid. When the rotating cylinder suffers the perturbation in working, the whirling motion occurs. An additional pressure, which acts on rotating cylinder internal surface, results from the relative rotation between liquid and rotor. This forms a complex fluid-solid coupling problems. Meanwhile, due to the interaction between rotating fluid and rotor, the phenomenon of rotor system instability happens at a certain RPM range. Therefore, it has very important significance to study dynamic stability of rotor filled with liquid.

The dynamic stability of rotor filled with liquid was first studied by Kollmann, Jlnnouchi, and Crandall [1-4]. The critical speed of Jeffcott single-disc rotor filled with liquid was calculated by Kollmann. The instability phenomenon was observed on the rotor partially filled with liquid in the experiment. Crandall analyzed the dynamic stability of a flexible rotor filled with liquid. After that, the dynamic characteristics of parallel whirling motion were analyzed for rotor filled with liquid by Wolf [5]. In the early study, the liquid of rotor chamber was assumed to be ideal and rotor body was rigid. For this reason, the dynamic stability was investigated for rotor filled with viscous liquid by Hendricks, Daich and Bar, Saito and Holm-Christensen [6-9]. Furthermore, rigid body angular displacement was introduced by Hendricks and Lichtenberg, the research on rotor filled with liquid was extended to three dimensions [10, 11]. Hiroshi Ota et al. made experiments on the vibration of rotor filled with liquid, the effect of viscosity and specific gravity on the unstable vibration were studied simultaneously [12]. Yuichi SATO proposed a dynamic absorber using a hollow rotor partially filled with liquid, a kind of technique was described to suppress structure vibration, and the effectiveness of proposed absorber was showed experimentally [13].Kuipers analyzed the stability of a flexibly mounted rotating cylinder partially filled with liquid. The rotor was assumed to rotate around its axis of symmetry with a constant angular velocity and the liquid was ideal in the chamber. The unstable region was derived and the effect of combined action of dissipative and gyroscopic terms was researched by him [14]. Zhang et al. further discussed the dynamic stability of rigid and flexible rotor filled with liquid [15].

In recent years, there are many related researches on the stability of rotor filled with liquid. The experimental investigation was made on the dynamic characteristics of an over-hung rigid centrifuge rotor partially filled with liquid by Zhu [16].The influence of fluid-fill ratio on rotor whirl frequency, the range of unstable region, and the rotor imbalance were researched. Furtherly, the stability of rotor filled with viscous incompressible fluid was analyzed by him. The dynamic stability was investigated on flexible rotating cylinder partially filled with liquid by Tao and Zhang, they proposed that the three dimensional flow can be simplified to two dimensional flow problem in the chamber [17]. Yukio Ishida et al. applied a method, which utilizes discontinuous spring characteristics, to suppress unstable vibrations in an asymmetrical rotor system, an asymmetrical shaft system, and a rotor system partially filled with liquid [18].N. V. Derendyaev et al. investigated the stability of a cylinder filled with different kinds of liquid, including viscous incompressible liquid, a stratified viscous incompressible liquid and a stratified non-uniform viscous incompressible liquid [19-22]. In 1988, Derendyaev proved that Andronov-Hopf bifurcations may occur in a rotor system filled with liquid [23]. In his work, the methodological mistakes made by Wolf Jr. et al. were analyzed [24].

Many researchers have investigated the vibration and fault diagnosis for a rotor system, but little information is available on the analytical solutions of the general whirling eigenequation [25-33]. In this paper, considering the complex form of the equation, a mathematical analysis method is proposed to obtain the solutions of the general whirling eigenequation. The effective of the method is verified by comparing with numerical method. Then, according to whirling eigenequation, the dynamic stability of a liquid-filled rotor system is analyzed, and the effects of fluid-fill ratio and mass ratio on unstable region and critical speed of rotor system are discussed.

2. General whirling eigenequation

The structure of rotor filled with liquid is shown in Fig. 1. Rotor is a hollow cylinder, the length of cylinder is $L$ , inner diameter is $a$ , mass is $m_{R}$ , centroid is $c$ . The liquid is filled in cylinder chamber, both ends of rotor are supported on bearings symmetrically through two flexible shafts, the distance between rotating shaft and liquid level is $b$ . Static coordinate $o$ - $x y z$ is established based on the axis of rotor under static state. When rotor rotates at an angular velocity $Ω$ , the rotation is assumed to be steady and unperturbed, rotor and liquid rotate synchronously, rotor axis $o$ coincides with liquid centroid $c$ , as shown in Fig. 1(b).

Fig. 1Structure of rotor filled with liquid

a)

b)

When rotor suffers perturbation in working, whirling motion is formed in chamber. We assume that the radial deviation of axis is $ε$ , rotor whirling frequency is $ω$ . Here, liquid centroid and rotor axis do not coincide with each other, rotor of centroid $c$ rotates with angular velocity $ω$ , which is shown in Fig. 2.

It is assumed that the liquid is ideal; in addition, $η (t)$ and $ξ (t)$ are relative perturbation displacements of the centroid $c$ in the direction of $η$ and $ξ$ . Then, the whirling kinetic equation of rotor centroid $c$ can be written as follows:

1

[\begin{array}{l} m_{R} \\ m_{R} \end{array}] \{\begin{array}{l} \ddot{η} \\ \ddot{ξ} \end{array}\} + [\begin{matrix} c & - 2 m_{R} Ω \\ 2 m_{R} Ω & c \end{matrix}] \{\begin{array}{l} \dot{η} \\ \dot{ξ} \end{array}\} + [\begin{matrix} k - m_{R} Ω^{2} & - c Ω \\ c Ω & k - m_{R} Ω^{2} \end{matrix}] \{\begin{array}{l} η \\ ξ \end{array}\} = \{\begin{array}{l} F_{η} \\ F_{ξ} \end{array}\},

where $c$ is external damping coefficient, $k$ is elastic support stiffness, $F_{η}$ and $F_{ξ}$ are force caused by fluid acts on the rotor, $Ω$ is self-rotation angular velocity.

Fig. 2Parallel whirling of a rotor filled with liquid

Based on fluid mechanics, Euler's equation of fluid in rotor chamber can be expressed as:

2

\{\begin{array}{l} \frac{\partial u}{\partial t} - 2 Ω v = - \frac{\partial G}{\partial r}, \\ \frac{\partial v}{\partial t} + 2 Ω u = - \frac{1}{r} \frac{\partial G}{\partial θ}, \end{array}

where the parameter $G$ is given by:

3

G = \frac{p_{1}}{ρ} + (a_{η} η + a_{ξ} ξ) r,

where $p_{1}$ is perturbation pressure, $a_{η}$ and $a_{ξ}$ are components of centroid $c$ absolute acceleration along the direction of $η$ and $ξ$ .

The continuity equation of fluid in rotor chamber can be written as:

4

\frac{u}{r} + \frac{\partial u}{\partial r} + \frac{1}{r} \frac{\partial v}{\partial θ} = 0 .

For solving the fluid-solid coupling equation of rotor whirling motion above, first and foremost, boundary conditions should be determined. We assume that fluid free surface is $r = b + ξ (θ, t)$ , here $ξ (θ, t)$ is perturbation infinitesimal displacement. Then boundary condition of fluid free surface can be represented as:

5

\frac{\partial p_{1}}{\partial t} |_{r = b} = - ρ Ω^{2} b \frac{\partial ξ}{\partial t} = - ρ Ω^{2} b u |_{r = b},

and fluid displacement boundary condition takes the form:

6

u |_{r = a} = 0 .

Based on Eqs. (1), (2) and (4), external damping coefficient is neglected, we can derive the general whirling motion eigenequation for rigid rotor arbitrarily filled with liquid by using boundary conditions Eqs. (5) and (6). The eigenequation is as follows:

7

(z^{2} + 1) {(s - z i)}^{2} + δ (2 s^{2} - 4 s z i - z^{2}) [1 + (μ + 1) z^{2}] = 0,

where $δ$ , $μ$ and $s$ are defined as:

δ = \frac{e^{2} - 1}{2}, μ = \frac{m_{c}}{m_{R}}, s = \frac{Ω}{ω_{0}},

where $m_{c}$ is the liquid mass, $ω_{0}$ is the nature frequency of a solid rotor. Then, we can expand Eq. (7) as follows:

8

(γ + μ) z^{4} + i 2 α s z^{3} + (γ - α s^{2}) z^{2} + i 2 (γ + 1) s z - (γ + 1) s^{2} = 0,

where $γ$ and $α$ can be shown as:

γ = \frac{a^{2} + b^{2}}{a^{2} - b^{2}}, α = γ + 2 μ + 1 .

3. Analytical solution for general whirling eigenequation

3.1. Solution procedure

Dynamic stability of rotor arbitrary filled with liquid is a quite complex problem. Furthermore, the general whirling eigenequation is a quartic complex coefficients equation. In order to investigate rotor dynamic stability, the eigenequation needs to be analyzed further. To this end, derivational process for the analytical solution of the eigenequation is given as below.

Eq. (8) takes the form:

9

- \frac{(γ + 1) s^{2}}{γ + μ} + i [\frac{2 (γ + 1) s}{γ + μ}] z + (\frac{γ - α s}{γ + μ}) z^{2} + i (\frac{2 α s}{γ + μ}) z^{3} + z^{4} = 0 .

Let $a = - \frac{(γ + 1) s^{2}}{γ + μ}$ , $b = i [\frac{2 (γ + 1) s}{γ + μ}]$ , $c = (\frac{γ - α s}{γ + μ})$ and $d = i (\frac{2 α s}{γ + μ})$ , Eq. (9) can be written in the simplified form:

10

a + b z + c z^{2} + d z^{3} + z^{4} = 0 .

We assume that $z_{1}$ , $z_{2}$ , $z_{3}$ and $z_{4}$ are four unknown roots of Eq. (10), let $- (z_{1} + z_{2}) = A$ , $z_{1} z_{2} = B$ , $- (z_{3} + z_{4}) = C$ and $z_{3} z_{4} = D$ , Eq. (10) can be written as:

11

(z - z_{1}) (z - z_{2}) (z - z_{3}) (z - z_{4}) = (z^{2} + A z + B) (z^{2} + C z + D) = 0 .

Then, Eq. (11) can be expanded as follows:

12

z^{4} + z^{3} (A + C) + z^{2} (B + D + A C) + z (A D + B C) + B D = 0 .

Let us suppose that eigenequation has two conjugate complex roots, or two real roots. There is an independent complex root inevitably among the others. Therefore, if $A$ and $B$ are real numbers, $C$ and $D$ become complex numbers certainly.

Let $a = a_{1} + i a_{2},$ $b = b_{1} + i b_{2},$ $c = c_{1} + i c_{2},$ $d = d_{1} + i d_{2},$ $C = C_{1} + i C_{2}$ and $D = D_{1} + i D_{2}$ , because Eq. (12) is equal to Eq. (10), the coefficients of two equations should be equal. According to complex number algebra, the relationships of coefficients real part for two equations are given by:

13

A + C_{1} = d_{1},

14

B + D_{1} + A C_{1} = c_{1},

15

A D_{1} + B C_{1} = b_{1},

16

B D_{1} = a_{1} .

In addition, we introduce an unknown quantity $P$ , which gives that:

17

P = B + D_{1} .

Substituting Eq. (17) into Eq. (14), we obtain:

18

P = c_{1} - A C_{1} .

By Eq. (14) and Eq. (16), we can derive:

19

2 a_{1} P = 2 a_{1} c_{1} - 2 A B C_{1} D_{1} .

Then we can square both sides of Eq. (15):

20

b_{1}^{2} = A^{2} D_{1}^{2} + 2 A B C_{1} D_{1} + B^{2} C_{1}^{2} .

For eliminating $A$ , $B$ , $C$ and $D$ , the Eq. (19) and Eq. (20) are added:

21

2 a_{1} P = 2 a_{1} c_{1} - b_{1}^{2} + A^{2} D_{1}^{2} + B^{2} C_{1}^{2} .

Square both sides of Eq. (13), and multiply by $a_{1}$ , we have:

22

a_{1} d_{1}^{2} = A^{2} D_{1} B + 2 A B C_{1} D_{1} + D_{1} B C_{1}^{2} .

Adding to Eq. (19), then we get:

23

2 a_{1} P = 2 a_{1} c_{1} - a_{1} d_{1}^{2} + A^{2} D_{1} B + D_{1} B C_{1}^{2} .

Add Eq. (21) to get:

24

\begin{array}{l} 4 a_{1} P = 4 a_{1} c_{1} - b_{1}^{2} - a_{1} d_{1}^{2} + A^{2} D_{1} (B + D_{1}) + (B + D_{1}) B C_{1}^{2} \\ = 4 a_{1} c_{1} - b_{1}^{2} - a_{1} d_{1}^{2} + (A^{2} D_{1} + B C_{1}^{2}) P . \end{array}

Here, removing $(A^{2} D_{1} + B C_{1}^{2})$ , by Eq. (13) and Eq. (15), we obtain:

25

b_{1} d_{1} = (A^{2} D_{1} + B C_{1}^{2}) + A C_{1} (B + D_{1}) .

Both sides are multiplied by $P$ , then the Eq. (25) becomes:

26

P b_{1} d_{1} = (A^{2} D_{1} + B C_{1}^{2}) P + A C_{1} P^{2} .

Subtract Eq. (24) to get:

27

P (b_{1} d_{1} - 4 a_{1}) = - (4 a_{1} c_{1} - b_{1}^{2} - a_{1} d_{1}^{2}) + A C_{1} P^{2} .

Both sides of Eq. (18) are multiplied by $P^{2}$ so that:

28

P^{3} = c_{1} P^{2} - A C_{1} P^{2} .

Add Eq. (27) to get:

29

P^{3} - c_{1} P^{2} + P (b_{1} d_{1} - 4 a_{1}) - a_{1} d_{1}^{2} - b_{1}^{2} + 4 a_{1} c_{1} = 0 .

Eq. (29) can be solved and the roots take the form:

P_{1} = c_{1}, P_{2} = 2 \sqrt{a_{1}}, P_{3} = - 2 \sqrt{a_{1}},

where:

a_{1} = - \frac{(γ + 1) s^{2}}{γ + μ}, b_{1} = 0, c_{1} = (\frac{γ - α s}{γ + μ}), d_{1} = 0 .

Then, $B$ and $D_{1}$ can be determined from Eq. (16) and Eq. (17), respectively:

B, D_{1} = \frac{P \pm \sqrt{P^{2} - 4 a_{1}}}{2} .

$A$ and $C_{1}$ are solved by Eq. (13), Eq. (14) and Eq. (17), respectively:

A, C_{1} = \frac{d_{1} \pm \sqrt{d_{1}^{2} - 4 c_{1} + 4 P}}{2} .

The values of $A$ , $B$ , $C_{1}$ , $D_{1}$ should satisfy Eq. (15). According to the values of $P$ , we can have a discussion as follows.

Case 1:

P_{1} = c_{1} .

In this case, two sets of values are obtained, respectively:

A = C_{1} = 0, B = \frac{(c_{1} + \sqrt{c_{1}^{2} - 4 a_{1}})}{2}, D_{1} = \frac{(c_{1} - \sqrt{c_{1}^{2} - 4 a_{1}})}{2},

and:

A = C_{1} = 0, B = \frac{(c_{1} - \sqrt{c_{1}^{2} - 4 a_{1}})}{2}, D_{1} = \frac{(c_{1} + \sqrt{c_{1}^{2} - 4 a_{1}})}{2} .

Case 2:

P_{2} = 2 \sqrt{a_{1}} .

In this condition, we get:

A = \sqrt{- c_{1} + 2 \sqrt{a_{1}}}, C_{1} = - \sqrt{- c_{1} + 2 \sqrt{a_{1}}}, B = D_{1} = \sqrt{a_{1}},

and:

A = - \sqrt{- c_{1} + 2 \sqrt{a_{1}}}, C_{1} = \sqrt{- c_{1} + 2 \sqrt{a_{1}}}, B = D_{1} = \sqrt{a_{1}} .

Case 3:

P_{3} = - 2 \sqrt{a_{1}} .

If $P$ takes this value, accordingly, the values of $A$ , $B$ , $C_{1}$ , $D_{1}$ take the form:

A, C_{1} = \pm \sqrt{- c_{1} - 2 \sqrt{a_{1}}}, B, D_{1} = - \sqrt{a_{1}} .

Due to $- c_{1} - 2 \sqrt{a_{1}} <$ 0, such situation does not exist. Only case 1 and case 2 are logical.

Now let us consider the imaginary parts of $C$ and $D$ . From Eq. (14), Eq. (15) and Eq. (16), we obtain:

30

\{\begin{array}{l} D_{2} = \frac{a_{2}}{B} = 0, \\ b_{2} - B C_{2} = 0, \\ c_{2} - A C_{2} = 0 . \end{array}

Because $A$ , $B$ , $C$ , $D$ can’t all be real, $C_{2} \neq$ 0, it follows that:

A = 0, C_{2} = \frac{b_{2}}{B} .

In conclusion, $A$ , $B$ , $C$ , $D$ can be solved, we have:

A = 0, B = c_{1} + \frac{\sqrt{c_{1}^{2} - 4 a_{1}}}{2}, C = \frac{i b_{2}}{B}, D = c_{1} - \frac{\sqrt{c_{1}^{2} - 4 a_{1}}}{2},

A = 0, B = c_{1} - \frac{\sqrt{c_{1}^{2} - 4 a_{1}}}{2}, C = \frac{i b_{2}}{B}, D = c_{1} + \frac{\sqrt{c_{1}^{2} - 4 a_{1}}}{2},

where:

a_{1} = - \frac{(γ + 1) s^{2}}{γ + μ}, c_{1} = \frac{γ - α s}{γ + μ}, b_{2} = \frac{2 (γ + 1) s}{γ + μ} .

The roots of Eq. (8) take the form:

31

\{\begin{matrix} z_{1}, z_{2} = \frac{- A \pm \sqrt{A^{2} - 4 B}}{2}, \\ z_{3}, z_{4} = \frac{- C \pm \sqrt{C^{2} - 4 D}}{2} . \end{matrix}

Finally, substituting the values for $A$ , $B$ , $C$ , $D$ into Eq. (31), four roots of Eq. (8) can be obtained.

3.2. Precision analysis

In Sec. 3.1, the method is utilized to solve general whirling eigenequation of rotor arbitrarily filled with liquid. To verify the precision of analytical results, a comparison is carried out between analytical solution and numeric results.

Let $μ =$ 0.206 and $γ =$ 2.6, here, $μ$ is mass ratio and $γ$ is fill parameter, eigenequation can be written in the form:

32

2.806 z^{4} + i 8.024 s z^{3} + (2.6 - 4.012 s^{2}) z^{2} + i 7.2 s z - 3.6 s^{2} = 0 .

If $s =$ 1, we can get the results of Eq. (32) by using analytical and numerical methods respectively. The comparison results are shown in Table 1.

Table 1Comparison results of two methods

Eigenvalue	Numerical results	Analytical results	Error (%)
$z_{1}$	0.95002 $i$	0.95326 $i$	0.33
$z_{2}$	–0.90932 $i$	–0.95326 $i$	4.8
$z_{3}$	–0.66415 $i$	–0.64935 $i$	–2.2
$z_{4}$	–2.23614 $i$	–2.17435 $i$	–2.3

From Table 1, we can observe that the analytical results accord with numerical results well, maximum error is only 4.8 %. Consequently, the derived analytical results can achieve a high accuracy, which shows the effective of the calculation method for solving general whirling motion eigenequation of rotor arbitrarily filled with liquid.

4. Dynamic stability analysis for rotor filled with liquid

4.1. System stability analysis

For researching dynamic stability of a rotor filled with liquid, we assume that the mass ratio $μ =$ 0.206 and the parameter $γ =$ 2.6, let $s =$ 0, 1, 2, 3, 4 and 5, respectively. Then the roots of eigenequation can be obtained. Fig. 3 shows the polar diagram of the root. From Fig. 3, it can be observed that the polar angles of roots are 0, $π$ and $- π$ , which explains that eigenvalues only contain imaginary part. Therefore, the solution for eigenequation can be represented as:

33

z = - i ξ .

Substituting Eq. (33) into Eq. (8), we have:

34

(γ + μ) ξ^{4} - 2 α s ξ^{3} + (α s^{2} - γ) ξ^{2} + 2 (γ + 1) s ξ - (γ + 1) s^{2} = 0 .

Based on Eq. (34), the plot of dimensionless whirling frequency $ξ$ versus dimensionless spinning frequency $S$ for $μ =$ 0.206 and $γ =$ 2.6 can be obtained, as shown in Fig. 4. We can see that there are two real roots and a pair of conjugate complex roots in the interval [ $a_{1}$ , $a_{2}$ ], which illustrates that the system is unstable. Moreover, point A denotes the rotor first-order critical speed in the Fig. 4 when $ξ = S$ , the critical spinning speed can be obtained as:

35

s_{c r}^{2} = \frac{1}{1 + μ} .

Point A is on the left side of $a_{1}$ , namely, instability rotating speed is greater than first order critical speed. Therefore, one can say that rotor passes critical speed first and then system instability is occurred. Let $ξ =$ 1, we obtain $s = 1 \pm \sqrt{2} / 2$ , which is corresponding to point $b_{1}$ and $b_{2}$ in Fig. 4. Thus, the estimation range of unstable region is:

36

\frac{1}{\sqrt{1 + μ}} < s < 1 + \frac{\sqrt{2}}{2} .

From Eq. (34), it can be found that the system stability is related to the fluid-fill ratio and the mass ratio. In the following section, we will discuss the influences of system parameters on the stability of a liquid-filled rotor system.

Fig. 3Polar diagram of the roots for the system eigenequation

a) $S =$ 0

b) $S =$ 1

c) $S =$ 2

d) $S =$ 3

e) $S =$ 4

f) $S =$ 5

Fig. 4Dimensionless whirling frequency ξ versus dimensionless spinning frequency S

4.2. The influence of parameters on the system stability

4.2.1. Fluid-fill ratio

In this paper, the function of external damping coefficient is neglected in parameters discussion. Let $x = 1 - b / a$ ; here $x$ is fluid-fill ratio, then we have:

37

γ = \frac{(a^{2} + b^{2})}{(a^{2} - b^{2})} = \frac{[1 + (1 - x)^{2}]}{[1 - (1 - x)^{2}]} .

We assume that the mass ratio is a constant, given that $μ =$ 0.206. Thus, the effects of the fluid-fill ratio on system stability are shown in Figs. 5-7.

Fig. 5The effects of the fluid-fill ratio x on system stability at μ=0.206

Fig. 6Variation of unstable region when fluid-fill ratio x is small

Fig. 7Variation of unstable region when fluid-fill ratio x is large

In Fig. 5, we can see that when fluid-fill ratio x is small, unstable region increases with $x$ . On the contrary, if the fluid-fill ratio is large, the unstable region will decrease gradually. Fig. 6 shows the variation of unstable region when fluid-fill ratio is small. It can be observed that when fluid-fill ratio is changed in the range of (0, 0.3], the lower and upper bounds of unstable region increase with $x$ . Due to the increasing speed of upper bound is greater than that of lower bound, the unstable region increases with fluid-fill ratio and moves towards high speed zone gradually.

Fig. 8Variation of the width of unstable region versus fluid-fill ratio x

Fig. 7 shows the variation of unstable region when fluid-fill ratio $x$ is large. When fluid-fill ratio $x$ is changed in the range of (0.3, 1], the lower bound of unstable region increases as $x$ increases. However, the upper bound increases firstly, and then decreases. As a result, unstable region becomes smaller with the increasing of fluid-fill ratio $x$ . When $x =$ 1, unstable region is almost disappeared, which indicates that the rotor system is stable. The width of unstable region changes with fluid-filled ratio $x$ is shown in Fig. 8. The results show that as fluid-fill ratio $x$ increases, the unstable region first expand, then narrow. If the chamber of rotor is fully filled with fluid, the instable phenomenon is disappeared.

4.2.2. Mass ratio

When we discuss the influence of mass ratio $μ$ on the system stability, fluid-fill ratio $x$ should be taken as a constant. Let fluid-fill parameter $γ =$ 2.6, the effects of mass ratio $μ$ on unstable region can be obtained, as shown in Fig. 9. It shows that with the growth of mass ratio $μ$ , the unstable region increases gradually.

Fig. 9The effects of the mass ratio μ on system stability for γ= 2.6

Fig. 10Variation of the width of unstable region versus mass ratio μ

The variation of the width of unstable region versus with mass ratio $μ$ is shown in Fig. 10. One can observe that when mass ratio $μ$ increases, the upper bound of unstable region has little changed, but lower bound decreases obviously. Therefore, unstable region presents the tendency of increase. Moreover, the unstable region moves towards low speed zone gradually.

4.2.3. The effects of parameters changed together on unstable region

In Section 4.2.1 and Section 4.2.2, we have analyzed the effects of individual parameter on unstable region. In this section, we will take into account the combined action of fluid-fill ratio $x$ and mass ratio $μ$ . The variations of unstable region are shown in Fig. 11 and Fig. 12.

Fig. 11 shows the effects of fluid-fill ratio $x$ on unstable region under different mass ratio $μ$ . It can be observed that as fluid-fill ratio $x$ increases, the unstable region first expands, then narrows. When the rotor chamber is fully filled with liquid, unstable region disappears, and the rotor system is stable. As mass ratio increases, unstable region moves to low speed zone gradually. From Fig. 11, it can be seen that the mass ratio $μ$ has little influence on the width of unstable region. However, the fluid-fill ratio x gives great effects, which illustrates that fluid-fill ratio is the primary factor to affect the system stability. Moreover, as fluid-fill ratio $x$ increases, the first critical speed of the system does not vary.

Fig. 12 gives the effects of mass ratio $μ$ on unstable region for different fluid-fill ratio $x$ . It shows that the unstable region increases with the growth of mass ratio $μ$ . The fluid-fill ratio $x$ has a great influence on the width of unstable region, as shown in Fig. 12. As fluid-fill ratio $x$ increases, the width of unstable region increases obviously. From Fig. 12, it can be clearly seen that the first order critical speed decreases with the increase in mass ratio $μ$ , which indicates that the first order critical speed is only determined by mass ratio $μ$ .

Fig. 11The effects of fluid-fill ratio x on unstable region under different mass ratio μ: a) μ= 0.206, b) μ= 0.4, c) μ= 0.6, d) μ= 0.8

a)

b)

c)

d)

Fig. 12The effects of mass ratio μ on unstable region under different fluid-fill ratio x: a) x= 0.1, b) x= 0.2, c) x= 0.3, d) x= 0.4

a)

b)

c)

d)

5. Conclusions

An analytical solution has been introduced to solve the general whirling eigenequation of a rotor filled with liquid. The analytical results of the eigenequation was obtained and the validity of the proposed method was illustrated in the paper. A comparison was carried out between analytical calculation and numerical method. The results show that the maximum error of the two results is only 4.8 %, which illustrates the reasonability of the present calculation method. On the basis of whirling eigenequation, the system stability was analyzed and the effects of system parameters on unstable region were discussed. The performed investigations show that the dynamic instability occurs at a particular bound of the rotating speed for a rotor filled with liquid. The first order critical speed, which is only related to the mass ratio $μ$ , is less than the instability spinning speed for the rotor system. Based on the present work, it could be found that the unstable region moves to the low speed zone as mass ratio $μ$ increases. The fluid-fill ratio $x$ is the prime factor affecting the width of the unstable region. When the rotor chamber is fully filled with liquid, unstable region disappears, and the system is stable.

References

Kollmann F. G. Experimentelle und theoretische Untersuchungen uber die kritischen Drehzahlen flüssigkeitsgefullter Hohlkorper. Forschung Auf Dem Gebiet Des Ingenieurwesens A, Vol. 28, Issue 4, 1962, p. 115-123.

Publisher
Jlnnouchi Y., Araki Y., Inoue J. The dynamic instability of a non-rotating pipe in a high speed rotor filled with liquid. Transactions of the Japan Society of Mechanical Engineers, Vol. 54, Issue 502, 1988, p. 1210-1216.

Publisher
Crandall S. H. Rotating and Reciprocating Machines. Handbook of Engineering Mechanics. New York, 1962.

Search CrossRef
Crandall S. H., Mroszczyk J. Whirling of a flexible cylinder filled with liquid. Proceedings of International Conference on Rotor Dynamics, Tokyo, 1986.

Search CrossRef
Wolf J. A. Whirl Dynamics of a Rotor Partially Filled with Liquid. Transaction of the ASME Journal of Applied Mechanics, Vol. 35, Issue 4, 1968, p. 676-682.

Publisher
Hendricks S. L., Morton J. B. Stability of a rotor partially filled with a viscous incompressible fluid. Journal of Applied Mechanics, Vol. 46, Issue 4, 1979, p. 913-918.

Publisher
Daich I. M., Bar I. L. Oscillations of a rotating solid body with a cavity partly filled with viscous fluid. Prikl Mekh, Vol. 9, Issue 5, 1973, p. 64-69.

Publisher
Saito S., Someya T. Self-excited vibration of a rotating hollow shaft partially filled with liquid. Journal of Mechanical Design, Vol. 102, Issue 1, 1980, p. 185-192.

Publisher
HolmChristensen O., Trager K. A Note of instability caused by liquid motions. Journal of Applied Mechanics, Vol. 58, Issue 3, 1991, p. 804-811.

Publisher
Hendricks S. L. Stability of a clamped-free rotor partially filled with liquid. Journal of Applied Mechanics, Vol. 53, Issue 1, 1986, p. 166-172.

Publisher
Lichtenberg G. Vibrations of an elastically mounted spinning rotor partially filled with liquid. Journal of Applied Mechanics, Vol. 104, Issue 2, 1982, p. 389-396.

Publisher
Ota H., Ishida Y., Sato A., Yamada T. Experiments on vibrations of a hollow rotor partially filled with liquid. Bulletin of JSME, Vol. 29, Issue 256, 1986, p. 3520-3529.

Publisher
Sato Y. Dynamic absorber using a hollow rotor partially filled with liquid. Transactions of the Japan Society of Mechanical Engineers, Vol. 54, Issue 3, 1988, p. 446-452.

Publisher
Kuipers M. On the stability of a flexibly mounted rotating cylinder partially filled with liquid. Applied Scientific Research, Vol. 13, Issue 1, 1964, p. 121-137.

Publisher
Zhang W., Tang J., Tao M. Dynamic stability of a rotor filled or partially filled with liquid. Journal of Applied Mechanics, Vol. 63, Issue 1, 1996, p. 101-105.

Publisher
ZhuChangsheng Experimental investigation into the instability of an over-hung rigid centrifuge rotor partially filled with fluid. Journal of Vibration and Acoustics, Vol. 124, Issue 4, 2002, p. 483-491.

Publisher
Tao M., Zhang W. Dynamic stability of flexible rotating cylinder partially filled with liquid. Journal of Fudan University, Vol. 38, Issue 2, 1999, p. 141-145.

Search CrossRef
Ishida Y., Liu J. Elimination of unstable ranges of rotors utilizing discontinuous spring characteristics: an asymmetrical shaft system, an asymmetrical rotor system, and a rotor system with liquid. Journal of Vibration and Acoustics, Vol. 132, Issue 1, 2010, p. 325-326.

Publisher
Derendyeav N. V. Sandalov V. M. To stability of steady-state rotation of cylinder partly filled with viscous incompressible liquid. Journal of Applied Mathematics and Mechanics, Vol. 46, Issue 4, 1982, p. 578-586.

Publisher
Derendyaev N. V. Stability of steady rotation of a cylinder filled with a stratified viscous incompressible liquid. Fluid Mechanics and Heat Transfer, Vol. 272, Issue 5, 1983, p. 1073-1076.

Search CrossRef
Derendyaev N. V., Senyatkin V. A. Stability conditions for the steady-state rotation of a cylinder filled with a stratified nonuniform viscous incompressible liquid. Journal of Applied Mechanics, Vol. 25, Issue 1, 1984, p. 30-39.

Publisher
Derendyaev N. V., Soldatov I. N. Study of stability and self-excited vibrations of a rotor with liquid on the basis of a discrete model. Proceedings of the VI Scientific Conference on Nonlinear Vibrations of Mechanical Systems, Nizhny Novgorod, 2002.

Search CrossRef
Derendyaev N. V. Andronov-Hopf bifurcation in the dynamics of a liquid-filled rotor system. Soviet Physics Doklady, Vol. 33, Issue 8, 1988, p. 592-594.

Search CrossRef
Derendyaev N. V. Nonconservative Problems of Dynamics of Rotor Systems Containing Fluid. Dissertation of Doctor of Science, Nizhny Novgorod, 1999.

Search CrossRef
Wu D., Rong C., Bing H. A novel two-stage hybrid swarm intelligence optimization algorithm and application. Soft Computing, Vol. 16, Issue 10, 2012, p. 1707-1722.

Publisher
Wu D., Xinhua Y., Li Z. An improved self-adaptive differential evolution algorithm and its application. Chemo Metrics and Intelligent Laboratory Systems, Vol. 128, Issue 15, 2013, p. 66-76.

Publisher
Wu D., Huimin Z., Jingjing L. An improved CACO algorithm based on adaptive method and multi-variant strategies. Soft Computing, Vol. 19, Issue 3, 2015, p. 701-713.

Publisher
Wu D., Huimin Z., Li Z. A novel collaborative optimization algorithm in solving complex optimization problems. Soft Computing, Vol. 21, Issue 15, 2017, p. 4387-4398.

Publisher
Wu D., Huimin Z., Xinhua Y. Study on an improved adaptive PSO algorithm for solving multi-objective gate assignment. Applied Soft Computing, Vol. 59, 2017, p. 288-302.

Publisher
Zhao H.M., Li D.Y., Deng W. Research on vibration suppression method of alternating current motor based on fractional order control strategy. Proceedings of the Institution of Mechanical Engineers Part E-Journal of Process Mechanical Engineering, Vol. 231, Issue 3, 2017, p. 786-799.

Publisher
Wu D., Rui Y., Huimin Z. A novel intelligent diagnosis method using optimal LS-SVM with improved PSO algorithm. Soft Computing, 2017, https://doi.org/10.1007/s00500-017-2940-9.

Publisher
Wu D, Shengjie Z, Huimin Z. A novel fault diagnosis method based on integrating empirical wavelet transform and fuzzy entropy for motor bearing. IEEE Access, 2018.

Search CrossRef
Huimin Z., Meng S., Wu D. A new feature extraction method based on EEMD and multi-scale fuzzy entropy for motor bearing. Entropy, Vol. 19, Issue 1, 2016, p. 14.

Publisher

Cited by

The Inertial Disturbances of Fluid Movement in the Chamber of a Liquid Autobalancer

Vasyl Lozynskyi | Thaer Shihab | Ilona Drach | Liubomyr Ropyak

(2024)

The Chaotic Property of BTA Deep-Hole Machining System under the Effect of Inner Cutting Fluid

Quanbin Zhang | Wu Zhao | Yamin Li | Shuangxi Jing | Francisco Beltran-Carbajal

(2023)

Dinamic evaluation of a gas turbine vibration to determine oil inside the axial compressor

Luiz Homero Lopes Martins | Roberto Firmento De Noronha

(2020)

About this article

Received

13 April 2018

Accepted

01 July 2018

Published

30 September 2018

SUBJECTS

Mechanical vibrations and applications

DOI

https://doi.org/10.21595/jve.2018.19886

Keywords

rotor filled with liquid

dynamic stability

parameters

unstable region

critical speed

Acknowledgements

The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The work is supported by the National Science Foundation of China (Grant No. 51775093) and the National Science Foundation of China (Grant No. 51335003).

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.

Previous article in issue Previous Next article in issue Next

Research article

2023 08 03

Evaluation of the flow state and static performance of smooth annular liquid seals

Zhongjie Yang, Jin Feng, Jingyu Wan, Xiaobin Yu, Zhangjin He, Wanfu Zhang

Research article

2020 09 30

Improving the rotordynamic stability of short labyrinth seals using positive preswirl

Wanfu Zhang, Qianlei Gu, Hao Cao, Yingfei Wang, Lu Yin

Research article

2020 03 31

Theoretical and experimental analysis of an unbalanced and cracked cardan shaft in the vicinity of the critical speed

Bernard Xavier Tchomeni, Alfayo Alugongo

Research article

2019 09 30

Numerical analysis of a vertical rotor rub system partially immersed in inviscid incompressible fluid

Desejo Filipeson Sozinando, Bernard Xavier Tchomeni, Alfayo Anyika Alugongo

B. Li, G. Wang, and H. Yuan, “An analytical investigation on the dynamic stability of a rotor filled with liquid,” Journal of Vibroengineering, Vol. 20, No. 6, pp. 2253–2267, Sep. 2018, https://doi.org/10.21595/jve.2018.19886

Copy Extrica

Copied to clipboard!

TY  - JOUR
DO  - 10.21595/jve.2018.19886
UR  - https://doi.org/10.21595/jve.2018.19886
TI  - An analytical investigation on the dynamic stability of a rotor filled with liquid
T2  - Journal of Vibroengineering
AU  - Yuan, Huiqun
AU  - Li, Bobo
AU  - Wang, Guangding
PY  - 2018
DA  - 2018/09/30
PB  - JVE International Ltd.
SP  - 2253-2267
IS  - 6
VL  - 20
SN  - 1392-8716
SN  - 2538-8460
ER  - 

Copy Ris

Copied to clipboard!

@article{Yuan_2018,
	doi = {10.21595/jve.2018.19886},
	url = {https://doi.org/10.21595/jve.2018.19886},
	year = 2018,
	month = {sep},
	publisher = {{JVE} International Ltd.},
	volume = {20},
	number = {6},
	pages = {2253--2267},
	author = {Huiqun Yuan and Bobo Li and Guangding Wang},
	title = {An analytical investigation on the dynamic stability of a rotor filled with liquid},
	journal = {Journal of Vibroengineering}
}

Copy Bibtex

Copied to clipboard!

[1]H. Yuan, B. Li, and G. Wang, “An analytical investigation on the dynamic stability of a rotor filled with liquid,” Journal of Vibroengineering, vol. 20, no. 6, pp. 2253–2267, Sep. 2018, doi: 10.21595/jve.2018.19886.

Copy IEEE

Copied to clipboard!

Yuan, Huiqun, Bobo Li, and Guangding Wang. “An Analytical Investigation on the Dynamic Stability of a Rotor Filled with Liquid.” Journal of Vibroengineering 20, no. 6 (September 30, 2018): 2253–67. https://doi.org/10.21595/jve.2018.19886.

Copy Chicago

Copied to clipboard!

An analytical investigation on the dynamic stability of a rotor filled with liquid

Abstract

Highlights

1. Introduction

2. General whirling eigenequation

3. Analytical solution for general whirling eigenequation

3.1. Solution procedure

3.2. Precision analysis

4. Dynamic stability analysis for rotor filled with liquid

4.1. System stability analysis

4.2. The influence of parameters on the system stability

4.2.1. Fluid-fill ratio

4.2.2. Mass ratio

4.2.3. The effects of parameters changed together on unstable region

5. Conclusions

References

Cited by

About this article

Related Articles