Abstract: Firstly, typical gradient-dependent nonlocal inelastic models were briefly reviewed. Secondly, based on the principle of ‘gradient-dependent energy dissipation', a 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.