During the service of rolling bearings,the rollingslip phenomenon occurs and residual stress is generated within the material,which makes the Hertz contact stress can not reflect the real stress state,and causes the Hertz contact stress distribution to deviate from experimental results.It is necessary to calculate the internal stress under rollingslip contact for obtaining the real stress field.A fast and simple calculation method of the internal stress under rollingslip contact was presented to replace the the timeconsuming finite element method.Based on the existing formulas,the stress field in the material under rollingslip contact was calculated by Matlab programming.The variation was compared of the internal stress field in twodimensional and threedimensional rollingslip contact under different friction coefficients.The results show that with the increasing of friction coefficient,the maximum principal shear stress becomes larger,the position is closer to the surface and the angle with rolling direction is smaller.The tangential friction force changes the magnitude and position of the maximum orthogonal shear stress on both sides of the contact point.As friction coefficient increases,the maximum orthogaonal shear stress on one side of the contact point becomes larger and the position becomes closer to the surface,while on the other side,the change shows an opposite trend.The method of internal stress calculation presented is convenient,the calculation results are analytical solutions,which is convenient to combine with other stress fields to obtain the real stress field inside the material.Through case analysis,it is revealed that the real stress field can explain the experimental phenomena better,which is significant to the study of microstructural evolution in bearing materials.