Finite element analysis of welding processes, which entail phase evolution, heat transfer and deformations, is considered in this paper. Attention focuses on numerical implementation of the thermo-elastic-plastic constitutive equation proposed by Leblond et al. [J. Mech. Phys. Solids 34(4) (1986a) 395; J. Mech. Phys. Solids 34(4) (1986b) 411] in consideration of the transformation plasticity. Based upon the multiplicative decomposition of deformation gradient, hyperelastoplastic formulation is borrowed for efficient numerical implementation, and the algorithmic consistent moduli for elastic-plastic deformations including transformation plasticity are obtained in the closed form. The convergence behavior of the present implementation is demonstrated via a couple of numerical examples. (c) 2004 Published by Elsevier Ltd.