This paper presents a validated three-dimensional thermomechanical model to investigate the effect of the welding sequence on weld induced residual stress distributions in a multipass welded piping branch junction. Three possible symmetrical welding sequences, i.e. one-welder, two-welder and four-welder sequence, which were expected to generate the least adverse stress distributions in actual welding circumstances, were proposed and evaluated. It is shown that the variations of the predicted maximum stresses along the fusion lines are not large for the three proposed welding sequences, although at both the branch and run cross section, the four-welder model predicts the highest peak residual stresses while the one-welder model gives the lowest values. As would be expected, the stress distributions predicted by the four-welder and two-welder models show more symmetrical features than those by the one-welder model. High residual stress is formed in the vicinity of the weld region irrespective of the sequence of the welding.