In this paper the complex rheological behaviors of rod-like polymer solutions in plane Couette flow are studied using micro-macro simulation. The micro-macro model is described by Doi theory with local inhomogeneous fluid. The finite volume method is used to solve the conservation equation on the macroscopic level and the Smoluchowski equation on the microscopic level in the numerical simulation. The numerical results predict several plane flow modes including some typical plane modes and two new complicated modes. Furthermore, simulation results show that the rheological phase diagrams for solutions of rigid rod molecules are mainly dependent on the Deborah number, the concentration of solutions and molecular scale, and the varying Deborah number have great influence on some micro-properties of tumbling period and orientation angle of flow-aligning.