A new method of reconstructing constitutive equations from microscopic simulations of kinetic models for polymeric solutions is presented. Microscopic simulations provide a set of data from which the constitutive equations are extracted (learned). In order to be consistent with non-equilibrium thermodynamics, we work within the Quasi-Equilibrium approximation.