This paper is the first part of the study on the accuracy improvement of numerical analysis of steel welds. The aim of this paper is to raise the accuracy of thermal analysis results, such as the shape and size of the weld cross section and the hardness distribution in HAZ(Heat-Affected Zone). It is known that the factors affecting on the accuracy are thermal properties, metallurgical properties and welding heat source model. It was found that the arbitrary distributed heat source model should be used to predict practical weld cross section shape and size. Also, in order to improve the prediction accuracy of HAZ hardness distribution, it was essential to consider 2 CCT(Continuous Cooling Transformation) diagrams in calculating volume fraction of transformed phases. One is the peak temperature being around melting temperature. The other is the peak temperature being around metallurgical transformation temperature.