Laser quenching has a short process cycle and high production efficiency, and it plays an important role in automobile, ship, machinery manufacturing, and other fields. The surface hardness of 40Cr steel laser quenching parts is uneven by unreasonable set on the quenching overlap zone. This will affect the microstructure of the quenched layer, resulting in cracking, corrosion during service, and other hazards that ultimately reduce reliability. Numerical simulations provide an effective way to quantitatively reveal the transient evolution of the multi-field coupling between temperature field, stress field, and phase transition field in quenching, which directly determines the extent of the overlap zone and quenching properties. The quenched phase transition layer profile is predicted to effectively determine the extent of the secondary tempering softening zone and solve the bottleneck problem of uneven surface hardness in quenching. The innovation of this paper is to establish a multi-field coupled numerical model of the 40Cr steel multi-track laser quenching process. The transient quenching temperature, phase transition hardening, and stress distribution were numerically calculated, and the size of the tempering zone under different overlapping rates was evaluated, revealing the internal coupling mechanism and correlation between multi-fields during the laser quenching. The quenching temperature, microstructure, and hardness distribution of 40Cr steel were tested by an infrared thermometer, Axio Vert A1 Zeiss microscope, Thermo ScientificTM Apreo scanning electron microscope, and Q10M microhardness tester, which verified the effectiveness of numerical simulation. The research can provide an important theoretical basis for optimizing quenching process parameters in production.