Notes on computational aspects of the fractional-order viscoelastic model

This paper deals with the computational aspect of the investigation of the relaxation properties of viscoelastic materials. The constitutive fractional Zener model is considered under continuous deformation with a jump at the origin. The analytical solution of this equation is obtained by the Laplace transform method. It is derived in a closed form in the terms of the Mittag-Leffler function. The method of numerical evaluation of the Mittag-Leffler function for arbitrary negative arguments which corresponds to physically meaningful interpretation is demonstrated. A numerical example is given to illustrate the effectiveness of this result.