In linear viscoelasticity, an important consideration is the Boltzmann causal integral equation sigma(t) = integral(-infinity)(t) G(t - tau)(gamma) over dot (tau) d tau, which defines how the stress sigma(t) at time t depends on the earlier history of the strain rate (gamma) over dot (t) via the relaxation modulus (kernel) G(t). The Boltzmann model of linear viscoelasticity is an appropriate model for materials that simultaneously exhibit viscous and elastic behaviour. In the formulation of a Boltzmann model, the key consideration is the choice of the relaxation modulus. In a wide variety of applications, including the modelling of the glassy state of dense matter, polymer dynamics, and bone and muscle rheology, the Kohlrausch function has proved to be more appropriate in modelling the associated relaxation and decay processes than the standard exponential function. It is therefore an appropriate practical choice for the relaxation modulus. In this paper, the relaxation modulus is modelled as the sum of Kohlrausch functions. It is established theoretically how the unknown parameters in such models can be recovered using moments of the measured stress, and of the corresponding known applied strain rates. For a single Kohlrausch relaxation modulus, the proposed methodology is quite simple to implement. Three moments of both the measured stress and the applied strain-rate are evaluated numerically. A simple algorithm then determines the Kohlrausch parameters as defined by these three moments. Simulations are performed to assess and confirm the utility of the numerical performance of these two steps. Crown Copyright (c) 2004 Published by Elsevier B.V. All rights reserved.