Aiming at modelling the mechanical behaviour of soil-rock mixtures accurately and efficiently, a coupled MPM-DEM formulation combining the material point method (MPM) and the discrete element method (DEM) is proposed. It is solved concurrently via the contact force linking the two individual methods. Specifically, the soil is modelled with MPM as continuums to avoid handling the contacts between fine particles. The rocks are modelled by DEM to capture the contact characteristics of rocks. This method is validated with ball impacting and block sliding tests first for the contact between material points and DEM particles. Its capability in describing the mechanics of soil-rock mixtures is thereafter proved by comparing the simulation results with pure DEM simulations of binary mixtures and laboratory tests of soil-rock mixtures. It is demonstrated that MPM-DEM can reproduce the stress–strain response of soil-rock mixtures and capture the influence of rock contents and rock sizes. In addition, a coarse-graining modelling scheme is implemented, i.e., representing the soil particles with fewer material points, which significantly increases the efficiency compared with pure DEM. Our proposed method provides a novel way to model soil-rock mixtures with reasonable computational efforts, which sheds light on simulating large-scale soil-rock mixtures in nature or engineering.