The traditional analysis of mechanical seal temperature field is often assumed that seal end faces are paralleled,without considering the influence of initial stress deformation on temperature field,there is a quite large difference between the temperature in practice and that of seal ring obtained from calculation.The finite element model of sternshaft mechanical seal considering the initial stress and deformation was set up.The temperature distribution of the end faces of mechanical face seals under specific conditions and the change trend of temperature along axial direction in the sealing ring were studied with indirect coupling method and by using Ansys,and the result was compared with that of traditional method.The results show that the presented numerical analysis method considers the effect of initial stress and deformation of end faces，the calculated maximum temperature in axial direction of seal end surfaces occurs in the location close to the actual contact inner diameter after the deformation of end surfaces,but not in the location close to the inner diameter calculated by the traditional method.As compared with the traditional method,the end surface pressure calculated by the presented numerical analysis method is bigger and the surface temperature is higer,especially in high temperature and pressure conditions when using a seal material with small elastic modulus.