It is very common for aircraft engines to have dual rotor or even triple rotor designs. Due to the complexity of having multiple rotor design, the transfer matrix methods have used in the past to deal with multiple rotor-bearing systems. However, due to transfer matrix method’s assumptions, sometimes resulted in numerical stability problems or root-missing problems. The purpose of this paper is to develop a systematic theoretical analysis of the dynamic characteristics of turbomachinery dual rotor-bearing systems. This dual rotor-bearing system analysis will start with a finite element (FEM) rotor-bearing system dynamic model, then using different methods to verify the analysis results including critical speed map and bearing stiffness. In an inertia coordinate system, a general model of continuous dual rotor-bearing systems is established based on a lagrangian formulation. Gyroscopic moment, rotary inertia, bending and shear deformations have been included in the model. From a point of view of the systematic approach, a solution of the finite element method is used to calculate the critical speeds by several different methods, which in turn can help to verify this dual rotor-bearing system approach. The effects of the speed ratio of dual rotors on the critical speed will be studied, which in turn can be used as one of the dual rotor design parameters. Also, both critical speeds are in effect functions of dual rotor speeds. Finally, the bearing stiffness between high speed and low speed shafts not only affect the critical speeds of the dual rotor system, but also affect the mode shapes of the system. Therefore, the bearing stiffness in between is of even greater importance in turbomachinery dual rotor or multiple rotor design.