Analytical analysis on damping characteristics of rotor system
Min Wu^{1} , Shengbo Yang^{2} , Wen Bao^{3} , Jinfu Yang^{4}
^{1, 2, 3}School of Energy Science and Engineering, Harbin Institute of Technology, Harbin, China
^{4}Institute of Engineering Thermophysics, Chinese Academy of Sciences, Beijing, China
^{1}Corresponding author
Journal of Vibroengineering, Vol. 19, Issue 3, 2017, p. 1697-1708.
https://doi.org/10.21595/jve.2017.18133
Received 25 December 2016; received in revised form 10 April 2017; accepted 12 April 2017; published 15 May 2017
JVE Conferences
A general vibration model of a flexible rotor system is established to investigate the influences of the damping characteristics on vibration behaviors. Based on the multi-scale method, the analytical solutions of steady-state and transient-state are derived under the positive and negative conditions of nonlinear damping. The physical significance of the coefficients and their influences on rotor behaviors are analyzed through theoretical analysis and numerical calculation. The experimental results elaborate the damping effect and verify the rationality of the model.
Keywords: flexible rotor system, general vibration model, analytical solution, nonlinear damping, vibration mechanism.
1. Introduction
Rotating machinery is widely used in generators, gas turbines, aero engines and other important power systems. It is the key component and core technology in the national defense industry and energy field such as aerospace, ship engineering and distributed energy. In the development towards high efficiency and high performance, due to its complicated structure, precise design and rigorous assembly, a high-speed flexible rotor system often occurs various faults and shows strong nonlinear behavior characteristics in practical operations. Thus, the stability mechanism of nonlinear rotor system becomes vital to meet the higher requirements of theoretical researches and engineering applications.
For a general rotor system, external damping is considered to be able to reduce the amplitude in the entire operating range and increase the system performance against interference. On the contrary, internal damping is considered to increase the amplitude in operating range, cause oil/gas film whirl, whip and other nonlinear phenomena, and even lead to system failure and halt. The internal damping is caused by the tangential force when stress center line and strain center line of rotor cross-section are no longer coincide. If external damping is insufficient to overcome the self-excited vibration caused by internal damping, the rotor system may become instable.
Newkirk [1] first discovered this phenomenon. Kimball [2] considered the instability was caused by internal friction. Timoshenko [3] did the quantitative analysis and derived the constitutive relation between stress and strain. Gunter [4] and Ehrich [5] also studied the effects of internal friction/damping on the stability of the rotor system respectively. Due to the limitations of linear theory in the analysis of nonlinear phenomena, the accuracy of linearization cannot meet the engineering requirements, so the theory of oil film instability develops from the eigenvalue criterion stability theory based on linear assumption, to the nonlinear stability theory based on nonlinear simulation such as energy method, spectrum analysis method, pumping method [6] and circumferential average velocity ratio method [7]. Tasker [8] used moving-block technique and the sparse time domain method to estimate the equivalent nonlinear damping from transient response data. Chandra [9] used the envelope of free vibration signal extracted by continuous wavelet transform approach to identify nonlinear external and internal damping, and the approximate analytical solution was obtained by Krylov-Bogobliubov method. The nonlinear dynamic stability research has received more and more attention [10].
In this paper, by introducing nonlinear damping force and nonlinear stiffness force, a general form of nonlinear vibration model is proposed for a high-speed flexible rotor system. The steady-state solution and transient-state solution are deduced respectively under positive and negative nonlinear damping conditions via multi-scale method. The effects of linear and nonlinear damping on rotor vibration characteristics are discussed through analytical analysis and numerical simulation, and are elaborated through experimental results.
2. General model establishment
Consider a general ideal single-span single-disc rotor of disc mass $m$, eccentricity $r$, rotational frequency $\mathrm{\Omega}$, as shown in Fig. 1. Shaft mass is ignored, and counterclockwise is positive. Though the rotor structure is simple, it can fully reflect the linear and nonlinear dynamic behaviors of the rotor system.
Fig. 1. Schematic of force analysis of a rotor system
For general liquid or gas bearing, in this paper consider only the bearing support force, ignoring the tangential force. Based on Newton’s second law, equation of motion in the $x$-axis direction can be described as:
where:
Generally, the rotor force can be divided into stiffness and damping force by the physical properties, and also can be divided into linear and nonlinear force by mathematical properties, namely:
where ${f}_{k}$ is elastic force, ${f}_{c}$ is damping force, ${f}_{l}$ is linear force, ${f}_{n}$ is nonlinear force, $\xi $ is linear scale factor, $\epsilon $ is nonlinear scale factor, $\mathrm{\Theta}=\epsilon /\xi $ is mix ratio. The nonlinear terms have been discussed in detail in previous work [11, 12]. Considering the nonlinear Hooke’s law [13-15] and nonlinear damping force [16, 17], the ${f}_{k}$ and ${f}_{c}$ are given as:
In the formulae above, the first order approximation are the traditional linear force expressions${f}_{k}=-kx$ and ${f}_{c}=-c\dot{x}$. If $n=0$ for ${f}_{c}$, then ${f}_{c}$ becomes dry friction which is irrelevant to velocity, so it is not considered. In this paper, the nonlinear forces are taken as the first two orders approximation:
where $k$ is linear stiffness coefficient, $\eta $ is nonlinear stiffness coefficient, $c$ is linear damping coefficient, $\gamma $ is nonlinear damping coefficient, ${\xi}_{k}=k$ is stiffness linear scale factor, ${\xi}_{c}=c$ is damping linear scale factor, ${\mathrm{\Theta}}_{k}=\eta /k$ is stiffness mix ratio, ${\mathrm{\Theta}}_{c}=\gamma /c$ is damping mix ratio, the positive sign before ${\mathrm{\Theta}}_{c}$ is for positive nonlinear damping and negative sign for negative nonlinear damping. By taken Eq. (2), Eq. (3) and Eq. (5) into Eq. (1), the general form of nonlinear vibration model for a rotor system can be obtained:
For briefness and clarity, some substitutions are given as:
where ${\omega}_{0}$ is natural frequency, ${\beta}_{0}$ is nonlinear stiffness amplification coefficient, $\mu $ is nondimensional linear damping coefficient, ${\alpha}_{0}$ is nonlinear damping amplification coefficient. Thus ${\xi}_{k}={\omega}_{0}^{2}m$, ${\mathrm{\Theta}}_{k}=\epsilon {\beta}_{0}^{2}/{\omega}_{0}^{2}$, ${\xi}_{c}=$ 2 μm, ${\mathrm{\Theta}}_{c}=\epsilon {\alpha}_{0}^{2}/2\mu $, and there are:
The equation is the nonlinear vibration model based on large deformation and large perturbation assumption for high-speed flexible rotor system. It is worth noting that, the positive sign before velocity square term indicates positive nonlinear damping, and the negative sign indicates negative nonlinear damping.
3. Deduction and analytical solutions
The quadratic approximate solution of Eq. (8) is deduced via multi-scale method:
Take Eq. (9) into Eq. (8), according to the power orders of $\epsilon $, and $D$ is a partial differential operator, there are:
The solution of Eq. (10) is:
where $C$ is an unknown complex function, the overline represents the conjugate function. Take Eq. (12) into Eq. (11):
where $cc$ represents the conjugate functions of each term at the right side of the equal sign. The coefficient of first term ${e}^{{r}_{1}{T}_{0}}$ must be zero to avoid duration term. At any point without acceleration, ${T}_{0}=\text{0}$, so ${e}^{-2\mu {T}_{0}}=\text{1}$. There is:
Assuming $C=a{e}^{j\phi}/2$, $C$ can be obtained:
where ${a}_{0}$, ${\phi}_{0}$ are the integration constants which can be determined by initial conditions. Take Eq. (15) into Eq. (13) so ${x}_{1}$ can be solved, then the solution $x$ is:
where:
${A}_{1}=-\frac{3\epsilon {\beta}_{0}^{2}}{4{\mathrm{\Omega}}^{2}r}{A}_{0}^{4},{A}_{2}=\frac{\epsilon {\alpha}_{0}^{2}}{2\sqrt{{\left({\lambda}^{2}-4\right)}^{2}+4{\delta}^{2}}}{A}_{0}^{2},{A}_{3}=-\frac{\epsilon {\beta}_{0}^{2}}{4{\mathrm{\Omega}}^{2}\sqrt{{\left({\lambda}^{2}-9\right)}^{2}+9{\delta}^{2}}}{A}_{0}^{3},$
${A}_{4}=-\frac{\epsilon {\alpha}_{0}^{2}}{2{\lambda}^{2}}{A}_{0}^{2},{A}_{01}=-\frac{3\epsilon {\beta}_{0}^{2}{a}_{0}^{2}}{2{\mathrm{\Omega}}^{2}r}{A}_{0}^{2},{A}_{40}=-\frac{1}{2}\epsilon {\alpha}_{0}^{2}{a}_{0}^{2},$
${A}_{11}=-\frac{\epsilon {\alpha}_{0}^{2}{a}_{0}\lambda}{\sqrt{4{\lambda}^{2}-{\delta}^{2}}+1}{A}_{0},{A}_{{11}^{\text{'}}}=-\frac{\epsilon {\alpha}_{0}^{2}{a}_{0}\lambda}{\sqrt{4{\lambda}^{2}-{\delta}^{2}}-1}{A}_{0},$
${A}_{21}=\frac{3\epsilon {\beta}_{0}^{2}{a}_{0}^{2}}{4{\mathrm{\Omega}}^{2}\sqrt{{\left(3{\lambda}^{2}+1-{\delta}^{2}+2\sqrt{4{\lambda}^{2}-{\delta}^{2}}\right)}^{2}+{\left(\delta \sqrt{4{\lambda}^{2}-{\delta}^{2}}+\delta \right)}^{2}}}{A}_{0},$
${A}_{{21}^{\text{'}}}=\frac{3\epsilon {\beta}_{0}^{2}{a}_{0}^{2}}{4{\mathrm{\Omega}}^{2}\sqrt{{\left(3{\lambda}^{2}+1-{\delta}^{2}-2\sqrt{4{\lambda}^{2}-{\delta}^{2}}\right)}^{2}+{\left(\delta \sqrt{4{\lambda}^{2}-{\delta}^{2}}-\delta \right)}^{2}}}{A}_{0},$
${A}_{12}=\frac{3\epsilon {\beta}_{0}^{2}{a}_{0}}{8{\mathrm{\Omega}}^{2}\left(\sqrt{4{\lambda}^{2}-{\delta}^{2}}+2\right)}{A}_{0}^{2},{A}_{{12}^{\text{'}}}=\frac{3\epsilon {\beta}_{0}^{2}{a}_{0}}{8{\mathrm{\Omega}}^{2}\left(-\sqrt{4{\lambda}^{2}-{\delta}^{2}}+2\right)}{A}_{0}^{2},$
${A}_{20}=\frac{\epsilon {\alpha}_{0}^{2}{a}_{0}^{2}\sqrt{{\left(-2{\lambda}^{2}+{\delta}^{2}\right)}^{2}+{\delta}^{2}\left(4{\lambda}^{2}-{\delta}^{2}\right)}}{4\sqrt{{\left(3{\lambda}^{2}-{\delta}^{2}\right)}^{2}+{\delta}^{2}\left(4{\lambda}^{2}-{\delta}^{2}\right)}},$
${A}_{30}=\frac{\epsilon {\beta}_{0}^{2}{a}_{0}^{3}}{4{\mathrm{\Omega}}^{2}\sqrt{{\left(8{\lambda}^{2}-3{\delta}^{2}\right)}^{2}+9{\delta}^{2}\left(4{\lambda}^{2}-{\delta}^{2}\right)}},{\theta}_{0}=-\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{t}\mathrm{g}\frac{\delta}{{\lambda}^{2}-1},{\theta}_{1}=2{\theta}_{0},$
${\theta}_{2}=-\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{t}\mathrm{g}\frac{2\delta}{{\lambda}^{2}-4}+2{\theta}_{0},{\theta}_{3}=-\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{t}\mathrm{g}\frac{3\delta}{{\lambda}^{2}-9}+3{\theta}_{0},$
${\theta}_{21}=-\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{t}\mathrm{g}\frac{\delta \sqrt{4{\lambda}^{2}-{\delta}^{2}}+\delta}{3{\lambda}^{2}+1-{\delta}^{2}+2\sqrt{4{\lambda}^{2}-{\delta}^{2}}}+{\theta}_{0}+2{\phi}_{0},$
${\theta}_{{21}^{\text{'}}}=-\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{t}\mathrm{g}\frac{\delta \sqrt{4{\lambda}^{2}-{\delta}^{2}}-\delta}{3{\lambda}^{2}+1-{\delta}^{2}-2\sqrt{4{\lambda}^{2}-{\delta}^{2}}}-{\theta}_{0}+2{\phi}_{0},$
${\theta}_{12}=2{\theta}_{0}+{\phi}_{0},{\theta}_{12\text{'}}=-2{\theta}_{0}+{\phi}_{0},$
${\theta}_{20}=-\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{t}\mathrm{g}\frac{\delta \sqrt{4{\lambda}^{2}-{\delta}^{2}}}{-2{\lambda}^{2}+{\delta}^{2}}-\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{t}\mathrm{g}\frac{\delta \sqrt{4{\lambda}^{2}-{\delta}^{2}}}{3{\lambda}^{2}-{\delta}^{2}}+2{\phi}_{0},$
${\theta}_{30}=-\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{t}\mathrm{g}\frac{3\delta \sqrt{4{\lambda}^{2}-{\delta}^{2}}}{8{\lambda}^{2}-3{\delta}^{2}}+3{\phi}_{0},$
where $\lambda $ is rotational speed ratio, $\delta $ is damping rotational speed ratio, ${A}_{i}$ is amplitude, ${\theta}_{i}$ is phase, $\omega $ is forced vibration frequency. Eq. (16) is the general solution response of the nonlinear vibration model, including steady-state solution and transient-state solution. The first five terms ${A}_{0}\mathrm{c}\mathrm{o}\mathrm{s}\left(\mathrm{\Omega}t+{\theta}_{0}\right)+{A}_{1}\mathrm{c}\mathrm{o}\mathrm{s}\left(\mathrm{\Omega}t+{\theta}_{1}\right)\pm {A}_{2}\mathrm{c}\mathrm{o}\mathrm{s}\left(2\mathrm{\Omega}t+{\theta}_{2}\right)+{A}_{3}\mathrm{c}\mathrm{o}\mathrm{s}\left(3\mathrm{\Omega}t+{\theta}_{3}\right)\pm {A}_{4}$ are steady-state solution, and the remaining terms are transient-state solution because they contain the attenuation term ${e}^{-\mu t}$. For the ± sign of Eq. (16), the positive sign corresponds to the positive nonlinear damping, and the negative sign corresponds to the negative nonlinear damping.
4. Analytical analysis and numerical calculation
4.1. Analytical analysis
The analysis of Eq. (17) can qualitatively elaborates the influences of the coefficients (such as linear stiffness $k$, linear damping $\mu $, nonlinear stiffness ${\beta}_{0}$, nonlinear damping ${\alpha}_{0}$) on responses. In Eq. (17), all the amplitudes ${A}_{i}$, phases ${\theta}_{i}$ and forced vibration frequency $\omega $ contain rotational speed ratio $\lambda $ and damping rotational speed ratio $\delta $. $\lambda $, defined as $\lambda ={\omega}_{0}/\mathrm{\Omega}=\sqrt{k/m}/\mathrm{\Omega}$, is a function of linear stiffness $k$. $\delta $, defined as $\delta =2\mu /\mathrm{\Omega}$, is a function of linear damping $\mu $. Amplitudes ${A}_{i}$ contain nonlinear stiffness ${\beta}_{0}$ and nonlinear damping ${\alpha}_{0}$, while phases ${\theta}_{i}$ do not contain the both, and frequency $\omega $ contains only ${\beta}_{0}$. It is indicated that both nonlinear coefficients have an impact on amplitudes but not phases, and only nonlinear stiffness has an impact on forced vibration frequency. The following will be a detailed analysis of coefficients’ effects on each parameter.
1. Phases ${\theta}_{i}$ only contain $\lambda \left(k\right)$ and $\delta \left(\mu \right)$, that is, phases only depend on linear stiffness $k$ and linear damping $\mu $, irrelevant to nonlinear stiffness ${\beta}_{0}$ and nonlinear damping ${\alpha}_{0}$. Nonlinear forces do not affect the phases change.
2. Forced vibration frequency $\omega $ contains $\lambda $, $\delta $ and ${\beta}_{0}$. Compared with the natural frequency of free vibration $\omega =\sqrt{{\omega}_{0}^{2}-{\mu}^{2}}+3\epsilon {\beta}_{0}^{2}{a}_{0}^{2}/8\sqrt{{\omega}_{0}^{2}-{\mu}^{2}}$ [11, 12], the natural frequency of forced vibration is:
where the additional term $6\epsilon {\beta}_{0}^{2}{A}_{0}^{2}/8\sqrt{{\omega}_{0}^{2}-{\mu}^{2}}$ is the natural frequency correction term caused by the unbalance mass excitation of forced vibration, and this is the difference with that of free vibration model. It can be seen that, if the system adds external excitation, a corresponding correction term will be added into the natural frequency, such as the unbalanced mass excitation can cause a natural frequency correction term of nonlinear stiffness ${\beta}_{0}$ and fundamental frequency amplitude ${A}_{0}$. The nonlinear damping ${\alpha}_{0}$ has no effect on the natural frequency.
3. Amplitudes ${A}_{i}$ contain $\lambda $, $\delta $, ${\beta}_{0}$ and ${\alpha}_{0}$. It is indicated that linear and nonlinear force will affect amplitudes. The specific effects of linear damping $\delta \left(\mu \right)$ and nonlinear damping ${\alpha}_{0}$ on the steady-state and transient-state amplitudes are analyzed in detail as follows.
4. For linear damping $\delta \left(\mu \right)$, the ${\delta}^{2}$^{}term appears in the denominators of the steady-state amplitudes ${A}_{0}-{A}_{4}$, so the steady-state amplitudes decrease as the absolute value of $\delta $ increases; the ${\delta}^{2}$^{}term appears in the denominators of the transient-state amplitudes and also in the numerator of ${A}_{20}$.
If linear damping $\delta \left(\mu \right)$ is positive, the ${e}^{-\mu t}$ is an attenuation term, so the transient-state amplitudes can be ignored for they are one magnitude smaller than that of steady-state. The vibration is dominated by steady-state, and the amplitudes decrease as the absolute value of $\delta $ increases. If $\delta \left(\mu \right)$ is negative, the ${e}^{-\mu t}$ becomes a divergence term, so the transient-state amplitudes are much larger than that of steady-state. The vibration is dominated by transient-state, and the amplitudes increase as the absolute value of $\delta $ increases. That is, positive $\delta $ can restrain vibration, negative $\delta $ can intensify vibration, and this is consistent with the phenomenon in engineering practice.
Moreover, this is also consistent with the classical theory that negative damping can lead to system instability, and in this paper a more detailed mechanism analysis is given by using analytical expressions of steady-state and transient-state. If $\delta $ is positive, it can restrain both steady-state and transient-state amplitudes, so the positive linear damping can restrain the vibration, which is unquestionable. If $\delta $ is negative, though it can still restrain steady-state amplitudes, the restraint effect is much smaller than the divergence effect on transient-state amplitudes, so the total amplitudes increase rapidly and the negative linear damping can cause instability. Of course, this theoretical analysis still needs further experimental verification.
5. For nonlinear damping ${\alpha}_{0}$, the ${\alpha}_{0}^{2}$ term appears only in the numerators of the steady-state amplitudes ${A}_{2},{A}_{4}$ and the transient-state amplitudes ${A}_{40}\text{,}$${A}_{11}\text{,}$${A}_{{11}^{\text{'}}}\text{,}$${A}_{20}\text{,}$ so the amplitudes increase as the absolute value of ${\alpha}_{0}$ increases. Thus, nonlinear damping can intensify vibration. However, the amplitude change caused by nonlinear damping is smaller than that of linear damping. The quantitative comparisons of response curves are in the following Section 4.2.
Apparently, by analyzing the parameter expressions Eq. (17), it can be seen that: (1) Phase modulation is related to linear stiffness $k$ and linear damping $\mu $. (2) Frequency modulation is related to nonlinear stiffness ${\beta}_{0}$, $k$, $\mu $. (3) Amplitude modulation is related to nonlinear damping ${\alpha}_{0}$, ${\beta}_{0}$, $k$, $\mu $. All the coefficients are constants, and the machanism analysis for damping on steady-state and transient-state solution needs to be further verified by experiments.
4.2. Numerical calculation
The theoretical calculation parameters are as follows: $m=$ 1 kg, $k=$ 1×10^{4} N/m, $\mu =$ 12 rad/s, ${\omega}_{0}=$ 100 rad/s, $r=$ 1×10^{-5} m, $\eta =\epsilon {\beta}_{0}^{2}m=$ 2.5×10^{12} N/m^{3}, $\gamma =\epsilon {\alpha}_{0}^{2}m=$ 9×10^{4} Ns^{2}/m^{2}. Using MATLAB software, the influence of the coefficients on the amplitudes, phases and responses of the general solution are calculated and analyzed.
Fig. 2 shows the total response of Fig. 2(a), steady-state response of Fig. 2(b) and transient-state response of Fig. 2(c, d) under the conditions of positive nonlinear damping (left) and negative nonlinear damping (right), where ${x}_{total}$, ${x}_{steady}$, ${x}_{transient}$ represent the general solution (i.e., the sum of steady-state and transient-state), steady-state solution and transient-state solution respectively, ${f}_{i}$ represents the response of corresponding term ${A}_{i}\times \mathrm{c}\mathrm{o}\mathrm{s}\left(\right)$ in which the subscript is the same. For the steady-state response, it is obviously that ${f}_{3}$ is the smaller one in Fig. 2(b). For the transient-state response, ${f}_{a0}$, ${f}_{11}$, ${f}_{{11}^{\text{'}}}$, ${f}_{{12}^{\text{'}}}$ are the larger ones in Fig. 2(c) while the rest smaller ones in Fig. 2(d). The transient-state response is about two orders of magnitude smaller than the steady-state response.
Fig. 2. a) The total, b) steady-state and c), d) transient-state response (left: positive nonlinear damping, right: negative nonlinear damping)
a)
b)
c)
d)
Compared the two responses under positive and negative nonlinear damping conditions, it can be seen that the response curves which are almost symmetrical about $x$-axis also have the difference of positive and negative, while the positive and negative of the displacement responses are only represent the difference in vibration direction, and the trend of curves remains the same.
Fig. 3 shows the amplitudes ${A}_{i}$ and phases ${\theta}_{i}$ of steady-state of Fig. 3(a) and transient-state of Fig. 3(b) under positive nonlinear damping condition. The curves under negative nonlinear damping condition have the same trend, so omitted. As can be seen from the figure, both amplitudes and phases change with rotational speed and have peak value at the critical speed ($\mathrm{\Omega}={\omega}_{0}=$ 100 rad/s), indicating that at critical point the amplitudes will reach extremum and the phases change are the integer multiple of 90 deg. The phases ${\theta}_{2}$, ${\theta}_{3}$ of steady-state have peak value at 1/2, 1/3 of speed range additionally, which is caused by the second and third harmonic generations.
Fig. 4 shows the influence of positive value and negative value change of the linear damping coefficient of Fig. 4(a) and nonlinear damping coefficient of Fig. 4(b) on responses, where the nonlinear damping coefficient is written in the form of $\gamma =\epsilon {\alpha}_{0}^{2}m$ for calculation convenience. Note that only negative linear damping coefficient $\mu $ corresponds to the steady-state response, while others correspond to the total response (i.e., the sum of steady-state and transient-state). It can be seen from Fig. 4 that the vibration decreases as the absolute value of linear damping increases (only steady-state vibration decreases as $\mu $ is negative), and the vibration increases as the absolute value of nonlinear damping increases. If the damping is doubled, the vibration change is less than twice caused by nonlinear damping, and about one order of magnitude caused by linear damping, so nonlinear damping’s impact is less than linear damping’s. If $\mu $ is negative, the total response is given in Fig. 5. It is obvious that the magnitude of transient-state response increases rapidly as the attenuation term becomes divergence term, and the steady-state response is negligible.
Fig. 3. a) Steady-state and b) transient-state
a)
b)
Fig. 4. a) Responses of linear and b) nonlinear damping coefficient (Only negative linear damping’s response is steady-state response, others are total responses.)
a)
b)
Fig. 5. Transient-state response is much larger than steady-state for negative linear damping
Therefore, it can be indicated that: (1) Positive linear damping can restrain vibration. (2) Negative linear damping’s restraint effect on steady-state is far less than the exponential divergence effect on transient-state, so it can intensify vibration. (3) Nonlinear damping can intensify vibration, but its influence is less than that of linear damping.
4.3. Experiment results
Set up a single-disc principle rotor test bed for verification, as shown in Fig. 6. Rotor material with $\phi $235 steel, disc mass 0.5 kg, disc radius 0.039 m, shaft radius 0.005 m, elastic modulus 200 GPa, Poisson ratio 0.3, bearing span 0.4335 m, disc distance from the left bearing of 0.13005 m. Displacement and phase signals are input to the DASP software for real-time monitoring and analysis [18, 19]. The experiment on the effect of the oil film force of sliding bearing is carried out to verify the damping effect on the vibration behaviors in the nonlinear vibration model proposed in this paper.
Fig. 7 shows the time-amplitude-frequency spectrum of Fig. 7(a), and the bifurcation diagram of speed-up process of Fig. 7(b). In Fig. 7(a), the rotor is first speed-up then speed-down, the maximum speed of 10237 rpm, where Res is resonance area, Sub1 is oil whirl area, Sub2 is oil whip area and Sub3 is double subharmonic area. In Fig. 7(b), typical shaft orbits are selected for each bifurcation section, namely, period-1 (5110 rpm), bifurcation (6878 rpm) and chaos (9739 rpm), which is shown in detail in Fig. 8.
Fig. 6. Single-disc principle rotor test bed
Fig. 7. a) Time-amplitude-frequency spectrum and b) bifurcation diagram of speed-up process
a)
b)
Fig. 8. Shaft orbits
For the vibration behaviors such as rotation and whirl in the experiment, the oil film force model of reference [20] is used to analyze the damping effect, which is also corresponding to the linear damping term $2\mu \dot{x}$ of Eq. (8). The oil film force model [20] is ${F}_{\pi}={f}_{1}{\omega}_{0}-{f}_{2}\mathrm{\Omega}+{f}_{3}\dot{\epsilon}$, where ${f}_{i}$ is positive coefficients, ${\omega}_{0}$, $\mathrm{\Omega}$, $\dot{\epsilon}$ is rotational frequency, whirl frequency and squeezing rate respectively. As can be seen from Fig. 8, the amplitude of shaft orbit increases from period-1 (5110 rpm) to bifurcation (6878 rpm) due to the negative linear damping effect which can reduce the loading capacity of the oil film. It is consistent with the theoretical analysis in this paper.
When the rotational speed is above 9000 rpm, the development path of chaos can be seen clearly from the figure, and the amplitude of shaft orbit increases gradually, such as 9739 rpm. It is because under the combined effects of negative linear damping and nonlinear damping, the subharmonic frequency of self-excited vibration generated by shaft and oil film is close to the natural frequency of the rotor system, which causes the frequency locking resonance, namely the oil film whip phenomenon. When the rotational speed increases in the whip area, the rotational frequency increases, so the positive linear damping increases, which lead to the amplitude decrease; the whirl frequency increases, so the negative linear damping increases, which lead to the amplitude increase. Therefore, the oil film whirl and whip characteristics of the experiment shows that the positive linear damping can restrain vibration, while the negative linear damping and the nonlinear damping can intensify vibration.
5. Conclusions
In this paper, by introducing linear and nonlinear force, a nonlinear vibration model is established for high-speed flexible rotor system, and the analytical solutions are deduced and analyzed. The conclusions are as follows.
1. The general form of nonlinear vibration model is constructed, and the analytical solutions of steady-state and transient-state are derived under the conditions of positive and negative nonlinear damping coefficient.
2. The influence and physical significance of the coefficients on parameters are analyzed. The phase modulation is related to linear stiffness $k$ and linear damping $\mu $. The frequency modulation is related to nonlinear stiffness ${\beta}_{0}$, $k$, $\mu $. The amplitude modulation is related to nonlinear damping ${\alpha}_{0}$, ${\beta}_{0}$, $k$, $\mu $.
3. Positive linear damping can restrain vibration, while negative linear damping and nonlinear damping can intensify vibration. Nonlinear damping has less influence on vibration than linear damping. The experimental results verify the damping effect and the rationality of the nonlinear vibration model.
Acknowledgements
The work is supported by National Science and Technology Projects (Grant No. 2012BAA11B02), and the support is gratefully acknowledged.
References
- Newkirk B. L. Shaft whipping. General Electric Review, Vol. 27, 1924, p. 169-178. [Search CrossRef]
- Kimball A. L. Internal friction as a cause of shaft whirling. Philosophical Magazine, Vol. 49, Issue 292, 1925, p. 724-727. [Publisher]
- Timoshenko S. Vibration Problem in Engineering. Second Edition, D. van Nostrand Company, New York, 1937. [Search CrossRef]
- Gunter E. J. The influence of internal friction on the stability of high speed rotors. Journal of Engineering for Industry, Vol. 89, Issue 4, 1967, p. 683-688. [Publisher]
- Ehrich F. F. Shaft whirl induced by rotor internal damping. Journal of Applied Mechanics, Vol. 31, Issue 2, 1964, p. 279-282. [Publisher]
- Crandall S. H. Heuristic Explanation of Journal Bearing Instability. Rotordynamic Instability Problems in High-Performance Turbomachinery, NASA CP-2250, 1982. [Search CrossRef]
- Muszynska A. Whirl and whip-rotor/bearing stability problems. Journal of Sound and Vibration, Vol. 110, Issue 3, 1986, p. 443-462. [Publisher]
- Tasker F., Chopra I. Nonlinear damping estimation from rotor stability data using time and frequency domain techniques. AIAA Journal, Vol. 30, Issue 5, 1992, p. 1383-1391. [Publisher]
- ChandraN. H., SekharA. S. Nonlinear damping identification in rotors using wavelet transform. Mechanism and Machine Theory, Vol. 100, 2016, p. 170-183. [Publisher]
- ZhaoJ., LinnettI., McleanL. Imbalance response and stability of eccentric squeeze-film-damped nonlinear rotor bearing systems. JSME International Journal Series B Fluids and Thermal Engineering, Vol. 37, Issue 4, 1994, p. 886-895. [Publisher]
- Wu M., Han D., Yang S., et al. Research on nonlinear vibration analytical model of rotor system. Journal of Aerospace Power, 2016, (in Press, in Chinese). [Search CrossRef]
- Wu M., Yang S., Han D., et al. Analytical model for nonlinear vibration of flexible rotor system. Journal of Vibroengineering, Vol. 18, Issue 8, 2016, p. 4980-4994. [Publisher]
- Liu B., Peng J. Nonlinear Dynamics. Higher Education Press, Beijing, 2004, (in Chinese). [Search CrossRef]
- Khalil H. Nonlinear Systems. Electronic Industry Press, Beijing, 2005, (in Chinese). [Search CrossRef]
- Zhou J., Zhu Y. Nonlinear Vibration. Xi’an Jiaotong University Press, Xi’an, 1998, (in Chinese). [Search CrossRef]
- Hu H. Nonlinear Vibration. Aviation Industry Press, Beijing, 2000, (in Chinese). [Search CrossRef]
- Chen Y. Modern Analytical Method in Nonlinear Dynamics. Science Press, Beijing, 1992, (in Chinese). [Search CrossRef]
- Yang J., Wu M. Single-disc rotor nonlinear vibration model and mechanism. Journal of Aerospace Power, Vol. 27, Issue 8, 2012, p. 1738-1745, (in Chinese). [Search CrossRef]
- Wu M., Yang J. The research on single-disc rotor nonlinear vibration model and mechanism. Mechatronics and Information Technology, Parts 1 and 2, Vol. 2, Issue 3, 2012, p. 954-959. [Search CrossRef]
- Yang J., Chen C., Yang S., et al. An analytic model of oil-film force on hydrodynamic journal bearing of finite length. ASME Turbo Expo 2008: Power for Land, Sea and Air, Parts A and B, 2008, p. 947-952. [Publisher]