In this work, we proposed a virtual laboratory based on full-field crystal plasticity simulation to track plastic anisotropy and to calibrate yield functions for multi-phase metals. The virtual laboratory, minimally, only requires easily accessible EBSD data for constructing the high-resolved microstructural representative volume element and macroscopic flow stress data for identifying the micromechanical parameters of constituent phases. An inverse simulation method based on global optimization scheme was developed for parameters identification, and a nonlinear least-squares method was employed to calibrate the yield functions. Various mechanical tests of an advanced high strengthening steel (DP980) sheet under different loading conditions were conducted to validate the virtual laboratory. Three well-known yield functions, the quadratic Hill48, Yld91, and Yld2004-18p, were selected as the validation benchmarks. All the studied functions, calibrated by numerous stress points under arbitrary loading conditions, successfully captured both the deformation and strength anisotropies. Furthermore, the full-field CP modeling well correlates the microscopic deformation mechanism and plastic heterogeneity to the macro-mechanical behavior of the sheet. The proposed virtual laboratory, which is readily extended with physically based CP model, could be a versatile tool to explore and predict the mechanical property and plastic anisotropy of advanced multi-phase metals.