In this paper a general solution method is presented to obtain the unbalance response orbit from the finite element based equations of motion of a gear-coupled two-shaft rotor-bearing system, whose shafts rotate at their different speeds from each other. Particularly, are proposed analytical solutions of the maximum and minimum radii of the orbit. The method has been applied to analyze the unbalance response of a 800 refrigeration-ton turbo-chiller rotor-bearing system having a bull-pinion speed increasing gear. Bumps in the unbalance response of the driven high speed compressor rotor system have been observed at the first torsional natural frequency due to the coupling effect of lateral and torsional dynamics. Further, the proposed analytical solutions have agreed well with those obtained by a full numerical approach. The proposed analytical solutions can be generally applied to obtain the maximum and minimum radii of the unbalance response orbits of dual-shaft rotor-bearing systems coupled by bearings as well.