Based on the unsteady Reynolds equation with Reynolds boundary, two-dimensional (2D) Reynolds equation is transformed into one-dimensional (1D) by taking the assumption of parabolic pressure distribution in axial direction in oil film field. Finite difference method was employed to solve 1D Reynolds equation, and the approximate pressure distribution was obtained in oil film field. Dynamic behaviors of a flexible rotor system with fixed-tilting pad journal bearings support were analyzed while the inertia of the pads was taken into consideration in the model. Imbalance responses of a symmetrical rotor-combination journal bearings (fixed-tilting pad journal bearings) system were investigated using Poincaré map and self-adaptive Runge-Kutta method. Numerical results reveal rich and complex nonlinear phenomena, such as periodic, quasi-periodic motion, etc.