From 49638df76574e1d48ab75cd70b555340b9353664 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 28 Aug 2014 08:06:45 +0000
Subject: [PATCH] [bugfix] Add missing skew part in finite-strain membrane
energy
This bug has now influence as long as mu_c is zero, but for
non-zero mu_c there was one term missing.
[[Imported from SVN: r9865]]
---
dune/gfe/cosseratenergystiffness.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dune/gfe/cosseratenergystiffness.hh b/dune/gfe/cosseratenergystiffness.hh
index eeba17e5..a7b31d83 100644
--- a/dune/gfe/cosseratenergystiffness.hh
+++ b/dune/gfe/cosseratenergystiffness.hh
@@ -204,7 +204,7 @@ public:
RT detU = U.determinant();
- return mu_ * sym(UMinus1).frobenius_norm2()
+ return mu_ * sym(UMinus1).frobenius_norm2() + mu_c_ * skew(UMinus1).frobenius_norm2()
+ (mu_*lambda_)/(2*mu_ + lambda_) * 0.5 * ((detU-1)*(detU-1) + (1.0/detU -1)*(1.0/detU -1));
}
--
GitLab