Abstract In this paper, an efficient numerical method consisting of the real mode component mode synthesis (CMS) model reduction, shooting method with parallel computing, and Floquet analysis was developed for nonlinear rotordynamics analysis of a flexible rotor supported by a four-lobe flexure pivot tilting pad journal bearing (FPTPJB) in load-on-pad (LOP) and load-between-pad (LBP) orientations in comparison to a fixed profile journal bearing (JB) of the same pad geometry. The method used the rotor's finite elements and bearing forces obtained from directly solving the Reynolds equation to determine the limit cycles and Hopf bifurcation types. For the investigated rotor and bearing parameters, the numerical results indicated that the onset speed of instability (OSI) of FPTPJB is considerably higher than that of JB of the same orientation. Also, FPTPJB in LOP orientation yielded substantially higher OSI than the LBP one, whereas the OSI of JB in LOBP orientation was higher than the LOP counterpart. Nonlinear calculation results indicated that all bearing types and orientations gave subcritical Hopf bifurcation. The FPTPJB in LOP orientation produced the largest stable operating region, whereas the JB in LBP configuration yielded the smallest one. The experiment showed subcritical Hopf bifurcation occurred at speed close to the calculated OSI in all cases except FPTPJB in LOP orientation that the OSI is higher than the maximum test rig speed. The whirling orbit had the same frequency as the first critical speed and was precessed in the shaft rotation direction.