Numerical Simulation of Crack Propagation in Concrete with Prefabricated Array Fractures Based on the Discrete Element Method
Concrete readily develops cracks under service loads, which poses severe risks to the overall safety of engineering structures. In this work, the discrete element method (DEM) integrated with PFC2D 5.0 numerical software is adopted to construct a mesoscale concrete numerical model containing pre-existing internal fractures, and uniaxial compressive loading simulations are subsequently carried out. Unlike previous studies that predominantly examined isolated fracture parameters, this work systematically investigates the coupled effects of fracture inclination angle, length, and quantity on crack propagation mechanisms at the mesoscale, and for the first time establishes a quantitative relationship between microcrack spatial distribution patterns and macroscopic mechanical degradation. Parametric analyses are performed to quantify the influences of fracture geometric characteristics, including fracture inclination angle (30°, 45°, 60°), fracture length (short, long and extra-long), fracture quantity (4, 8 and 16), as well as the comparison between intact and fractured concrete specimens. The fracture quantities of 4, 8, and 16 are selected to represent low, medium, and high levels of initial defect density within the concrete matrix, corresponding to approximately 1%, 2%, and 4% of the total specimen area, respectively, thereby enabling a systematic investigation into the progressive deterioration of mechanical performance with increasing internal damage severity. The whole evolution process of crack initiation, crack propagation and ultimate failure patterns of concrete is systematically explored. Numerical results reveal that specimens with larger fracture angles exhibit higher compressive strength yet generate abundant newly formed microcracks, whereas low-angle prefabricated fractures are prone to triggering abrupt brittle failure. Specimens embedded with shorter fractures achieve superior mechanical strength and develop denser, more intensive microcrack distributions; in contrast, long pre-existing fractures drastically degrade compressive strength while limiting the generation of secondary cracks. Reducing the number of internal defects simultaneously improves compressive strength and expands the coverage range of the induced fracture network. Specimens with 16 prefabricated fractures deliver the weakest mechanical performance, owing to the excessively high initial defect density inside the matrix. In comparison with fractured samples, intact concrete without pre-set fractures achieves better comprehensive performance in terms of compressive strength, deformation compatibility and uniform microcrack development. A core conclusion drawn from this study is that the total quantity of microcracks cannot serve as a direct indicator to evaluate the damage degradation degree of concrete; instead, the spatial distribution pattern of microcracks dominates the deterioration level. Evenly scattered microcrack populations maintain relatively high residual strength, whereas the concentrated coalescence of microcracks into continuous penetrating macrocracks leads to an abrupt decline in structural load-carrying capacity. The findings of this research can provide theoretical references for stability evaluation and safety diagnosis of defective concrete structures in practical engineering.