I am trying to implement fractional derivatives based viscoelastic model in Abaqus using a UMAT subroutine. I am running into convergence issues. The general equation I am using is:
DSTRESS(I) = DSTRESS(I) + DDSDDE(I,J)*DSTRAN(J)-(ah*th(i,j)+ad*td(i,j))*fsumstress(j)+(bh*th(i,j) +...