The paper presents an extension to the spheropolyhedra method for the simulation of granular materials comprising particles of general shapes with bonding. A bonding, cement, or cohesion model for particles sharing common faces is introduced. The bonding force is elastic and has a strain-based breaking threshold for modelling fracture. An initial study is conducted based on the Brazilian tensile test to check how the parameters of the proposed model affect the principal variables measured in this test. Afterwards, solid cubic blocks are then subjected to a triaxial test to explore the mathematical macroscopic failure model. It is found that the peak strength envelope is the product of the superposition of frictional and fracture failure mec...