In this paper a new boundary-only based damage modeling is presented. Despite the damage occurred inside the problem domain, the solution procedure requires boundary-only discretization of the problem. A few additional points are located at damage places. The Eshelby equivalent inclusion theory is coupled with the direct boundary integral equation method to predict the damage patterns. The relevantment like damage stiffness matrix which is obtained directly in a condensed form at degrees of freedom along the problem boundary. A nonlinear incremental–iterative solution procedure is implemented. Several numerical examples are presented to demonstrate the validity of the proposed technique.