This study is inspired by the recent development of “virtual material testing laboratory” in which the main equipment are full field crystal plasticity models. Ample examples have demonstrated its applications to sheet forming operations. In those applications, the mechanical anisotropy originated from the crystallographic texture can be adequately described, such as r-values and earing. Formability is also another very important property in sheet metal forming, which yet have not been equipped in these virtual laboratories. Though theoretical models for formability can be dated back to 1800s, all popular models at the moment suffer respective limitations. In this study, we make an attempt to the prediction of formability limit diagram by crystal plasticity finite element method, with the objective of overcoming the limitations posed by the current theoretical models. The focus is placed on the texture dependence of the formability limit diagram.