Abstract
Viscoelastic, or composite materials that are hereditary deformable, have been characterised by exponential and weakly singular kernels in a hereditary equation. An exponential kernel is easy to be numerically implemented, but does not well describe complex vibratory behaviour of a hereditary deformable system. On the other hand, a weakly singular kernel is known to describe the complex vibratory behaviour, but is nontrivial to be numerically implemented. This study presents a numerical formulation for solving a hereditary equation with a weakly singular kernel. Recursive algebraic equations, which are numerically solvable, are formulated by using the Galerkin method enhanced by a numerical integration and elimination of weak singularity. Numerical experiments showed that the present approach with a weakly singular kernel is well fitted into a realistic vibratory behaviour of a hereditary deformable system under dynamic loads, as compared to the same approach with an exponential kernel.