This work is devoted to the fracture processes numerical modeling of elastic-brittle bodies taking into account the statistical distribution of subregions strength values and stress concentration. The novel formulation of the boundary value problem and its solution algorithm are developed. The loading diagrams obtained in computational simulations, the corresponding growth curves of the destroyed elements relative number and the damaging process kinetics are analyzed. The presence of the postcritical deformation stage at the macro-level is noted. The influence of the strength properties distribution range and the depth of the concentrator on the maximum load value and the damage evolution is determined. The significant influence of the finite elements’ properties distribution on the fracture processes modeling results is concluded. © 2023, Frat. Integrita Strutr. All rights reserved.