Fisher et al. (2004) and Ghoreishy (2012) conducted viscoelastic modeling using the Prony series, and in this study, modeling was carried out according to Fisher’s research methodology.From the generalized Maxwell model of the linear viscoelastic model, the Prony series can be obtained using Equation (13).The Prony series coefficients were determined in the following manner: First, frequency domain data in the range of 1 to 40 Hz (limited by the range of the DMA instrument) were collected under isothermal conditions at temperatures between -60 and 80 °C (increments of 5 °C). Subsequently, the time-temperature superposition principle was used to construct a master curve in a wide frequency range using the Williams-Landel-Ferry (WLF) Figure (Williams et al., 1955) as shown in Equation (14) and determine the shift factor aT (the horizontal translation factor). Tr is the reference temperature, which can be assigned by an operator, and the material constants were fitted, as shown in Table 6. Finally, these experimental reference curves were fitted with an 8-term Prony series using the linear least squares solver DYNAMFIT (Bradshaw and Brinson 1997), and the parameters are listed in Table 7.