The welding residual stress of beamto-column joints was simulated by finite element software in order to study the distribution of welding residual stress in beamto-column joints，and the experimental verification was carried out． Firstly，a stochastic thermodynamic model for the whole welding process of beamto-column joints was established to analyze the welding temperature field of beamto-column joints． Secondly，the distribution of welding residual stress of beamto-column joints was analyzed by indirect thermalstructural coupled analysis． Finally，the welding residual stress of beamto-column joints was measured by the blind hole method and compared with the simulation． The results show that due to the far distance between the upper and lower flange welds and the existence of process holes， the interaction between the welds is very small，and the welding residual stresses of the two welds are basically the same; owing to the close distance between the web welds and the great interaction between the two welds， the distribution of the residual stresses of the two welds is inconsistent; in the weld toe area of beamto-column joints， the distribution of residual stresses in the upper and lower flanges of beams is basically the same， the equivalent residual stresses are“V”type distribution， the horizontal， longitudinal and normal residual stresses are basically“M”type distribution，while the stress distributions of the two welds in the web of beams are obviously different，and the reason is that the two welds of web have mutual influence during welding．