From 9a5ecb7f3f149ea37d2f52178685c78f9342d5a3 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 29 Aug 2012 11:11:33 +0200
Subject: [PATCH] Soft body; new tolerance
---
src/one-body-sample.parset | 12 ++++++------
src/one-body-sample.py | 2 +-
2 files changed, 7 insertions(+), 7 deletions(-)
diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset
index 695ac2ca..8d24747b 100644
--- a/src/one-body-sample.parset
+++ b/src/one-body-sample.parset
@@ -17,7 +17,7 @@ enable_timer = false
refinements = 4
[body]
-E = 1e8
+E = 2e7
# "Most steels and rigid polymers when used within their design limits
# (before yield) exhibit values of about 0.3"
# -- http://en.wikipedia.org/wiki/Poisson's_ratio
@@ -25,12 +25,12 @@ nu = 0.3 # Needs to be strictly below 0.5; the closer we get, the
# tougher the problem. But we do have convergence!
[solver]
-tolerance = 1e-10
+tolerance = 1e-11
[solver.tnnmg]
maxiterations = 1000000
-fixed_point_iterations = 1000
-fixed_point_tolerance = 1e-12
+fixed_point_iterations = 1000000
+fixed_point_tolerance = 1e-11
[solver.tnnmg.linear]
maxiterations = 1000000
@@ -60,7 +60,7 @@ acceptFactor = 1.0
requiredResidual = 1e-12
[boundary.friction]
-normalstress = 5000000
+normalstress = 10000000
# "Most dry materials in combination have friction coefficient values
# between 0.3 and 0.6"
# -- http://en.wikipedia.org/wiki/Friction#Coefficient_of_friction
@@ -84,7 +84,7 @@ model = Dieterich # Ruina
a = 0.0075
b = 0.015
# If L is too large, reaching a steady state takes too long
-L = 1e-8
+L = 1e-3
# The second-order a and b terms in [..] are very important in
# determining whether sliding is stable or unstable and (in appli-
diff --git a/src/one-body-sample.py b/src/one-body-sample.py
index 8f4759d4..cd062556 100644
--- a/src/one-body-sample.py
+++ b/src/one-body-sample.py
@@ -14,7 +14,7 @@ class neumannCondition:
## FIXME: If we change time discretisation, we currently change the problem
class dirichletCondition:
def __call__(self, x):
- return 3e-3 / 16
+ return 2e-3
# return 0
fst = 3e-4
snd = 5e-4
--
GitLab