We investigate compressive failure of heterogeneous materials on the basis of a continuous progressive damage model. The model explicitely accounts for tensile and shear local damage and reproduces the main features of compressive failure of brittle materials like rocks or ice. We show that the size distribution of damage-clusters, as well as the evolution of an order parameter, the size of the largest damage-cluster, argue for a critical interpretation of fracture. The compressive failure strength follows a normal distribution with a very small size effect on the mean strength, in good agreement with experiments.Understanding how materials break is a fundamental problem that has both theoretical and practical relevance. The topic has received considerable renewed attention during the last few decades because of the limitations of the classical Griffith theory for heterogeneous media [1][2][3]. The practical applications are numerous, from material and structural design, to the important problem of size effects on strength (e.g. [4]). The two key components that make material failure challenging to understand are long-range interactions and material disorder.Traditionally, Weibull and Gumbel distributions associated with the weakest-link approach have been widely used to describe the strength of brittle materials. These distributions naturally arise from extreme-value statistics of initial defect distributions based on the assumptions that [5] (i) defects do not interact with one another, (ii) failure of the whole system is dictated by the activation of the largest flaw (the weakest-link hypothesis), and finally (iii) the material strength can be related to the critical defect size. These assumptions are reasonable for materials with relatively weak disorder loaded under tension, but do not hold for heterogeneous materials with broad distribution of initial disorder or for loading conditions stabilizing crack propagation, such as compression. In these cases there is experimental evidence that a considerable amount of damage can be sustained before failure [6]. Under these conditions, failure is the culmination of a complex process involving the nucleation, propagation, interaction and coalescence of many microcracks [7]. Stress states observed under various natural conditions, ranging from soil and rock mechanics to earthquake physics, suggest the importance of compressive failure.A cornerstone for the understanding of breakdown of disordered media has been lattice models of fracture in * lucas.girard@geo.uzh.ch which networks with prescribed bond failure thresholds are subject to increasing external loads [2]. In these models, failure is described on a qualitative level as the interplay between disorder and elasticity. When strong disorder is considered, these models suggest that fracture strength does indeed not follow a Weibull or a Gumbel distribution but a log-normal distribution [8,9]. Similar strength distributions have been obtained for different model types (fuses and springs), in 2D and 3D, suggesting ...