Firstly, typical gradient-dependent nonlocal inelastic models were briefly gradient-dependent constitutive model for plasticity coupled with isotropic damage was presented in the framework of continuum thermodynamics. Numerical scheme for calculation of Laplacian term of damage field with the numerical results obtained by FEM calculation was proposed. Equations have been presented on the basis of Taylor series for both 2-dimensional and 3-dimensional cases, respectively. Numerical results have indicated the validity of the proposed gradient-dependent model and corresponding numerical scheme.