Most of ferrous b.c.c weld materials may experience martensitic transformation during rapid cooling after welding. It is well known that volume expansion due to the phase transformation could influence on the relaxation of welding residual stress. To apply this effect practically, it is a prerequisite to establish a numerical model which is able to estimate the effect of phase transformation on residual stress relaxation quantitatively. For this purpose, the analysis is carried out in two regions. i.e., heating and cooling, because the variation of material properties following a phase transformation in cooling is different in comparison with the case in heating, even at the same temperature. The variation of material properties following phase transformation is considered by the adjustment of specific heat and thermal expansion coefficient, and the distribution of residual stress in analysis is compared with that of experiment by previous study. consequently, in this study, simplified numerical procedures considering phase transformation, which based on a commercial finite element package was established through comparing with the experimental data of residual stress distribution by other researcher. To consider the phase transformation effect on residual stress relaxation, the transition of mechanical and thermal property such as thermal expansion coefficient and specific heat capacity was found by try and error method in this analysis.