In glass compression molding, most current modeling approaches of temperature-dependent viscoelastic behavior of glass materials are restricted to thermo-rheologically simple assumption. This research conducts a detailed study and demonstrates that this assumption, however, is not adequate for glass molding simulations over a wide range of molding temperatures. In this paper, we introduce a new method that eliminates the prerequisite of relaxation functions and shift factors for modeling of the thermo-viscoelastic material behavior. More specifically, the temperature effect is directly incorporated into each parameter of the mechanical model. The mechanical model parameters are derived from creep displacements using uniaxial compression experiments. Validations of the proposed method are conducted for three different glass categories, including borosilicate, aluminosilicate, and chalcogenide glasses. Excellent agreement between the creep experiments and simulation results is found in all glasses over long pressing time up to 900 seconds and a large temperature range that corresponds to the glass viscosity of log (eta) = 9.5 - 6.8 Pas. The method eventually promises an enhancement of the glass molding simulation.