This study is a numerical analysis of the industrial quenching process for leaf springs developed by the CAVEO company. The leaf chosen for this study is of a parabolic profile made of EN-51CrV4 steel (AISI 6150). The aim of this study is to set up a numerical model to predict thermal, metallurgical, and mechanical behavior of a leaf spring from exit of the heating furnace to exit of the quenching bath going through a cambering operation. This study would therefore allow the company to switch from a development scheme based on experiments using physical prototypes tested on the production line to a new scheme based on virtual prototypes using numerical simulation. The development of the numerical model using the finite element method is carried out using the ABAQUS/Implicit solver coupled with two user subroutines Phase and UMAT. The first one have been developed to compute microstructure evolution and the second one to define the constitutive law taking into account phase transformations. This model helps us to follow the spatio-temporal evolutions of temperature and microstructure in the leaf, as well as the variation of the leaf deflection during the process. The proposed numerical model is supported by an experimental protocol based on infrared thermographic images, Rockwell-C hardness measurements, metallographic observations, and deflection measurements. Indeed, the results of the proposed thermo-mechanical and metallurgical model are closed to experimental results.