A stochastic meshless method is presented for solving boundary-value problems in damage mechanics under linear elastic conditions that involves random material properties. The material properties are formulated as a function of isotropic damage and damage is modeled as a homogeneous random field. A meshless formulation is developed to predict stochastic structural response. Unlike the finite element method, the meshless method requires no structured mesh, since only a scattered set of nodal points is required in the domain of interest. There is no need for fixed connectivities between nodes. In conjunction with the meshless equations, classical perturbation expansions are derived to predict second-moment characteristics of response. Numerical examples based on one- and two-dimensional problems are presented to examine the accuracy and convergence of the stochastic meshless method. A good agreement is obtained between the results of the proposed method and Monte Carlo simulation. Since mesh generation of complex structures can be a far more time-consuming and costly effort than the solution of a discrete set of equations, the meshless method provides an attractive alternative to finite element method for solving stochastic mechanics problems.