The evolution of the phase transformation and the resulting internal stresses and strains in me- tallic parts during quenching were modeled numerically. The numerical simulation of the metal quenching process was based on the metallo-thermo-mechanical theory using the finite element method to couple the temperature, phase transformation, and stress-strain fields. The numerical models are presented for the heat treatment and kinetics of the phase transformation. The finite element models and the phase transition kinetics accurately predict the distribution of the microstructure volume fractions, the temperature, the distor- tion, and the stress-strain relation during quenching. The two examples used to validate the models are the quenching of a small gear and of a large turbine rotor. The simulation results for the martensite phase vol- ume fraction, the stresses, and the distortion in the gear agree well with the experimental data. The models can be used to optimize the quenching conditions to ensure product quality. The evolution of the phase transformation and the resulting internal stresses and strains in me tali parts during quenching were modeled numerically. The numerical simulation of the metal quenching process based on the metallo-thermo-mechanical theory using the finite element method to couple the The numerical models are presented for the heat treatment and kinetics of the phase transformation. The finite element models and the phase transition kinetics precisely predict the distribution of the microstructure volume fractions, the temperature, the The two examples used to validate the models are the quenching of a small gear and of a large turbine rotor. The simulation results for the martensite phase vol ume fraction, the stresses, and the distortion in the gear agree well with the experimental data. The models can be used to optimize the quench ing conditions to ensure product quality.