sigma = 2*mu*symm(gradU)  + lambda*(I*tr(gradU));