A constitutive model, accounting for multiphase and multiscale coupling, is proposed for the water retention domain and the stress–strain response of compacted clayey soils. The model is based on a conceptual interpretation of the microfabric evolution of compacted soils along generalised hydromechanical paths, detected by means of mercury intrusion porosimeter tests. Multiphase coupling is provided by the mutual interaction between the mechanical and the hydraulic states. Multiscale coupling is introduced by a measure of the size of the aggregates, which influences both the retention and the stress–strain response, in the phenomenological constitutive equations. Model capabilities are verified by comparison with relevant experimental data from laboratory tests on compacted Boom clay and other selected experimental data on different compacted clayey soils.