Abstract-Based on the concepts of statistical mesoscopic damage mechanics, the rupture of a heterogeneous medium is investigated in terms of numerical simulations of a network model, subjected to simple shear loading. The heterogeneities are simulated by varying the sizes and fracture strains of the elements of the network. Progressive damage is governed by a damage field equation and a dynamic function of damage (DFD). From the damage field equation, a criterion for damage localization can be derived, and the DFD can be extracted from the simulations of the network. Importantly, the DFD intrinsically governs the damage localization. Both stress-free and periodic boundary conditions for the network are examined. It is found that damage localization may be the underlying mechanism of eventual rupture and thus could be used as a possible precursor of earthquake rupture.