From 0fd0cfca4cc79356f801e01222c434eeb59a7844 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Mon, 19 Dec 2011 16:25:34 +0000
Subject: [PATCH] Use the long form for the quadratic energy stiffness. I
tested it: when kappa=1 is does the same as the short form
[[Imported from SVN: r8308]]
---
dune/gfe/cosseratenergystiffness.hh | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/dune/gfe/cosseratenergystiffness.hh b/dune/gfe/cosseratenergystiffness.hh
index 8753f5ee..3d7c5de8 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) {
--
GitLab