To accurately predict temperature field changes of the subsurface of the workpiece, a heat source model of the contact zone in the real grinding process is established using the temperature matching method according to the principle of anti-heat source and the actually measured grinding temperature. The finite element method is utilized to simulate the change of the grinding temperature field of the workpiece, and make a comparison with the triangular heat source model. The grinding temperature field at different depths of the workpiece in the grinding zone is measured using the K-type thermocouple. The results indicate that the heat source model according to the temperature matching method has a good consistency with the measured values at different depths on the subsurface. The relative errors are between 1.8% and 8%. Comparing with the triangular heat source model, the distribution accuracy of the predicted temperature field is improved almost 2 times.