A stochastic meshless method is presented for solving boundary-value problems in damage mechanics under elastic-plastic conditions that involves random material properties. Material is assumed to have an initial damage, which follows a lognormal filed. An isotropic unified damage law, formulated by Lemaitre is used in this study. A meshless formulation based on element free Galerkin method (EFGM) is developed to predict stochastic structural response. A scaled matrix approach is used for applying the essential boundary conditions in EFGM. The proposed method is based on perturbation technique. First order perturbation technique for damage growth in an elastic plastic analysis is formulated. Newton-Raphson algorithm is used to solve for material nonlinearity. A numerical example is solved to study the effect of random initial damage in the structural response and further damage growth. Monte Carlo technique is used to validate the proposed method. © 2008 by ASME.