We derive the governing equations for the overall behaviour of linear viscoelastic composites comprising two families of elastic inclusions, subphases and/or fibres, and an incompressible Newtonian fluid interacting with the solid phases at the microscale. We assume that the distance between each of the subphases is very small in comparison to the length of the whole material (the macroscale). We can exploit this sharp scale separation and apply the asymptotic (periodic) homogenization method (AHM) which decouples spatial scales and leads to the derivation of the new homogenised model. It does this via upscaling the fluid-structure interaction problem that arises between the multiple elastic phases and the fluid. As we do not assume that the fluid flow is characterised by a parabolic profile, the new macroscale model, which consists of partial differential equations, is of Kelvin-Voigt viscoelastic type (rather than poroelastic). The novel model has coefficients that encode the properties of the microstructure and are to be computed by solving a single local differential fluid-structure interaction (FSI) problem where the solid and the fluid phases are all present and described by the one problem. The model reduces to the case described by Burridge and Keller (1981) when there is only one elastic phase in contact with the fluid. This model is applicable when the distance between adjacent phases is smaller than the average radius of the fluid flowing in the pores, which can be the case for various highly heterogeneous systems encountered in real-world (e.g., biological, or geological) scenarios of interest.
Keywords: fluid–structure interaction; homogenization; viscoelasticity.