thermoelasticity;  contact issues;  complementary equations;  Finite element solution

By defining the sliding yield function and the flow potential function of the contact surface, the complementary equation of the thermal contact boundary is obtained. Based on the principle of virtual work and the contact element technology, the finite element equation of the three-dimensional elastic contact problem under thermal stress is derived, and the contact pseudo-load is introduced into the equation to reflect the influence of the contact state change on the load planting. In the process of shading solving, the strategy of accelerating every two iterations is adopted, which significantly improves the solution efficiency.