The present study was an attempt for systematic data conversion between FDM and FEM in order to evaluate the thermal stress distribution during quenching process. It has been generally recognized that FDM is efficient in flow and temperature analysis and FEM in that of stress. But it induced difficulty and tedious work in analysis that one uses both FDM and FEM to take their advantages because of the discrepancy of nodes between analysis tools. So we proposed field data conversion procedure from FDM to FEM in 3-dimensional space, then applied this procedure to analysis of quenching process. The simulation procedure calculates the distributions of temperature and microstructure using FDM and microstructure evolution equations of diffusion and diffusionless transformation. FEM was used for predicting the distributions of thermal stress. The present numerical code includes coupled temperaturephase transformation kinetics and temperature-microstructure dependent material properties. Calculated results were compared with previous experimental data to verify the method, which showed good agreements.