Soils exhibit non-linear stress–strain behaviour, even at relatively low strain levels. Existing Winkler models for suction caisson foundations cannot capture this small-strain, non-linear soil behaviour. To address this issue, this paper describes a new non-linear elastic Winkler model for the uniaxial loading of suction caissons. The soil reaction curves employed in the model are formulated as scaled versions of the soil response as observed in standard laboratory tests (e.g. triaxial or simple shear tests). The scaling relationships needed to map the observed soil-element behaviour onto the soil reaction curves employed in the Winkler model are determined from an extensive numerical study employing three-dimensional finite-element analysis. Key features of the proposed Winkler model include: computational efficiency, wide applicability (it can be used for caisson design in clay, silt or sand) and design convenience (the required soil reaction curves can be determined straightforwardly from standard laboratory test results). The proposed model is suitable for small and intermediate caisson displacements (corresponding to fatigue and serviceability limit state conditions) but it is not applicable to ultimate limit state analyses.