The nonlocal integral approaches are established methods known for their effectiveness in dealing with the challenges of ill-posed behaviour and mesh sensitivity during the strain softening process, especially within classical continuum theories. The present work extends the strain difference-based nonlocal approach to continuum damage mechanics (CDM) theories for predicting the failure behaviour of quasi-brittle materials. The strain difference approach is adopted for two reasons: (i) it maintains physical realism and indirectly incorporates nonlocal effects while computing both the damage variable and the damage energy release rate, and (ii) the constitutive model maintains a symmetry nature in constructing the nonlocal part of the stiffness matrix, reducing the computational cost in strain-softening problems. The state and damage evolution equations for the present nonlocal approach are formulated within the thermodynamic framework. Numerical examples involving 2D and 3D cases are solved: (i) to study the influence of nonlocal parameter ‘α’ on the damage distribution while capturing the post-peak behaviours, (ii) to demonstrate the effectiveness of the present methodology by obtaining mesh-independent results and accurately predicting damage propagation paths, and (iii) to compare the experimental results and numerical structural response for coarse mesh, highlighting both computational efficiency and reduction of mesh-bias dependency. The results obtained by the strain difference-based nonlocal damage model agree well with the experimental results in the literature.