A meso-scale analysis method using the natural element method is proposed for the analysis of material damage of brittle microcracking solids. The microcracking is assumed to occur along Voronoi edges in the Voronoi diagram generated using the nodal points as the generators. The mechanical effect of microcracks is considered by controlling the material constants in the neighborhood of the micorcracks. The macro elastic moduli of isotropic solids containing a number of randomly distributed microcracks are calculated considering the effect of microcrack closure to demonstrate the validity of the proposed method.