From ad66c3a5732a9cc896765e6163275fae8ab97615 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 14 Mar 2012 14:33:58 +0100
Subject: [PATCH] New problem (with velocity weakening)
---
src/one-body-sample.parset | 8 ++++----
src/one-body-sample.py | 6 +++---
2 files changed, 7 insertions(+), 7 deletions(-)
diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset
index efe421eb..cf4ef441 100644
--- a/src/one-body-sample.parset
+++ b/src/one-body-sample.parset
@@ -74,16 +74,16 @@ model = Ruina
[boundary.friction.state]
evolve = true
# log(alpha(0))
-initial = 0
+initial = -2.42037
[boundary.friction.ruina]
# "For rocks, typical values of A and B range from 0.005 to 0.015"
# -- Ronaldo I. Borja and Craig D. Foster:
# Continuum mathematical modeling of slip weakening in geological systems
-a = 0.015
-b = 0.005
+a = 0.0075
+b = 0.015
# If L is too large, reaching a steady state takes too long
-L = 1e-4
+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 81d09658..3ab9b22d 100644
--- a/src/one-body-sample.py
+++ b/src/one-body-sample.py
@@ -14,9 +14,9 @@ class neumannCondition:
class dirichletCondition:
def __call__(self, x):
# return 0
- fst = 5e-4
- snd = 2e-4
- trd = 5e-4
+ fst = 3e-4
+ snd = 5e-4
+ trd = 3e-4
if x < 1.0/3:
return fst
elif x < 2.0/3:
--
GitLab