The thermomechanical coupling of rocks refers to the interaction between the mechanical and thermodynamic behaviors of rocks induced by temperature changes. The study of this coupling interaction is essential for understanding the mechanical and thermodynamic properties of the surrounding rocks in underground engineering. In this study, an improved temperature-dependent linear parallel bond model is introduced under the framework of a particle flow simulation. A series of numerical thermomechanical coupling tests are then conducted to calibrate the micro-parameters of the proposed model by considering the mechanical behavior of the rock under different thermomechanical loadings. Good agreement between the numerical results and experimental data are obtained, particularly in terms of the compression, tension, and elastic responses of granite. With this improved model, the thermodynamic response and underlying cracking behavior of a deep-buried tunnel under different thermal loading conditions are investigated and discussed in detail.