The interaction between thermomechanical parameters and microstructure evolution is so intense that it must be considered during the finite element method (FEM) simulation of the hot plastic working process, for materials that are difficult to deform. Taking the microstructure evolution into account, a novel type of constitutive relationship has been put forward for the IN718 alloy.The microstructure evolution model was first established for the dominant microstructure evolution processes. Then the microstructure evolution models and the method to determine the local flow stress of the corresponding microstructure for current thermomechanical parameters and deformation history were presented. Once the local flow stresses of different structures and their volume contributions were defined, the apparent flow stress of the material could be determined as the weighted sum of the local flow stresses and volume contributions. To validate the proposed method, a thermoviscoplastics and microstructure evolution coupled analysis for a forging process of a critical 1N718 disk forging was performed. The predicting results were in close agreement with the experimental data.