In this paper, based on Simufact Welding finite element analysis software, a numerical simulation of the temperature and residual stress distribution of the weldolet-header multi-layer multi-pass welding process is carried out, and the simulation results are verified through experiments. The experimental results are in good agreement with the numerical simulation results, which proves the validity of the numerical simulation results. Through the results of the numerical simulation, the influence of the welding sequence and interlayer temperature on the temperature and residual stress distribution at different locations of the saddle-shaped weld was studied. The results show that the temperature and residual stress distribution on the header and weldolet are asymmetric, and the high-stress area of the saddle-shaped welded joint always appears at the saddle shoulder or saddle belly position. When the interlayer temperature is 300 °C, the peak residual stress reaches a minimum of 428.35 MPa. Adjusting the welding sequence can change the distribution trend of residual stress. There is no high-stress area on the first welding side of the two-stage welding path-2. The peak values of residual stresses for continuous welding path-1 and two-stage welding path-2 are 428.35 MPa and 434.01 MPa, respectively, which are very close to each other.