It is necessary to solve the thermo-mechanical coupled contact problem in structural mechanical analysis, such as the dam structural analysis. Because of the complexity of coupled models, it is very difficult to solve the discretized system in structural mechanical analysis. At the same time, for real applications such as the dam structural analysis, the simulation domain has complex structures, which result in a number of mesh elements more than 10 8 for high resolution simulations. Therefore, the scale of the discretized system is very large. In this paper, the discretized themo-mechanical coupling contact system on hundreds million unstructured meshes is parallel solved by the Newton-Krylov method, in which the efficiency of the Krylov methods is strongly dependent on the preconditioning. An efficient preconditioning method is constructed for the themo-mechanical coupling contact problem. Three steps are used to construct the preconditioner. Firstly, for the mechanical problem, the mechanical effect is analyzed for the dam structural analysis, and a preconditioner is constructed for the elasticity