Numerical method for thermoelastic contact of coated materials with thermal imperfections based on a modified conjugate gradient method

Thermal imperfections may arise at the coating-substrate interface during manufacturing processes such as fusion welding or additive layer deposition. The evaluation of thermoelastic contact behavior in coating-substrate systems is essential for reliability assessment. Combining the Discrete Convolution-Fast Fourier Transform (DC-FFT) algorithm, a modified Conjugate Gradient Method (CGM) is developed to establish a thermoelastic contact model for an elastic sphere sliding on coated materials with low-conductivity (LC) and high-conductivity (HC) imperfection interfaces. The modified CGM demonstrates good convergence and computational efficiency in handling the heat partition under various magnitudes of thermal imperfections in contrast with conventional CGM. The contact model is also validated through comprehensive analysis of thermoelastic responses and heat partition behavior. Numerical results based on the present contact model reveal that thermal imperfections significantly influence the profile of temperature, pressure and stress components distribution. Specifically, LC imperfection causes the obvious jumping behaviors of temperature and in-plane stress at the coating-substrate interface, and HC imperfection leads to significant decay rate of temperature and in-plane stresses within coating layer. Furthermore, the existence of imperfections will affect the sensitivity of temperature rise to system parameters compared with thermally perfect condition.