This paper explores the computational modeling of nonlocal strain, damage, and fracture in concrete, considering the isolated contribution of two random, spatially variable properties related to the fracture process: Young’s modulus (E) and tensile strength (ft). Applying a continuum damage model, heterogeneous specimens of concrete with random and spatially varying E or ft were found to produce substantial differences when contrasted with traditional homogeneous (non-random) specimens. These differences include variable and uncertain strain and damage, wandering of the failure paths, and differing (sometimes lower) peak forces, i.e. increased probabilities of failure in the heterogeneous specimens. It is found that ft variability contributes more (from 1.7 to up to 4 times more, depending on the parameter) to the overall performance variability of the concrete than E variability, which has a comparatively lower contribution. Performance is evaluated using (1) force-displacement response, (2) individual, average, and standard deviation maps of non-local strain and damage, (3) fracture paths and strain and damage values along the fractures. The modeling methodology is illustrated for two specimen geometries: a square plate with a circular hole, and an L-shaped plate. The computational results correlate well with reported experimental data of fracture in concrete specimens.