Compressibility and tablet forming ability (compactibility) of bimodal mixtures of differently sized granules formed from microcrystalline cellulose were studied experimentally and numerically with the discrete element method (DEM). Compression data was analysed using the Kawakita equation. A multi-body contact law that accounts for contact dependence resulting from plastic incompressibility/geometric hardening was used in the DEM simulations. The experimental Kawakita a and 1/b parameters both depended non-monotonically on composition (weight fraction of large particles). For the a parameter, this dependence was explained by variations in the porosity of the initial granule beds; for the 1/b parameter, other factors were found to be of importance as well. The numerical results generally compared favourably with the experiments, demonstrating the usefulness of the DEM at high relative densities, provided that a suitable multi-particle contact model is used. For all mixtures, the tensile strength of the formed tablets increased with increasing applied pressure. The tensile strength generally decreased with increasing fraction of large particle, and this decrease was more rapid for large differences in particle size. A possible interpretation of these findings was proposed, in terms of differences in lateral support of small particles in the vicinity of large particles.