This paper introduces a three-dimensional (3D) semi-explicit probabilistic numerical model for simulating crack propagation within the framework of the finite element method. The model specifically addresses macrocrack propagation using linear volume elements. The criteria governing the macrocrack propagation is based on the softening behavior observed in concrete under uniaxial tension. This softening behavior corresponds to a dissipated cracking energy that is equal to the mode I critical fracture energy (GIC) used in the Linear Elastic Fracture Mechanics theory (LEFM). The probabilistic nature of this model revolves around the random distribution of two mechanical properties: tensile strength (ft) and fracture energy, which varies based on the volume of finite elements. The scattering of the fracture energy increases as the volume of finite elements decreases in order to consider the strong heterogeneity of the material. This work primarily aims to estimate the relationship between the standard deviation of GIC and the volume of finite elements through the development of the numerical model. For this purpose, an inverse analysis is conducted based on a fracture mechanical test simulation. This test involves macrocrack propagation in a large Double Cantilever Beam (DCB) specimen with a crack length exceeding two meters. The proposed inverse analysis procedure yields highly significant results, indicating that the numerical model effectively evaluates both crack length and crack opening during propagation.