From 29890363c84fab7f828981a6fa966dcf438777da Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Tue, 22 Nov 2011 14:08:45 +0100
Subject: [PATCH] Bug fix: s -> eta*s
---
dune/tectonic/nicefunction.hh | 12 ++++++++++--
1 file changed, 10 insertions(+), 2 deletions(-)
diff --git a/dune/tectonic/nicefunction.hh b/dune/tectonic/nicefunction.hh
index 81625152..fdcdc0e4 100644
--- a/dune/tectonic/nicefunction.hh
+++ b/dune/tectonic/nicefunction.hh
@@ -51,18 +51,26 @@ class RuinaFunction : public NiceFunction {
y *= coefficient;
}
+ /*
+ (leaving some terms aside): with s > rho
+
+ 1/eta d/dx [ a * (s log s - s) + mu s ] where s = eta x
+ = 1/eta [ a * (log (eta x) * eta) + eta mu ]
+ = a * log(eta x) + mu
+ */
double virtual leftDifferential(double s) const {
if (eta * s <= rho)
return 0;
- return coefficient * normalStress * (a * std::log(s) + mu);
+ return coefficient * normalStress * (a * std::log(eta * s) + mu);
}
+ /* see above */
double virtual rightDifferential(double s) const {
if (eta * s <= rho)
return 0;
- return coefficient * normalStress * (a * std::log(s) + mu);
+ return coefficient * normalStress * (a * std::log(eta * s) + mu);
}
private:
--
GitLab