A semimicroscopic derivation is presented of equations of motion for the density and the flow velocity of concentrated systems of entangled polymers. The essential ingredient is the transient force that results from perturbations of overlapping polymers due to flow. A Smoluchowski equation is derived that includes these transient forces. From this, an equation of motion for the polymer number density is obtained, in which body forces couple the evolution of the polymer density to the local velocity field. Using a semimicroscopic Ansatz for the dynamics of the number of entanglements between overlapping polymers, and for the perturbations of the pair-correlation function due to flow, body forces are calculated for nonuniform systems where the density as well as the shear rate varies with position. Explicit expressions are derived for the shear viscosity and normal forces, as well as for nonlocal contributions to the body force, such as the shear-curvature viscosity. A contribution to the equation of motion for the density is found that describes mass transport due to spatial variation of the shear rate. The two coupled equations of motion for the density and flow velocity predict flow instabilities that will be discussed in more detail in a forthcoming publication.