diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset index 695ac2ca305d70ff9fd1fc518db8f56c443d2595..8d24747b92ab9d56da4e700f5f601f8647117500 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 8f4759d49f5b927db3b535e0058728724c4b0b4d..cd0625563cbd731b6e1b7a60ba68a2b6222d6641 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