We present results of molecular dynamics simulations of linear polyelectrolytes in solution. The fundamental model for polyelectrolytes in solution is studied, Specifically, simulations are performed for multichain systems of a flexible chain model of charged polymers. The full Coulomb interactions of the monomers and counterions are treated explicitly. Experimental measurements of the osmotic pressure and the structure factor are reproduced. The simulations reveal a new picture of the chain structure based on calculations of the structure factor, persistence length, end-to-end distance, etc. We present a detailed discussion of the chain structure and a comparison with present theories. In contrast to the predicted dilute limit of rodlike chains, we find that the chains have significant bending at very low densities. Furthermore, the chains contract significantly before they overlap. We also show that counterion condensation dramatically alters the chain structure.