It would be convenient to have a model, albeit approximate, of particle-laden materials (suspensions) that would not need large amounts of computing and/or experimentation to implement for design purposes. There are now adequate models of the pure matrix fluid behaviour, but there are no such models for suspensions with large particles (non-colloidal suspensions). One of the obstacles has been the single-minded devotion to shearing motions of suspensions; experience with the matrix modelling has shown that it is not possible to formulate widely usable models if only shear is considered. Here some new results of axially symmetric elongational tests on suspensions are compared with shearing data. Some suggestions for modelling these and other observations based on using strain rate and strain in a modified Reiner-Rivlin constitutive equation are presented. The model generally works quite well, but it does not predict the positive storage modulus seen in small and medium amplitude oscillatory shear flows.