This paper is concerned with the development of a quantitative analysis method that can extract the mixed linear and non-linear dynamical response of a viscoelastic process in composite materials. A tractable set of simultaneous equations with well behaved coefficients can be generated by taking time series moments of a suitably truncated Volterra series expansion. This moment hierarchy is a set of inhomogeneous non-linear integral equations, based on a vector multidimensional convolution form of the Volterra functional series expansion. The hierarchy developed is used to analyse the time dependent viscoelastic properties of fibre reinforced composite materials. Estimates of the temporal response of the measured stress to the measured mechanical force were then used to predict the out of sample stress field values. It is demonstrated that the response functions provided a good representation of the viscoelastic process. Generally speaking, the response functions estimated from the data are used to determine the dynamic and steady state transport coefficients, which can, be used to develop either an empirical field theory of the phenomena or alternatively be used in the design process.