By using the DAMASK simulation package we developed a new approach to predict the evolution of anisotropic yield functions by coupling large scale forming simulations directly with crystal plasticity-spectral based virtual experiments, realizing a multi-scale model for metal forming.