This paper addresses a method for numerical simulation in the pressing process of hot glass. Updated Lagrangian finite element formulations are employed for the flow and energy equations to accommodate moving meshes. The model is assumed axi-symmetric and creep flow is assumed due to the high viscosity. Commercial software ANSYS is used to solve the coupled flow and energy equations. Moving contact points as well as free surface during the pressing are effectively calculated and updated by utilizing API functions of CAD software Unigraphics. The mesh distortion problem near the wall is overcome by automatic remeshing, and the temperatures of the new mesh are conveniently interpolated by using a unique function of ANSYS. The developed model is applied to the pressing process of TV glasses. In conclusion, the presented method shows that the pressing process accompanying moving boundary can be simulated by effectively combining general purpose software without resorting to special dedicated codes.