The purpose of the present work is to introduce a crystallization law into Rem3D, a 3D code written in C++ and dedicated to the injection molding of polymers. We kept the basic hypotheses of Avrami's model and cast the kinetic equations into a differential system that is solved numerically. The variation of the density of potential nuclei with temperature is taken into account. Furthermore, the distribution of mean spherulite sizes can be deduced from the calculations. The second part of the paper is an experimental study of crystallization in well-controlled conditions (2D, isothermal or constant cooling-rate). It establishes a procedure for the determination of the nucleation and growth parameters used in the theoretical model, and gives a first validation of this model. Finally, the crystallization equations are introduced into Rem3D, in order to assess the feasibility of our new approach. Some typical results concerning the evolution of the transformed volume fraction in injection-molded parts are presented.