A method for solving the coupling relationship between temperature, stress and phase-transformation in quenching process was presented, and the FEM software for quenching process was developed. End-quenching of P20 steel, 2-D transient heat transfer problem without latent heat and quenching process of 1080 steel were simulated using this software by established several suitable FEM models. The numerical analysis values of temperature and the experimental results of phase-transformation volumes, hardness distribution and stress were contrasted with FEM simulation results. The contrast between the experimental results and the simulation ones shows that the simulation results of FEM software are consistent with the experimental ones or the numerical analysis values.