Thermoelastic vibration is the main cause of brake jitter and noise, including significant fluctuations in temperature, thermal stress, structure size and braking moment. In order to improve the reliability of thermoelastic vibration results, a multi-physical field coupling modeling method is proposed in this paper, which can realize the simultaneous calculation of temperature field, stress field and air flow field. Firstly, the theoretical analysis of the whole braking coupling process and feasibility analysis is carried out. Secondly, ABAQUS solid model and FLUENT fluid model are established respectively, and MPCCI is used to share the parameters of the two models, including temperature, node displacement and heat transfer coefficient. Finally, different paths are constructed to study the fluctuation of thermoelastic parameters in different directions. The research conclusion of thermoelastic vibration can provide important basis for optimization of brake NVH (Noise Vibration and Harshness).