From 14eda0a08c78506a46933bccf4988fb4171c6867 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 14 Mar 2012 14:39:38 +0100
Subject: [PATCH] Print the coefficient of friction
---
src/one-body-sample.cc | 19 +++++++++++++++++++
src/one-body-sample.parset | 5 +++--
2 files changed, 22 insertions(+), 2 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 2067be0e..e39e754c 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -282,6 +282,25 @@ int main(int argc, char *argv[]) {
print_evolution<SingletonVectorType, VectorType>(
frictionalNodes, *s4_new, u4, functions.get("neumannCondition"),
run, h * run, octave_writer);
+
+ if (parset.get<bool>("printCoefficient"))
+ for (size_t i = 0; i < frictionalNodes.size(); ++i) {
+ if (!frictionalNodes[i][0])
+ continue;
+
+ double const V = u4_diff[i].two_norm();
+ double const state = (*s4_new)[i];
+ auto const L = parset.get<double>("boundary.friction.ruina.L");
+ auto const a = parset.get<double>("boundary.friction.ruina.a");
+ auto const b = parset.get<double>("boundary.friction.ruina.b");
+ auto const eta = parset.get<double>("boundary.friction.eta");
+ auto const mu = parset.get<double>("boundary.friction.mu");
+
+ std::cout << mu + a *std::log(V * eta) +
+ b * (state - std::log(eta * L)) << std::endl;
+
+ break;
+ }
}
u4 += u4_diff;
diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset
index cf4ef441..1f660220 100644
--- a/src/one-body-sample.parset
+++ b/src/one-body-sample.parset
@@ -2,7 +2,8 @@
timesteps = 600
verbose = false
-printEvolution = true
+printCoefficient = true
+printEvolution = false
printFrictionalBoundary = false
printProgress = false
@@ -13,7 +14,7 @@ enable = false
iterations = 50
[grid]
-refinements = 6
+refinements = 4
[body]
E = 1e8
--
GitLab