Theoretical and bifurcation analysis of a flexible rotor supported by gas-lubricated bearing system with porous bushing

Cheng-Chi Wang1 , Jui-Pin Hung2

1, 2Graduate Institute of Precision Manufacturing, National Chin-Yi University of Technology, Taichung, Taiwan

1Corresponding author

Journal of Vibroengineering, Vol. 18, Issue 3, 2016, p. 1934-1940. https://doi.org/10.21595/jve.2015.16752
Received 6 October 2015; received in revised form 9 December 2015; accepted 20 December 2015; published 15 May 2016

Copyright © 2016 JVE International Ltd. 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.
Creative Commons License
Abstract.

The bifurcation behaviors of a flexible rotor supported by gas-lubricated bearing system with porous bushing are analyzed by a novel numerical method combining the finite difference method and differential transformation method. The results obtained by proposed method are verified with traditional finite difference method, and the analytical results by these two methods are consistent and also in good agreement. Furthermore, the dynamic orbits, power spectra, bifurcation diagram, Poincaré maps and maximum Lyapunov exponent are used to confirm the changes of rotor behavior as the rotor mass is increased. The results show that the rotor center reveals complex dynamic behaviors including periodic, sub-harmonic and chaotic motions. Especially, the rotor bearing system behaves chaos over the ranges of rotor mass 10.7 ms< 13.8 kg. The current results provide an effective means of the gas bearing systems and further understanding of the nonlinear dynamic behavior of bearing systems characterized by different rotor masses.

Keywords: bifurcation, porous bushing, chaotic.

1. Introduction

Gas-lubricated bearing system with porous bushing (GBSP) is different from other air bearing system due its particularly attractive features such as low-noise rotation, high load capacity, improved damping properties and zero friction. GBSP is simpler in construction and cheaper than the externally pressurized bearings. The comparison with oil bearings and GBSP reveals that GBSP generates less heat and provides higher accuracy. GBSP are also applied in high speed spindles, and machine tools.

In 1964, Sneck and Yen [1] considered one-dimensional flow for radial direction in a porous medium and developed a perturbation solution for a finite journal bearing. Then, they verified theoretical solutions experimentally in 1965 [2].

In 2001, Sinha et al. [3] analyzed externally pressurized conical bearings with the porous constant gap when the slider rotated with uniform angular velocity. The system governing equations including coupled momentum and energy equations are analyzed by using finite difference method and various bearing characteristics are also determined. The results show that the inlet pressure decreases remarkably for highly porous surfaces, resulting in reduced load capacity of the bearing. But, the torque of system with respect to variation of permeability remains unaffected.

For the rotor dynamic study, the nonlinear dynamic behavior of rigid rotor supported by noncircular air journal bearing system including two, three and four-lobe bearings are solved by Rashidi et al. [4, 5] The Reynolds’ equation for air pressure distribution is analyzed by using finite element method and the rotor dynamic equations are calculated by Runge-Kutta method. The numerical results show that bifurcation phenomenon with different bearing numbers and rotor masses including periodic and non-periodic motions at different operational situations. It’s found that bearing number and rotor mass are the major parameters for the bearing system.

In 2010 to 2013 [6-9], Wang et al. studied various kinds of gas bearing and analyzed the bifurcation and dynamics behaviors avoiding the occurrence of fatigue or unexpected situations. For example, as united gas-lubricated bearings, the rotor center behaves a complex dynamic behavior under different values of the rotor mass and bearing number. The maximum Lyapunov exponent is also applied to distinguish different motions comprising periodic, sub-harmonic, quasi-periodic and chaotic responses. Furthermore, the dynamic behavior of noncircular air bearing system is studied by Wang et al. and a hybrid method including differential transformation method and finite difference method are proposed and used to solve the governing equations with various gas viscosity and rotational speed.

The remainder of this study is organized as follows. Section 2 develops a mathematical model of GBSP. Section 3 determines the Reynolds’ equation to obtain the gas pressure distribution by a hybrid method combining the finite difference method (FDM) and the differential transformation method (DTM). The numerical solutions are then compared with those obtained by traditional FDM. Section 4 presents the simulation results obtained using the proposed hybrid method for the vibrations of the rotor center for various rotor masses.

2. Governing equations

The GBSP system is used to form a closed cavity around the porous bushing with an appropriate casing and the porous material acts as a restrictor. The compressed gas is supplied and fed at a constant pressure, P^f, and the gas flows through the porous bushing into the bearing clearance. As the supply pressure drops to P^, it exhausts into the atmosphere at pressure Pa. The GBSP model incorporates the ideal assumptions including the gas flow is assumed isothermal, the side flow of gas of the bearing is neglected, and the gas viscosity is constant. It is considered that a GBSP system comprising a perfectly-balanced flexible rotor supported symmetrically on two identical GBSPs mounting in turn on rigid pedestals. Due to the rotor is assumed as balancing perfectly and the GBSP is symmetric about its central axes, the system is confined to a single bearing supporting a rotor with two degrees of translatory oscillation in the transverse plane.

The gas pressure distribution between the rotor and the bushing is modeled by the non-dimensionalized Reynolds’ equation as follows:

(1)
2 ς Λ β P ^ f τ - P ^ f θ ϕ τ = 2 P ^ f z ^ 2 ,
(2)
Λ ϕ 2 P ^ f 2 z ^ 2 z ^ = 1 + 2 Λ 1 - 2 ς ϕ τ P ^ h ^ θ + 4 ς Λ P ^ h ^ τ
            = θ h ^ 3 P ^ 2 θ + D L 2 h ^ 3 2 P ^ 2 z ^ 2 ,

where P^ is the non-dimensionalized pressure, P^f is the non-dimensionalized pressure in porous media, h^ is the non-dimensionalized thickness of gas film, Λ is the bearing number and θ, z^ are the non-dimensionalized coordinates.

3. Simulation analysis by a novel numerical method

The simulation analysis is completed by two different numerical methods including the proposed method in this paper and traditional FDM. Accordingly, the FDM is applied for Eqs. (1) and (2) to discretize the θ and z^ directions by the central-difference scheme and then the time domain τ by the implicit-back-difference scheme. Also, a uniform mesh size is used. In order to compare the calculation results and increase the accuracy of system, a novel numerical method is proposed in this study. It is commenced by applying the DTM to discretize the Reynolds’ equation given in Eqs. (1) and (2) with respect to time and thus Eqs. (1) and (2) become as Eqs. (3)-(6). Then, The FDM is used to discretize Eqs. (3)-(6) with respect to the θ and z^ directions. Note that Eqs. (3)-(6) are discretized using the second-order-accurate central-difference scheme for both the first and the second derivatives. Substituting Eqs. (5)-(6) into Eqs. (3)-(4) yields Eqs. (7), (8):

(3)
2 ς Λ β P ^ f τ - 2 ς Λ β ϕ τ P ^ f θ = 2 P ^ f z ^ 2 ,
(4)
Λ ϕ 2 P ^ f 2 z ^ 2 z ^ = 1 + 2 Λ P ^ θ h ^ + 2 Λ h ^ θ P ^ - 4 Λ ς ϕ τ P ^ θ h ^
            - 4 ς Λ ϕ τ h ^ θ P ^ + 4 Λ ς P ^ τ h ^ + 4 ς Λ h ^ τ P ^
          = J θ Q θ + J 2 Q 2 θ + D L 2 J 2 Q z ^ 2 ,

where:

(5)
J k = h ^ 3 = h ^ h ^ h ^ = l = 0 k h ^ i , j k - l m = 0 l h ^ i , j l - m h ^ i , j m ,
(6)
Q k = P ^ 2 = P ^ P ^ = l = 0 k P ^ i , j k - l P ^ i , j l .

The governing equations become:

(7)
2 ς Λ β k + 1 H ~ P ^ f i , j k + 1
            - 2 ς Λ β l = 0 k k + 1 H ~ P ^ f i + 1 , j k - l - P ^ f i - 1 , j k - l 2 Δ θ ϕ i , j l + 1
            = P ^ f i , j + 1 ( m ) - 2 P ^ f i , j ( m ) + P ^ f i , j - 1 ( m ) Δ z ^ 2 ,
(8)
Λ ϕ P ^ f 2 i , j + 1 m - 2 P ^ f 2 i , j m + P ^ f 2 i , j - 1 m Δ z ^ 2 z ^ = 1
            + 2 Λ l = 0 k h ^ i , j k - l P ^ i + 1 , j l - P ^ i - 1 , j l 2 Δ θ + 2 Λ l = 0 k P ^ i , j k - l h ^ i + 1 , j l - h ^ i - 1 , j l 2 Δ θ
            - 4 Λ ς l = 0 k h ^ i , j k - l m = 0 l k + 1 H ~ P ^ i + 1 , j m - P ^ i - 1 , j m 2 Δ θ ϕ i , j l - m
            - 4 Λ ς l = 0 k P ^ i , j k - l m = 0 l k + 1 H ~ h ^ i + 1 , j m - h ^ i - 1 , j m 2 Δ θ ϕ i , j l - m
            - 4 Λ ς l = 0 k k + 1 H ~ h ^ i , j k - l P ^ i , j l + 1 + 4 Λ ς l = 0 k k + 1 H ~ P ^ i , j k - l h ^ i , j l + 1
            = l = 0 k J i + 1 , j k - l - J i - 1 , j k - l 2 Δ θ Q i + 1 , j l - Q i - 1 , j l 2 Δ θ
            + l = 0 k J i , j k - l Q i , j + 1 l - 2 Q i , j l + Q i , j - 1 l Δ I 2
            + D L 2 l = 0 k J i , j k - l Q i , j + 1 l - 2 Q i , j l + Q i , j - 1 l Δ z ^ 2 .

The calculation of iterative procedure can be integrated as follows. First, the new values of acceleration, velocity, and displacement of rotor are obtained by following a time increment Δτ. Second, from first step the rotor center displacements can then be determined and the corresponding change of the gas film gap (h^) can be analyzed. Then the new value of h^ is substituted into governing dequations to give the new pressure distribution of the gap between the shaft and the bushing. Third, the pressure distribution obtained from second step is integrated to calculate the internal force caused by the gas-lubricated film. Finally, the displacement and velocity values calculated in the first step, the gas pressure computed in the second step, and the internal force analyzed in the third step are taken as the new initial conditions. Applying this new set of conditions, the calculation procedure returns to the first step to calculate the changes in the GBSP system during the time interval Δτ2Δτ. These data including the orbital paths and velocity of the rotor center are used to generate bifurcation diagrams and Poincaré maps.

4. Numerical results

Table 1 presents the Poincaré maps obtained by the FDM and DTM&FDM for the orbits of the rotor center. It is observed that a good agreement by DTM&FDM method exists between the two sets of results at different rotor mass values. It also compares with different values of the time step. It can be seen that the orbits are in agreement to approximately 4 decimal places for the different time steps especially obtained by DTM&FDM method as shown in Table 1 and 2.

Table 1. Comparison of rotor center orbits calculated by FDM and DTM&FDM methods, respectively

Displacement conditions
X2(nT)
Time step: 0.001
Time step: 0.01
FDM
m r =   7.2 kg
ω =   1980 rad/s
0.008823820432
0.008961092741
DTM&FDM
0.008825114313
0.008824802704
FDM
m r =   10.1 kg
ω =   1980 rad/s
0.059034858544
0.058062840042
DTM&FDM
0.059136319514
0.059136084996

Table 2. Poincaré maps with different time increments and rotor mass by DTM&FDM

m r = 7.2 kg
m r = 10.1 kg
τ
X2(nT)
Y2(nT)
τ
X2(nT)
Y2(nT)
π /300
0.008823943
–0.060185357
π /300
0.059136649
–0.075228092
π /600
0.008823325
–0.060185087
π /600
0.059136640
–0.075228387

Fig. 1. Dynamic orbits of rotor center

 Dynamic orbits of rotor center

a) mr= 9.7 kg

 Dynamic orbits of rotor center

b)mr= 9.96 kg

 Dynamic orbits of rotor center

c)mr= 10.8 kg

 Dynamic orbits of rotor center

d)mr= 13.5 kg

Figs. 1 and 2 show that the orbits and power spectra of the rotor center are regular at low values of the rotor mass (mr= 9.7 and 9.96 kg), but become irregular at mr= 10.8 kg. This irregular behavior diverge as non-symmetric and non-periodic motion. When the rotor mass is increased at mr= 13.5 kg, the system converges to periodic motion.

Fig. 2. Power spectra of rotor center

 Power spectra of rotor center

a)mr= 9.7 kg

 Power spectra of rotor center

b)mr= 9.96 kg

 Power spectra of rotor center

d)mr= 10.8 kg

 Power spectra of rotor center

e)mr= 13.5 kg

Fig. 3. Bifurcation diagrams versus rotor mass

 Bifurcation diagrams versus rotor mass

a) X2(nT)

 Bifurcation diagrams versus rotor mass

b) Y2(nT)

Fig. 3 presents the bifurcation diagrams and shows the rotor center displacement of GBSP against the rotor mass. It is observed that qualitatively different behavior is varied at different values of mr within the range 6.0 to 14.0 kg. Fig. 4 present the Poincaré maps at mr= 9.7, 9.96, 10.8, and 13.5 kg, respectively. The dynamic motion of the rotor center is T-periodic at lower values of the rotor mass, i.e. mr< 9.7 kg. However, the periodic motion loses its stability at mr= 9.7 kg and is replaced by sub-harmonic motion of 2T as shown in Fig. 3 and Fig. 4(a). This sub-harmonic is persisted over the interval 9.7 mr< 9.96, but when the mass is increased to mr= 9.96 kg, the sub-harmonic motion is changed to a T-periodic motion in the x- and y-directions shown in Fig. 3 and Fig. 4(b).

When the mass is increased to mr= 10.8 kg, system is changed to chaotic motion and maintained over the interval 10.8 mr< 13.5 kg. The chaotic situation is shown in Fig. 3 and proved in Fig. 4(c). However, the chaotic motion transferred its motion at a rotor mass of 13.5 kg and is bifurcated into a 2T-periodic motion shown in Fig. 3 and Fig. 4(d). Then, the sub-harmonic of 2T motion is behaved over the interval 13.5 mr 14.0 kg. It can be seen that two discrete points in the Poincaré maps at mr= 9.7 and 13.5 kg mean 2T sub-harmonic motion and non-order discrete points at mr= 10.8 kg reveal chaos.

In order to further verify the occurrence of chaotic motion, the maximum Lyapunov exponent is applied and it is shown that a positive value (as mr equals 10.8) presents chaos over the intervals 10.8 mr< 13.5 kg shown in Fig. 5(a) and equals zero which means system behaves periodic motion. From the discussions above, it is evident that the behavior of the rotor center is dependent on the rotor mass. Table 3 summarizes the motions performed by the rotor center for rotor mass values in the range 6.0 mr 14.0 kg.

Fig. 4. Poincaré maps of rotor center trajectories

 Poincaré maps of rotor center trajectories

a)mr= 9.7 kg

 Poincaré maps of rotor center trajectories

b)mr= 9.96 kg

 Poincaré maps of rotor center trajectories

c)mr= 10.8 kg

 Poincaré maps of rotor center trajectories

d)mr= 13.5 kg

Fig. 5. Maximum Lyapunov exponents of system at different values of rotor mass

 Maximum Lyapunov exponents of system at different values of rotor mass

a) mr= 10.8 kg

 Maximum Lyapunov exponents of system at different values of rotor mass

b)mr= 9.7 kg

Table 3. Variation of rotor center response with rotor mass over interval 6.0 mr 14.0 kg

Rotor mass
[6.0, 9.7)
[9.7, 9.96)
[9.96, 10.8)
[10.8, 13.5)
[13.5, 14.0)
Behavior
T
2T
T
Chaotic
2T

5. Conclusions

In this paper, the bifurcation behaviors of a flexible rotor supported by gas-lubricated bearing system with porous bushing (GBSP) are analyzed and the chaotic behavior has also been found by utilizing a hybrid numerical scheme comprising the differential transformation method and the finite difference method in this paper. The system orbits, power spectra, bifurcation diagrams, Poincaré maps and maximum Lyapunov exponent have revealed the presence of a complex dynamic behavior comprising periodic, sub-harmonic and chaotic responses of the rotor center. From Table 1 the results obtained by the FDM and proposed hybrid method for the orbits of the rotor center prove that a good agreement exists between different sets of results. For the results of nonlinear dynamic motion, this paper provide a fully understanding of GBSP systems characterized by different rotor masses. Specifically, the results have shown that system exists chaotic motion over the ranges of rotor mass 10.8 mr< 13.5 kg under rotational speed of 1980 rad/s.

Acknowledgements

The financial support of this research by Ministry of Science and Technology, under the Projects Nos. MOST 103-2622-E-167-015-CC3 and 103-2221-E-167-011 are greatly appreciated.

References

  1. Sneck H. J., Yen K. T. The externally pressurized, porous wall, gas-lubricated journal bearing. Part 1. ASLE Transactions, Vol. 7, 1964, p. 288-298. [Publisher]
  2. Sneck H. J., Elwell R. C. The externally pressurized, porous wall, gas-lubricated journal bearing. Part 2. ASLE Transactions, Vol. 8, 1965, p. 339-345. [Search CrossRef]
  3. Sinha P., Chandra P., Bhartiya S. S. Thermal effects in externally pressurized porous conical bearings with variable viscosity. ACTA Mechanica, Vol. 149, 2001, p. 215-227. [Publisher]
  4. Rashidi R., Karami Mohammadi A., Bakhtiarinejad F. Preload effect on nonlinear dynamic behavior of a rigid rotor supported by noncircular gas-lubricated journal bearing systems (three and four lobe). Nonlinear Dynamics, Vol. 60, 2009, p. 231-253. [Publisher]
  5. Rashidi R., Karami Mohammadi A., Bakhtiarinejad F. Bifurcation and nonlinear dynamic analysis of a rigid rotor supported by two-lobe noncircular gas-lubricated journal bearing system. Nonlinear Dynamics, Vol. 61, 2010, p. 783-802. [Publisher]
  6. Wang C. C., Kuo C. L. Nonlinear dynamic analysis of a relatively short spherical gas journal bearing system. Journal of Mechanical Science and Technology, Vol. 24, Issue 8, 2010, p. 1565-1571. [Publisher]
  7. Wang C. C. Bifurcation and nonlinear dynamic analysis of united gas-lubricated bearing system. Computers and Mathematics with Applications, Vol. 64, Issue 5, 2012, p. 729-738. [Publisher]
  8. Wang C. C., Wang C. C. Bifurcation and nonlinear dynamic analysis of noncircular aerodynamic journal bearing system. Nonlinear Dynamics, Vol. 72, Issue 1, 2013, p. 477-489. [Publisher]
  9. Wang C. C. Bifurcation analysis of bearing number in ultra-short gas bearing system. Smart Science, Vol. 1, Issue 1, 2013, p. 18-24. [Search CrossRef]

Cited By

Artificial Intelligence Review
Aneesh G. Nath, Sandeep S. Udmale, Sanjay Kumar Singh
2021
Petroleum Engineering
E.A. Kozhukhov, E.A. Petrovsky, K.A. Bashmur, O.A. Kolenchukov
2020
Journal of the Brazilian Society of Mechanical Sciences and Engineering
Chenhui Jia, Zhiwu Cui, Shijun Guo, Ming Qiu, Wensuo Ma
2019
Journal of Vibroengineering
Xiulong Chen, Yonghao Jia
2018
Advances in Mechanical Engineering
Cheng-Chi Wang, Tsui-Er Lee
2017
Journal of Vibroengineering
Rui Zhu, Yanru Zhang, Jianxing Ren, Hongguang Li, Qingkai Han
2017