Among other critical conditions in rotor systems the large non-linear vibration excited by bearing non-linearity causes the rotor failure. For reducing this catastrophic failure and predictive detection of this phenomenon the analysis of orthotropic bearing non-linearity in rotor system using higher order frequency response functions (HFRFs) is conducted and is shown to be theoretically feasible as that of non-rotating structures. The complex HFRFs based on the Volterra series are newly developed for the process and investigated their features by using the simple forms of the FRFs associated with the forward and the backward modes.