In the last decades simulation tools for Brownian dynamics of polymers have attracted more and more interest. Here we present a mathematically consistent finite element approach to the simulation of Brownian polymer dynamics. The viscous damping forces are accounted for by an anisotropic friction model. By comparison with theoretical predictions and experimental data we demonstrate the reliability and efficiency of this method.