A numerical approach to interparticle breakage is applied using the rock failure process analysis code, RFPA2D. A 2D particle assembly in a container is first numerically simulated to obtain the fringe patterns of stress fields that resemble the photoelastic test. Then, in addition, the interparticle breakage of the particle assembly in a chamber is conducted. The chamber consists of a steel container and a steel platen for transferring the load, and contains 15 particles of arbitrary sizes and irregular shapes. A plane strain condition is assumed. The particle bed is loaded under form conditions, in which the size reduction and the applied force are a function of the displacement. The numerical results indicate that, during the crushing process, three principal regimes appear: (i) the elastic deformation regime, where each particle deforms elastically; (ii) the fragmentation regime, where the particle assembly is crushed in a particle-by-particle fashion; and (iii) the assembly hardening regime, where the densified assembly recovers a significant stiffness. The dominant mode of failure is at first splitting, which is more or less parallel to the loading direction, and then progressive crushing, which mainly depends on the confinement from the chamber walls. The analysis of the load–displacement curves of the assembly obtained from the simulations reveals a high undulating load plateau, which suggests a macro-ductile behaviour.