We present and discuss results obtained by an extensive nonequilibrium molecular dynamics computer simulation study of polymer solutions under shear, where the chain consists of N beads connected by a finitely extendable nonlinear elastic (FENE) spring force and the solvent is explicitly taken into account. Various scaling laws are extracted from the data which allow one to predict the qualitative-and to certain extent also quantitative-structural and rheological behavior of polymer solutions under good solvent conditions. For most quantities, the results drawn from simulation are compared with experimental data and theoretical predictions which are based on similar models (e.g., harmonic bond potentials or Brownian dynamics methods). Specifically, and in contrast to common theoretical approaches, the simulation yield information about a set of different but characteristic relaxation times, which determine the rheological and structural behavior (for example, flow birefringence, structure factor, rotational dynamics) separately, the difference either resulting from the underlying static or dynamic nature or from relaxation processes which act on different length scales.