From 52fff792b16dfadca84a5b94ed38e0ed1a4b6fd2 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 3 Jul 2013 16:13:17 +0000
Subject: [PATCH] Bugfix in nonquadraticMembraneEnergyGradient
[[Imported from SVN: r9290]]
---
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 29503a15..4a216dce 100644
--- a/dune/gfe/cosseratenergystiffness.hh
+++ b/dune/gfe/cosseratenergystiffness.hh
@@ -646,8 +646,8 @@ nonquadraticMembraneEnergyGradient(typename TargetSpace::EmbeddedTangentVector&
for (size_t v_i=0; v_i<7; v_i++)
for (size_t i=0; i<3; i++)
for (size_t j=0; j<3; j++)
- embeddedLocalGradient[v_i] += 2 * (detU - 1 - (1.0/detU -1)/(detU*detU)) * adjU[j][i] * dU_dv[i][j][v_i];
-
+ embeddedLocalGradient[v_i] += mu_*lambda_/(2*mu_+lambda_) * (detU - 1 - (1.0/detU -1)/(detU*detU)) * adjU[j][i] * dU_dv[i][j][v_i];
+
}
--
GitLab