The mechanical properties of rock can be weakened under the influence of high temperatures. To describe the mechanical behavior of rock under the action of high temperature more accurately, based on the Zienkiewicz–Pande yield criterion, the damage variable Dc which accounts for the coupling between high temperature and load is introduced. According to plastic potential theory and plastic flow law, the iterative incremental method for a high-temperature and load-coupled damage constitutive model in Flac3D is deduced in detail and compiled into the corresponding dynamic link library file (.dll file). By modifying the shape function to degenerate into the Mohr–Coulomb constitutive model, an elastic–plastic analysis of an ideal circular tunnel is performed, and a comparison is made between calculation results obtained from the built-in Mohr–Coulomb constitutive model in Flac3D, proving the correctness of the secondary development program. Finally, numerical simulations are conducted to study the effects of high-temperature damage using rock uniaxial compression tests, and the model’s validity is established by comparing it with previous experimental results.