diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset index f068672174c90c7e3e47f9be052b8e3432caab2d..6abe52255cd0e00f70217d4c35e2cc5c98551df5 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 d9653f3148563ef8f1ec17bcb7c7515c59af004a..4f1beca0afea4711a4581177942a0a08c99041a1 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