Mechanical behavior of flexible risers can be challenging due largely to its complex design generating strong nonlinear problems. Nonlinear material properties, as one of them, from polymer layers dominate the overall viscoelastic responses of flexible risers which may play an inevitable role on the global analysis in deepwater application. An alternative to predict the viscoelastic behavior comprising of the time domain and the frequency domain has been proposed recently by the authors (Liu and Vaz, 2016). Given the fact that polymeric material properties are temperature-dependent and that the temperature profiles in flexile risers vary continuously in both axial and radial direction, the temperature of the internal hydrocarbons must affect the viscoelastic responses. However, such phenomenon dose not draw much attention in previous studies. This paper presents an improved model for overcoming some drawbacks in the proposed model involving assumption of steady temperature distribution in polymer layer and no gap appearance between the adjacent layers. The computing method of model is developed by using a step by step test approach. Consequently, some important parameters like equivalent axial stiffness, contact pressure or gap between the near layers, and force-deformation relationship can be observed. Parametric studies are conducted on the axisymmetric viscoelastic behavior of flexible risers to study the role of input temperature and loading frequency. Results show that equivalent axial stiffness given by the improved model is smaller than before. It can also be found that the gap between metal layer and polymer layer appear easily and increases as time goes on.