Research Article  Open Access
Qi Xu, Yuanqing Luo, Hongliang Yao, Lichao Zhao, Bangchun Wen, "Eliminating the FluidInduced Vibration and Improving the Stability of the Rotor/Seal System Using the InerterBased Dynamic Vibration Absorber", Shock and Vibration, vol. 2019, Article ID 1746563, 9 pages, 2019. https://doi.org/10.1155/2019/1746563
Eliminating the FluidInduced Vibration and Improving the Stability of the Rotor/Seal System Using the InerterBased Dynamic Vibration Absorber
Abstract
A theoretical research on eliminating the instability vibration and improving the stability of the rotor/seal system using the inerterbased dynamic vibration absorber (IDVA) is presented in this paper. The modified Jeffcott rotor and Muszynska nonlinear seal force models are employed. The proposed IDVA is the damping element of the traditional dynamic vibration absorber (DVA) replaced by one of the six configurations of the inerter. The instability threshold speed of the system is obtained by applying the Routh–Hurwitz stability criterion. The quantum particle swarm optimization (QPSO) method is utilized to optimize the parameters of the IDVA. The numerical method is applied to investigate the nonlinear dynamic responses and stability. The results show that the IDVA can effectively improve the stability and reduce the instability vibration of the rotor/seal system. Furthermore, the performance of the IDVA is more effective than that of the traditional DVA.
1. Introduction
The nonlinear seal force between the rotor and seal is prone to the instability vibration in the rotating machine. The instability vibration is one kind of the selfexcited vibrations and adversely causes the amplitude to increase dramatically. Hence, the instability vibration must be controlled.
Great efforts for the instability vibration control had been made by Bently and Muszynska [1]. The nonlinear behavior and stability had been analysed via the Bently/Muszynska (B/M) nonlinear seal force model. The expression of the instability threshold was defined, and the “cross stiffness” component was the main reason for the instability. The fluid circumferential average velocity ratio is considered a key parameter for the instability vibration control. It can be realized using the various seals to disrupt the fluid flow [2–4]. The antiswirl flow in the seal was a typical method to eliminate the instability [5]. It was widely used in engineering. Zhang et al. [6] computed the dynamic and flow characteristics of the honeycomb seal and corresponding labyrinth seal in the rotor, and the results showed that the instability vibration could be reduced due to the additional aerodynamic force generated by the seal. Then, they presented two types of the labyrinth seals in the rotor system: the straightthrough and corresponding interlaced labyrinth seal. The dynamic and flow characteristics of the rotor/seal system were analysed by the numerical method, and it showed that the straightthrough labyrinth seal could reduce the rotor instability vibration [7]. Jia et al. [8] investigated the dynamic leakage rate of the rotorlabyrinth seal system with the variable rotating speed. The similar conclusion that the labyrinth seal could improve the rotor stability was obtained. Zhang et al. [9] studied dynamic characteristics of the rotor/bearing/labyrinthseal system. The research indicated that the shorter seal length could improve the stability of the rotor system, but the sealing effect was decreased.
Recently, the DVA has begun to attract attention to eliminate the rotor vibration. Hu and He [10] designed a rotor DVA to online control the rotor critical speed vibration and discussed the effect of the installation position on the suppression performance. Heidari and Monjezi [11] investigated the vibration control of the unbalanced rotor using the virtual passive DVA with the optimal parameters. The permanent magnet DVAs with the negative and tunable stiffness were presented to reduce the rotor unbalance vibration by Yao et al. [12, 13]. Besides, the nonlinear energy sink (NES, a kind of the nonlinear stiffness DVA) and centrifugal pendulum vibration absorber (CPVA, a kind of the torsion DVA) have been applied for the antiresonance and torsional vibration suppression of the rotor, respectively. For example, Taghipour et al. [14] studied the antiresonance for the nonlinear bearing of the rotor with the DVA and NES. Tehrani and Dardel [15] investigated the DVA and NES to prevent the contact between the rotor and stator. The NES was also employed for the vibration control of the unbalanced hollow rotor [16] and continuous rotorbliskjournal bearing system [17]. Shi et al. [18] and Nishimura [19] et al. described the torsional vibration reduction in the rotor using the CPVA.
Now that the inerter, which is proposed for the vibration control of the vehicle suspension [20], has already been applied in the DVA, called the IDVA. Hu and Chen [21] investigated the parameter optimization to improve the performance of the IDVA. The extended fixedpoint theory was used for the parameter optimization of the IDVA [22, 23]. Lazarek et al. [24] designed a tuned mass damper with the variable inerter. The IDVA has been utilized to the vibration control of the cable [25, 26], the offshore wind turbine [27], tall building [28], and so on.
The traditional DVA could eliminate the instability vibration in the rotor/seal system [29, 30]. The fluidinduced instability elimination of the rotor with the IDVA is proposed in this paper. First, the model of the rotor/seal system with the IDVA is established. Second, the expression of the instability threshold speed is obtained. Third, the numerical calculation is employed for the nonlinear dynamic and stability analysis. Then, the performance comparison between the IDVA and traditional DVA is discussed. Finally, some conclusions are drawn to summarize the effects of the fluidinduced vibration elimination and stability region extension on the rotor/seal system by the IDVA.
2. Model of Rotor/Seal System with IDVA
Figure 1shows the model of the rotor/seal system with the IDVA. The instability vibration usually happens in the turbomachinery, such as the centrifugal compressor. Normally, the running speed of the kind of the rotor system is between the first and secondorder critical speeds. Therefore, the modified Jeffcott rotor model and B/M nonlinear seal force model are employed for the rotor/seal system model [1, 29]. The rotor model is employed for considering the firstorder vibration shape. The nonlinear dynamic behaviors and stability of the real rotor/seal system can be described by the models. The IDVA is added in the rotor in the perpendicular directions by the rolling bearing in order to avoid rotation with the rotor [10]. The DVA and IDVA are depicted in Figure 2. The damping element of the DVA is replaced by one of the six configurations of the inerter: c_{1}, c_{2}, c_{3}, c_{4}, c_{5}, and c_{6} [21]. Taking c_{3} as an example in the main part of this paper, the equations of motion for the rotor/seal system with the IDVA are as follows:where m, d, and k are the rotor generalized mass, lateral external damping, and lateral isotropic stiffness, respectively. m_{e} and r denote, respectively, the eccentric mass and eccentricity. ω is the rotating speed. m_{a} and k_{a} denote, respectively, the mass and stiffness coefficients of the IDVA. x, y, x_{a}, and y_{a} are the horizontal and vertical displacement generalized coordinates of the rotor and IDVA, respectively. f_{xSeal} and f_{ySeal} are the components of the nonlinear seal force in the x and y directions, respectively. The expressions arewhere m_{f}, d_{f}, and k_{f} denote, respectively, the seal fluid inertia, seal fluid film radial damping, and seal fluid stiffness coefficients. λ is the seal fluid circumferential average velocity ratio. The expressions of the nonlinear terms arewhere d_{0}, k_{0}, λ_{0}, n, and b are the parameters of the nonlinear seal force. These parameters of the nonlinear force are described in detail by Bently and Muszynska [1, 31]. r_{f} is the radial clearance of the seal.
(a)
(b)
f_{xIDVA} and f_{yIDVA} denote, respectively, the force components of the IDVA that acted on the rotor in the x and y directions. The expressions are as follows:where b, c, and k_{2} are the inertance, damping, and stiffness coefficients of the configuration of the inerter. x_{c} and y_{c} denote, respectively, the generalized displacements between the inertance and damping element in the x and y directions. x_{k} and y_{k} denote, respectively, the generalized displacements between the inertance and stiffness element in the x and y directions.
To define the dimensionless parameters:
The dimensionless transform is as following:and define the generalized coordinates (Taking the x coordinate as an example):
Conveniently, the dimensionless generalized coordinates are still expressed as x, y, x_{a}, y_{a}, x_{c}, y_{c}, x_{k}, and y_{k}, respectively. The original equations become the dimensionless equations:where
The dimensionless components of the nonlinear seal force:where
The dimensionless force components of the IDVA acted on the rotor:where
3. Stability Analysis
The unbalance force is assumed zero, the nonlinear terms and seal fluid inertia coefficient are neglected as they are minimal for the stability analysis [32]. The original dimensionless equations (8) can be rewritten:
The characteristic equation can be obtained:where μ is the eigenvalue and q is the degree of freedom of equation (14). a_{0}, a_{1}, a_{2}, … , a_{2q} are the coefficients of the characteristic equation.
According to the RouthHurwitz stability criterion [33], the instability threshold speed ω_{i} of the rotor system can be obtained from the following expressions:where H is the Hurwitz determinant, and:
The rotor/seal system is stable as ω < ω_{i} and unstable as ω > ω_{i}.
4. Numerical Simulations and Discussion
The numerical calculation is employed due to the complicated nonlinear terms of equation (8) and implicit equation (16). The parameters of the rotor/seal system are obtained from reference [29] in order to compare with the previous research. The rotor m = 1 kg, ζ = 0.025, ω_{r} = 200 rad/s, r_{f} = 10^{−3} m and ε_{e} = 0.01. The nonlinear seal force ε_{f} = 10^{−6}, d_{0} = 100 Nm/s, k_{0} = 1 N/m, λ_{0} = 0.48, n = 2 and b = 0.5. Assuming that the dimensionless rotating speed s = ω/ω_{r}, the frequency ratio σ_{a} = ω_{a}/ω_{r} and σ_{b} = ω_{b}/ω_{r}. The QPSO algorithm is utilized to obtain the optimal IDVA parameters. The target function is the rotor instability threshold speed ω_{i}. The mass ratio ε_{a} = 0.01 and ε_{b} is a constant. Moreover, the damping ratio ζ_{a}, dimensionless frequency ratio σ_{a} and σ_{b} are the optimization parameter. Correspondingly, the parameter optimization ranges are 0 <ζ_{a} < 1, 0 <σ_{a} ≤ 2 and 0 <σ_{b} ≤ 2, respectively. The flow chart of the QPSO method is shown in Figure 3.
The optimization results are shown in Table 1. It can be seen that the maximum of the instability threshold speed is ω_{i} = 2.759 as ε_{b} = 0.03, ζ_{a} = 0.05, σ_{a} = 0.98, and σ_{b} = 0.95. The convergence diagram of the optimization as ε_{b} = 0.03 is in Figure 4. It shows that the optimization is convergent as the number of the iterations is greater than 15.

Therefore, the parameters of the IDVA are selected as follows: ε_{b} = 0.03, ζ_{a} = 0.05, σ_{a} = 0.98, and σ_{b} = 0.95. The changes of the coefficients in equation (15) and Hurwitz determinants to estimate the stability are presented in Figure 5. It can be seen that a_{0}, a_{1}, a_{2}, … , a_{12} are positive (Figure 5(a)); H_{1}, H_{2}, … , H_{8} are also positive (Figure 5(b)), H_{9} and H_{10} are also positive at s ≤ 2.76, and H_{11} is equal to zero when s is between 2.755 and 2.76 (Figure 5(c)). The coefficients and Hurwitz determinants satisfy the Routh–Hurwitz stability criterion. The bifurcation diagram, amplitudefrequency response curve, frequency spectrum, orbit, and Poincaré map are used to investigate the nonlinear behavior by the Newmark time integration method [34]. The algorithm is unconditionally convergent as the values of parameters β ≥ 1/4 and δ ≥ 1/2 [35].
(a)
(b)
(c)
The bifurcation diagram and amplitudefrequency response curve of the rotor/seal system in the x coordinate are shown in Figure 6. The blue and red solid dots and lines are the bifurcation and amplitudefrequency response curve of the rotor/seal system without and with the IDVA, respectively. It can be observed in Figure 6(a) that the vibration of the rotor/seal system is synchronous to one period at s < 2.31. Only one point is correspondingly shown in the bifurcation diagram for every rotating speed. Then, the vibration changes and bifurcates at s = 2.31. The vibration is doubleperiodic, multiperiodic, or quasiperiodic at s ≥ 2.31 and can be determined by using the Poincaré map. s = 2.31 is the instability threshold speed of the rotor/seal system. The variances of the vibration in the rotor/seal system as adding the IDVA are as follows: the vibration of the system is synchronous to one period at s < 2.805; then the vibration changes and bifurcates at s = 2.805; the vibration is doubleperiodic, multiperiodic, or quasiperiodic at s ≥ 2.805 and can be determined by using the Poincaré map; and s = 2.805 is the instability threshold speed of the rotor/seal system. Comparing with the rotor/seal system without the IDVA, the instability threshold speed increases by about 21.43%. The instability vibration between s = 2.31 and s = 2.805 becomes the synchronous vibration. The instability vibration is eliminated completely in this range. Simultaneously, the instability threshold speed increases by about 1.67% comparing with Table 1. The result illustrates that the instability threshold obtained by the linear equation (16) is less than that obtained by the nonlinear equation (8). Whether or not to add the IDVA, the amplitude of the system increases rapidly when the rotating speed exceeds the instability threshold speed as shown in Figure 6(b). But the instability amplitude can be reduced by the IDVA at s ≥ 2.805. The instability vibration is eliminated partially in this range.
(a)
(b)
The speed s = 2, s = 2.5, and s = 3 are picked, respectively, for the dynamic behavior analysis in the range of s < 2.31, 2.31 < s < 2.805, and s > 2.805 as shown in Figures 7–9. The blue and red solid dots and lines are the response of the rotor/seal system without and with the IDVA, respectively. Figure 7 presents the response is synchronous to one period, and there is less change after adding the IDVA as the speed s = 2 is away from the resonance region. The response is fiveperiodic when s = 2.5 as shown in Figure 8. The instability vibration frequency is close to the firstorder natural frequency of the rotor. The amplitude of the instability frequency is much higher than that of the rotating frequency. Then the response becomes one period, and the instability frequency disappears after adding the IDVA. As the speed continues to increase, Figure 9 represents the response is threeperiodic when s = 3. The response becomes quasiperiod, and the instability frequency is eliminated partially after adding the IDVA.
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
5. Comparison with the Previous Work
The traditional DVA has been applied to eliminate the instability vibration in the rotor/seal system in reference [29]. The comparison results of the instability threshold speed are shown in Table 2. It can be seen that the instability threshold speed of the rotor with the IDVA is larger than that of the rotor with the traditional DVA except the configuration of the inerter c_{2}. The maximum instability threshold speed is obtained from the configuration of the inerter c_{3}.

The comparison results of the bifurcation diagram and amplitudefrequency response curve are shown in Figure 10. The black and red solid dots and lines are the response of the rotor/seal system with the traditional DVA and IDVA, respectively. It can been seen that the instability threshold speed of the rotor with the IDVA (c_{3}) is larger than that of the rotor with the traditional DVA, and the instability amplitude of the rotor with the IDVA (c_{3}) is less than that of the rotor with the traditional DVA. The result still illustrates that the linear model appears more conservative in evaluating instability.
(a)
(b)
The comparison results of the relationship between the instability threshold speed and frequency ratio σ_{a} and the instability threshold speed and damping ratio ζ_{a} are shown in Figure 11. The black dashed line represents the instability threshold speed of the rotor. Comparing the rotor with the traditional DVA, Figure 11(a) shows that the range of increase in instability threshold speed is almost constant as the frequency ratio σ_{a} changes, while Figure 11(b) shows that the range of increase in instability threshold speed is decreased as the damping ratio ζ_{a} changes. And the maximum of the instability threshold speed increases by about 3.13%.
(a)
(b)
6. Conclusions
The fluidinduced vibration elimination and stability margin extension of the rotor/seal system using the IDVA is proposed in this paper. The rotor/seal system dynamic model is established, and the six configurations of the IDVA are added, respectively, in the system model. The instability threshold speed of the system is obtained by applying the Routh–Hurwitz stability criterion. The QPSO algorithm is utilized to obtain the optimal IDVA parameters. The numerical method is applied to investigate the nonlinear dynamic responses and stability of the rotor/seal system without and with the IDVA. The results are compared with the previous work. The conclusions are drawn as follows.(1)The IDVA can eliminate the fluidinduced instability vibration and improve the instability threshold speed of the rotor/seal system.(2)Comparing with the traditional DVA, the effect of the IDVA on the instability vibration elimination and stability improvement is better, but the range of increase in the instability threshold speed is decreased as the damping ratio ζ_{a} changes.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors would like to gratefully acknowledge the National Natural Science Foundation of China (grant no. U1708257) and the Fundamental Research Funds for Central Universities (no. N180313009) for the financial support for this study.
References
 D. E. Bently and A. Muszynska, Rotordynamics, CRC Press, Taylor & Francis Group, Boca Raton, USA, 2005.
 D. Childs, D. Elrod, and K. Hale, “Annular honeycomb seals: test results for leakage and rotordynamic coefficients; comparisons to labyrinth and smooth configurations,” Journal of Tribology, vol. 111, no. 2, pp. 293–301, 1989. View at: Publisher Site  Google Scholar
 P. Varney and I. Green, “Steadystate response of a flexibly mounted stator mechanical face seal subject to dynamic forcing of a flexible rotor,” Journal of Tribology, vol. 139, no. 6, Article ID 062201, 2017. View at: Publisher Site  Google Scholar
 M. Arghir, N. ManhHung, D. Tonon, and J. Dehouve, “Analytic modeling of floating ring annular seals,” Journal of Engineering for Gas Turbines and PowerTransactions of the ASME, vol. 134, no. 5, Article ID 052507, 2012. View at: Publisher Site  Google Scholar
 D. Sun, S. Wang, and Y. Ai, “Theoretical and experimental research on the performance of antiswirl flow for the static and dynamic characteristics of seals,” Journal of Mechanical Engineering, vol. 52, no. 3, pp. 101–109, 2016. View at: Publisher Site  Google Scholar
 H. Zhang, Q. Zheng, G. Yue, and J. Gao, “Numerical analysis of flows and aerodynamic forces in honeycomb and labyrinth seals,” Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 227, no. 9, pp. 1965–1979, 2013. View at: Publisher Site  Google Scholar
 H. Zhang, X. Jia, X. Pan, B. Jiang, and Q. Zheng, “Interaction between rotor and annular seals: interlaced and straightthrough labyrinth seals,” Journal of Propulsion and Power, vol. 32, no. 6, pp. 1483–1493, 2016. View at: Publisher Site  Google Scholar
 X. Jia, H. Zhang, Q. Zheng, S. Fan, and Z. Tian, “Investigation on rotorlabyrinth seal system with variable rotating speed,” International Journal of Turbo & JetEngines, vol. 36, no. 1, pp. 19–29, 2019. View at: Publisher Site  Google Scholar
 E. Zhang, Y. Jiao, and Z. Chen, “Dynamic behavior analysis of a rotor system based on a nonlinear labyrinthseal forces model,” Journal of Computational and Nonlinear Dynamics, vol. 13, no. 10, Article ID 101002, 2018. View at: Publisher Site  Google Scholar
 H.L. Hu and L.D. He, “Online control of critical speed vibrations of a singlespan rotor by a rotor dynamic vibration absorber at different installation positions,” Journal of Mechanical Science and Technology, vol. 31, no. 5, pp. 2075–2081, 2017. View at: Publisher Site  Google Scholar
 H. Heidari and B. Monjezi, “Vibration control of imbalanced Jeffcott rotor by virtual passive dynamic absorber with optimal parameter values,” Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 232, no. 23, pp. 4278–4288, 2018. View at: Publisher Site  Google Scholar
 H. Yao, Z. Chen, and B. Wen, “Dynamic vibration absorber with negative stiffness for rotor system,” Shock and Vibration, vol. 2016, Article ID 5231704, 13 pages, 2016. View at: Publisher Site  Google Scholar
 H. Yao, T. Wang, B. Wen, and B. Qiu, “A tunable dynamic vibration absorber for unbalanced rotor system,” Journal of Mechanical Science and Technology, vol. 32, no. 4, pp. 1519–1528, 2018. View at: Publisher Site  Google Scholar
 J. Taghipour, M. Dardel, and M. H. Pashaei, “Vibration mitigation of a nonlinear rotor system with linear and nonlinear vibration absorbers,” Mechanism and Machine Theory, vol. 128, pp. 586–615, 2018. View at: Publisher Site  Google Scholar
 G. G. Tehrani and M. Dardel, “Mitigation of nonlinear oscillations of a Jeffcott rotor system with an optimized damper and nonlinear energy sink,” International Journal of NonLinear Mechanics, vol. 98, pp. 122–136, 2018. View at: Publisher Site  Google Scholar
 C. Guo, M. A. AlShudeifat, A. F. Vakakis, L. A. Bergman, D. M. McFarland, and J. Yan, “Vibration reduction in unbalanced hollow rotor systems with nonlinear energy sinks,” Nonlinear Dynamics, vol. 79, no. 1, pp. 527–538, 2015. View at: Publisher Site  Google Scholar
 S. Bab, S. E. Khadem, M. Shahgholi, and A. Abbasi, “Vibration attenuation of a continuous rotorbliskjournal bearing system employing smooth nonlinear energy sinks,” Mechanical Systems and Signal Processing, vol. 84, pp. 128–157, 2017. View at: Publisher Site  Google Scholar
 C. Shi, S. W. Shaw, and R. G. Parker, “Vibration reduction in a tilting rotor using centrifugal pendulum vibration absorbers,” Journal of Sound and Vibration, vol. 385, pp. 55–68, 2016. View at: Publisher Site  Google Scholar
 K. Nishimura, T. Ikeda, and Y. Harata, “Localization phenomena in torsional rotating shaft systems with multiple centrifugal pendulum vibration absorbers,” Nonlinear Dynamics, vol. 83, no. 3, pp. 1705–1726, 2016. View at: Publisher Site  Google Scholar
 Y. Hu, M. Z. Q. Chen, and Y. Sun, “Comfortoriented vehicle suspension design with skyhook inerter configuration,” Journal of Sound and Vibration, vol. 405, pp. 34–47, 2017. View at: Publisher Site  Google Scholar
 Y. Hu and M. Z. Q. Chen, “Performance evaluation for inerterbased dynamic vibration absorbers,” International Journal of Mechanical Sciences, vol. 99, pp. 297–307, 2015. View at: Publisher Site  Google Scholar
 E. Barredo, A. Blanco, J. Colín et al., “Closedform solutions for the optimal design of inerterbased dynamic vibration absorbers,” International Journal of Mechanical Sciences, vol. 144, pp. 41–53, 2018. View at: Publisher Site  Google Scholar
 X. Wang, X. Liu, Y. Shan, Y. Shen, and T. He, “Analysis and optimization of the novel inerterbased dynamic vibration absorbers,” IEEE Access, vol. 6, pp. 33169–33182, 2018. View at: Publisher Site  Google Scholar
 M. Lazarek, P. Brzeski, and P. Perlikowski, “Design and identification of parameters of tuned mass damper with inerter which enables changes of inertance,” Mechanism and Machine Theory, vol. 119, pp. 161–173, 2018. View at: Publisher Site  Google Scholar
 X. Shi and S. Zhu, “Dynamic characteristics of stay cables with inerter dampers,” Journal of Sound and Vibration, vol. 423, pp. 287–305, 2018. View at: Publisher Site  Google Scholar
 J. Luo, J. Z. Jiang, and J. H. G. Macdonald, “Cable vibration suppression with InerterBased absorbers,” Journal of Engineering Mechanics, vol. 145, no. 2, Article ID 04018134, 2019. View at: Publisher Site  Google Scholar
 Y. Hu, J. Wang, M. Z. Q. Chen, Z. Li, and Y. Sun, “Load mitigation for a bargetype floating offshore wind turbine via inerterbased passive structural control,” Engineering Structures, vol. 177, pp. 198–209, 2018. View at: Publisher Site  Google Scholar
 A. Giaralis and F. Petrini, “Windinduced vibration mitigation in tall buildings using the tuned massdamperinerter,” Journal of Structural Engineering, vol. 143, no. 9, Article ID 04017127, 2017. View at: Publisher Site  Google Scholar
 Q. Xu, J. Niu, H. Yao, L. Zhao, and B. Wen, “Fluidinduced vibration elimination of a rotor/seal system with the dynamic vibration absorber,” Shock and Vibration, vol. 2018, Article ID 1738941, 15 pages, 2018. View at: Publisher Site  Google Scholar
 Q. Xu, J. Niu, H. Yao, L. Zhao, and B. Wen, “Nonlinear dynamic behavior and stability of a rotor/seal system with the dynamic vibration absorber,” Advances in Mechanical Engineering, vol. 11, no. 1, Article ID 1687814018819578, 2019. View at: Publisher Site  Google Scholar
 J. Hua, S. Swaddiwudhipong, Z. S. Liu, and Q. Y. Xu, “Numerical analysis of nonlinear rotor—seal system,” Journal of Sound and Vibration, vol. 283, no. 3–5, pp. 525–542, 2005. View at: Publisher Site  Google Scholar
 C.C. Fan and M.C. Pan, “Active elimination of oil and dry whips in a rotating machine with an electromagnetic actuator,” International Journal of Mechanical Sciences, vol. 53, no. 2, pp. 126–134, 2011. View at: Publisher Site  Google Scholar
 Y. Huang, Z. Tian, R. Chen, and H. Cao, “A simpler method to calculate instability threshold speed of hydrodynamic journal bearings,” Mechanism and Machine Theory, vol. 108, pp. 209–216, 2017. View at: Publisher Site  Google Scholar
 W. Kim, J. Y. Lee, and J. Chung, “Dynamic analysis for a planetary gear with timevarying pressure angles and contact ratios,” Journal of Sound and Vibration, vol. 331, no. 4, pp. 883–901, 2012. View at: Publisher Site  Google Scholar
 H. Ma, H. Li, H. Niu, R. Song, and B. Wen, “Numerical and experimental analysis of the firstand secondmode instability in a rotorbearing system,” Archive of Applied Mechanics, vol. 84, no. 4, pp. 519–541, 2014. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Qi Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.