The residual stress distribution obtained by coupling the laser metal deposition residual stress and the milling residual stress is one of the crucial factors affecting the quality and performance of hybrid manufacturing parts. However, it is challenging to explore directly through experiments. A three-dimensional thermal-mechanical fully coupled model for hybrid manufacturing of 316L stainless steel was established, taking into account the multi-physics during laser metal deposition and milling. The accuracy of the hybrid manufacturing model is validated from multiple dimensions using temperature, deformation, and milling force. A height deviation normalization method is proposed to eliminate the influence of printing height deviation caused by randomly fluctuated process parameters on the measuring deformation of the part. The minimum root-mean-square relative error of model temperature is 4 %, the minimum error of printing height deviation is 8 %, and the minimum relative error of milling force in the feed direction is 26 %. The simulated milling force of hybrid manufacturing shows that the factors affecting the milling force rank as the axial depth of cut, the radial depth of cut, milling speed, and spindle speed. The simulated residual stress distribution zone along the depth direction of hybrid manufacturing divides into the work hardening layer, the elastic recovery zone, and the elastic stability zone. The hybrid manufacturing residual stress equals the milling residual stress in the work hardening layer, superimposes linearly with the product of the milling residual stress and depth in the elastic recovery zone, and equals the residual stress of laser metal deposition in the elastic stability zone.