The interaction model between the shield and the soil was simplified into the elastic half-space hole problem in order to analyze the contact stress between shield shell and soil during shield tunneling. The shield tunneling was simplified into radial, vertical and horizontal displacement boundary modes. Then a method of the additional contact stress between the shield shell and soil was proposed by using the complex variable function method, and the method was verified by three-dimensional numerical simulation. The sensitivity analysis of relevant parameters was conducted by using the method. Results show that the additional contact stress appears a number of gradually increasing stress peaks with the increase of the relative displacement of the machine and soil, and the corresponding peak stresses are approximately the same when the displacements in different modes are the same. When the displacement boundary condition is the same, the smaller Poisson's ratio and the larger elastic modulus of soil are, the greater the additional contact stress at the extreme point is, while the burial depth has little effect on the mechanical - soil additional contact stress. The distribution function of additional stress coefficient around the shield shell was defined, and the balance equations of force and moment increment of the shield were established. The calculations show that the horizontal posture adjustment load, vertical posture adjustment load, yaw moment of the shield machine are proportional to the change of attitude angle, and the longitudinal posture adjustment load is less affected by the change of attitude angle.