diff --git a/dune/gfe/cosseratenergystiffness.hh b/dune/gfe/cosseratenergystiffness.hh
index 8753f5eed8c49c571c583833c41d7371f01fca65..3d7c5de88a6c0d3f40156020470021c761aafc27 100644
--- a/dune/gfe/cosseratenergystiffness.hh
+++ b/dune/gfe/cosseratenergystiffness.hh
@@ -323,8 +323,8 @@ energy(const Entity& element,
         
         // Add the local energy density
         if (gridDim==2) {
-            energy += weight * thickness_ * quadraticMembraneEnergy(U);
-            //energy += weight * thickness_ * longQuadraticMembraneEnergy(U);
+            //energy += weight * thickness_ * quadraticMembraneEnergy(U);
+            energy += weight * thickness_ * longQuadraticMembraneEnergy(U);
             energy += weight * thickness_ * curvatureEnergy(DR);
             energy += weight * std::pow(thickness_,3) / 12.0 * bendingEnergy(R,DR);
         } else if (gridDim==3) {