The linear isothermo-viscoelastic constitutive equation is established according to the principle of viscoelastic mechanics. Given the boundary conditions of the temperature field, the linear themo-viscoelastic constitutive equation is established acording to the analysis of the thermorheologically simple. The stress analysis model is constructed on the base of some reasonable hypotheses which consider the restraint conditions of mold and the characteristics of injection molding in the post-filling stage. The mathematical model is calculated by the finite difference method. The results can help to predict the warpage of plastic products.