From 5dbb98a32e6f90233eb20280fe80367615ca94bb Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Sun, 26 Aug 2012 20:23:48 +0200
Subject: [PATCH] Solve a different problem
---
src/one-body-sample.parset | 6 +++---
src/one-body-sample.py | 1 +
2 files changed, 4 insertions(+), 3 deletions(-)
diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset
index f0686721..6abe5225 100644
--- a/src/one-body-sample.parset
+++ b/src/one-body-sample.parset
@@ -1,5 +1,5 @@
# -*- mode:conf -*-
-timesteps = 500
+timesteps = 1000
verbose = false
printCoefficient = false
@@ -60,7 +60,7 @@ acceptFactor = 1.0
requiredResidual = 1e-12
[boundary.friction]
-normalstress = 0.1 # laursen depends a lot more on this
+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-3 # 1e-3 for the Dirichlet case, 1e-5 - 1e-7 for the Neumann case
+L = 1e-8 # 1e-3 for the Dirichlet case, 1e-5 - 1e-7 for the Neumann case
# 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 d9653f31..4f1beca0 100644
--- a/src/one-body-sample.py
+++ b/src/one-body-sample.py
@@ -13,6 +13,7 @@ class neumannCondition:
class dirichletCondition:
def __call__(self, x):
+ return 3e-3
# return 0
fst = 3e-4
snd = 5e-4
--
GitLab