In my code I solve for a 2D displacement field u and then proceed to calculate strain
eps = ( grad(u)+transpose(grad(u)) )
e11 = eps[0, 0]
e22 = eps[1, 1]
eii = eps[1,1]+eps[0,0]
Plots for e11 and e22 are like this:
e11
e22
which look right according to the problem. But if I plot eii then the positive part goes away
somehow,
eii
Am I doing it wrong? Please help.
Regards
Deb.