diff --git a/dune/gfe/cosseratenergystiffness.hh b/dune/gfe/cosseratenergystiffness.hh
index 5d2f0a979da6c7a006f5b241415dcd937e30c8f1..8753f5eed8c49c571c583833c41d7371f01fca65 100644
--- a/dune/gfe/cosseratenergystiffness.hh
+++ b/dune/gfe/cosseratenergystiffness.hh
@@ -161,7 +161,7 @@ public:
         Dune::FieldMatrix<double,dim-1,dim-1> sym2x2;
         for (int i=0; i<dim-1; i++)
             for (int j=0; j<dim-1; j++)
-                sym2x2[i][j] = 0.5 * (U[i][j] + U[j][i]);
+                sym2x2[i][j] = 0.5 * (U[i][j] + U[j][i]) - (i==j);
 
         result += mu_ * sym2x2.frobenius_norm2();