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.

In this paper, a meshfree method called adaptive CTM–RPIM is developed to model geotechnical problems with large deformation. The developed adaptive CTM–RPIM is a combination of the Cartesian transformation method (CTM), the radial point interpolation method (RPIM) and the alpha shape method. To reduce the requirement for meshes, the CTM is adopted to transform domain integrals into line integrals, and the RPIM is applied to construct interpolation functions. The alpha shape method, which is capable of capturing severe boundary evolution due to large deformations, is then introduced into the CTM–RPIM to form the adaptive CTM–RPIM. The accuracy of CTM–RPIM is first verified by considering a cantilever beam under small deformation, where the influence of key parameters on the simulation results is explored. Afterward, the ability of the adaptive CTM–RPIM to handle large deformation problems is demonstrated by simulating cantilever beams with large deformations for which analytical solutions are available. Finally, the ability of the proposed method to model the geotechnical large deformations is illustrated from both quasi-static and dynamic aspects, where a slope failure problem and a footing bearing capacity problem are modeled to evaluate the stability of geotechnical structures; and a 2-D soil collapse experiment using small aluminum bars is simulated to show the method capability in describing the soil flows. These benchmark examples demonstrate that the adaptive CTM–RPIM is a numerical method with broad application prospects for modeling large deformation problems in geotechnical engineering.