How to implement gradient-enhanced model with damage plasticity?

Dear all,

I am doing works on the numerical implementation of gradient enhanced model with damaged plasiticity. In my model, the nonlocal formulation is carried on the internal variables of concrete, which is denoted as κ={κt, κc}. In the local CDM framework, the yield surface F(σ_,κ), the plastic potential φ(σ_,κ) and the damage evolution law d(κ) are all functions of the internal variables. But I wonder whether I should use the local κ or the nonlocal κ_ to describe F and φ in the plastic part? If I use the nonlocal κ_ for F and φ in the plastic corrector, the value of κ_ that obtained from interpolation of the nodal internal variables would be changed after plastic corrector. However, in my mind, within each iteration step, the value of κ_ should remain the same since it's just a parameter obtained from interpolation... But if I use local κ for F and φ in the plastic corrector, and then use nonlocal κ_ in the damage part, I cannot get the correct results since the damage computed from κ_ always seems to be relatively small... I hope you could understand my meaning and give some advices.  Thank you and best regards,Cheng Ma