A meso-scale analysis method using the natural element method, which is a kind of meshless method, is proposed for the analysis of material damage of brittle microcracking solids such as ceramic materials, concrete and rocks. 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 microcracks. The macro elastic moduli of anisotropic as well as isotropic solids containing a number of randomly distributed microcracks are calculated in order to demonstrate the validity of the proposed method.