i found a way actually If you are really interested please write your email in a pm to me.

I am trying to model a self similar analysis procedure and in my case I have to start calculating from the most left element and calculate the integra...

I found an interesting error and fixed it now the algorithm works okey but i found out that my DDSDDE values are negative in the phase transition regi...

So I got the DDSDDE before and after here are they why are they negative DDSDDE1 * 82788.4615384615 * * * *35480.7692307692 * * * *35480.7692307692 ...

I did it as you said and this time abaqus has some errors like , ***WARNING: THE SYSTEM MATRIX HAS 16 NEGATIVE EIGENVALUES. IN AN EIGENVALU...

thnx again i am using return mapping algorithm for my yield function and then i get updated stress and jacobian values regarding to new phase transfor...