Elimination of memory from the equations of motion of hereditary   viscoelasticity for increased efficiency of numerical integration
  A method of eliminating the memory from the equations of motion of linear viscoelasticity is presented. Replacing the unbounded memory by a quadrature over a finite or semi-finite interval leads to considerable reduction of computational effort and storage. The method applies to viscoelastic media with separable completely monotonic relaxation moduli with an explicitly known retardation spectrum. In the seismological Strick-Mainardi model the quadrature is a Gauss-Jacobi quaddrature. The relation to fractional-order viscoelasticity is shown 
Andrzej Hanyga
astro-ph.IM astro-ph.EP math.NA