We present a three-dimensional (3D) numerical model to investigate complex fracture behavior using cohesive elements.An efficient packing algorithm is employed to create the mesoscale model of heterogeneous capsule-based self-healing concrete.Spherical aggregates are used and directly generated from specified size distributions with different volume fractions.Spherical capsules are also used and created based on a particular diameter,and wall thickness.Bilinear traction-separation laws of cohesive elements along the boundaries of the mortar matrix,aggregates,capsules,and their interfaces are pre-inserted to simulate crack initiation and propagation.These pre-inserted cohesive elements are also applied into the initial meshes of solid elements to account for fracture in the mortar matrix.Different realizations are carried out and statistically analyzed.The proposed model provides an effective tool for predicting the complex fracture response of capsule-based self-healing concrete at the meso-scale.