A mathematical model was developed for the characterization of quasistatic nonlinear viscoelastic behaviour of large arteries with activated smooth muscle. Arteries were considered to belong to the class of viscoelastic materials with fading memory and the first order term of the integral polynomial constitutive equations of Pipkin and Rogers satisfying the nonlinear superposition principle was modified to handle responses to increasing and decreasing loads independently. The two creep functions contained by the obtained one dimensional constitutive equation were determined experimentally from the series of creep and recovery tests of increasing amplitude performed on isolated canine iliac arteries following the activation of the vascular smooth muscle by normal dose of norepinephrine in vitro. Utilizing the constitutive equation of the arterial wall and the tabulated values of creep functions successive stress-strain hysteresis loops of various constant stress rates were simulated by digital computer. The computed hysteresis loops demonstrated the main characteristics, such as the weak and asymmetric rate-sensitivity of the experimentally observed hysteresis of arteries qualitatively well, thus allowing certain conclusions on the mechanism of quasistatic viscoelastic behaviour of vascular smooth muscle.