Coarse-grained molecular dynamic simulations are performed for lyotropic liquid-crystalline (LC) solutions of semiflexible rod-like molecules, which are realized by harmonic-stretching potential and angle-bending potential. The phase behavior and molecular dynamic information are investigated. When the concentration of rod-like molecules increases, the solution transits from isotropic to binary phase and then to anisotropic phase. The influences of chain length, temperature, solvent solubility, and chain rigidity are studied. The simulation for the phase behavior successfully reproduces the predictions of Flory lattice theory on the rod-like molecule solution. The decrement of the chain rigidity increases the concentration of phase transition, and the apparent persistence length is increases by molecular alignments with LC formations at higher concentrations.