The bifurcation analysis of bearing number in hybrid porous gas-lubricated bearing(HPGB)system are analyzed by a hybrid numerical method combining the finite difference method and differential transformation method.The numerical results are obtained and verified by the traditional finite difference method and the current analytical results are found to be in good agreement.Furthermore,the results reveal the changes which take place in the dynamic behavior of the bearing system as the bearing numbers are increased.From the dynamic responses of the rotor center,they reveal complex dynamic behaviors including periodic,sub-harmonic motion and chaos.The proposed method and results provide an effective means of gaining insights into the nonlinear dynamics of HPGB systems.