With the consideration of slip deformation mechanism and various slip systems of body centered cubic (BCC) metals, Taylor-type and finite element polycrystal models were embedded into the commercial finite element code ABAQUS to realize crystal plasticity finite element modeling, based on the rate dependent crystal constitutive equations. Initial orientations measured by electron backscatter diffraction (EBSD) were directly input into the crystal plasticity finite element model to simulate the development of rolling texture of interstitial-free steel (IF steel) at various reductions. The modeled results show a good agreement with the experimental results. With increasing reduction, the predicted and experimental rolling textures tend to sharper, and the results simulated by the Taylor-type model are stronger than those simulated by finite element model.'Conclusions are obtained that rolling textures calculated with 48 {110} <111>+ {112} <111>+ {123} <111> slip systems are more approximate to EBSD results.