diff --git a/dune/gfe/cosseratenergystiffness.hh b/dune/gfe/cosseratenergystiffness.hh
index 4d1b8060775f947938f59270dfa156bbbec56947..7c95109e2dd76f98509918c56eaed0d714884019 100644
--- a/dune/gfe/cosseratenergystiffness.hh
+++ b/dune/gfe/cosseratenergystiffness.hh
@@ -365,9 +365,9 @@ public:
     RT curvatureEnergy(const Tensor3<field_type,3,3,3>& DR) const
     {
 #ifdef DONT_USE_CURL
-        return mu_ * std::pow(L_c_ * DR.frobenius_norm2(),q_/2.0);
+        return mu_ * std::pow(L_c_ * L_c_ * DR.frobenius_norm2(),q_/2.0);
 #else
-        return mu_ * std::pow(L_c_ * curl(DR).frobenius_norm2(),q_/2.0);
+        return mu_ * std::pow(L_c_ * L_c_ * curl(DR).frobenius_norm2(),q_/2.0);
 #endif
     }