Heterogeneous materials are involved in numerous natural phenomena and industrial applications. Therefore, it is crucial to understand how macroscopic mechanical properties of such materials emerge from underlying micro-structure. To address this topic, numerical simulation approaches are valuable tools,especially when discretization scale is smaller than micro-structure scale. In such a way, the influence of structural parameters can be extensively studied.Among available numerical methods, bond-based peridynamics is a non-local approach which relies on breakable elastic bonds between distant material points within a neighborhood of finite size. As peridynamics is a nonlocal approach, the influence of local mesh anisotropy on crack patterns is significantly reduced, compared with local approaches like lattice elements method.In present study, the bond-based peridynamic approach was implemented in 2D and parallelized by message exchange (MPI). The scalability of the calculation tool and the mesh convergence were tested. Finally, three cases are studied:1) the probability of rupture of porous materials; 2) the evolution of the mechanical properties of a cohesive granular medium for different cementitious matrix contents; 3) the evolution of the damage in various phases of a cellular material as a function of cell wall/interface toughness ratio.